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FOREWORD 


This  report  was  prepared  by  The  Boeing  Company,  Vertol 
Division  of  Philadelphia,  Pennsylvania,  for  the  Aerospace 
Dynamics  Branch,  Vehicle  Dynamics  Division,  Air  Force  Flight 
Dynamics  Laboratory,  Wright-Patterson  Air  Force  Base,  Ohio, 
under  Contract  F33615-71-C-1310.  This  research  is  part  of  a 
continuing  effort  to  develop  new  and  improved  techniques  for 
defining  dynamic  and  aeroelastic  phenomena  for  rotor/propeller- 
powered  V/STOL  flight  vehicles  under  the  Air  Force  Systems 
Command's  exploratory  development  program.  This  contract  was 
initiated  under  Project  1370,  "Dynamic  Problems  in  Military 
Flight  Vehicles,"  Task  137005,  "Prediction  and  Control  of 
Flight  Vehicle  Vibration."  Mr.  A.  R.  Basso  of  the  Aerospace 
Dynamics  Branch  was  the  Project  Engineer. 

The  final  report  is  presented  in  three  volumes.  The  first 
volume  contains  a  state-of-the-art  review  of  stability  and 
blade  vibratory  loads  in  V/STOL  aircraft.  The  second  volume 
contains  the  development  of  the  analytical  methods,  the  corre¬ 
lation  of  analytical  results  with  experimental  data,  and  the 
results  of  parametric  investigations.  The  third  volume  con¬ 
tains  a  user's  guide  to  the  digital  computer  programs  including 
input  and  output  formats.  The  third  volume  is  not  being 
distributed;  however,  it  is  available  upon  request  from  the 
’Air  Force  Flight  Dynamics  Laboratory /FYS ,  Wright-Patterson 
Air  Force  Base,  Ohio  45433. 

Mr.  H.  R.  Alexander  was  The  Boeing  Company,  Vertol  Division 
Project  Engineer. 

This  report  covers  work  conducted  from  February  1971  through 
February  1972.  The  manuscript  was  released  by  the  authors  in 
February  1972  for  publication  as  an  AFFDL  Technical  Report. 

This  Technical  Report  has  been  reviewed  and  is  approved. 

WALTEIr  J .  l&I&TOW 
Assistant  for  Research  and 
Technology 

Vehicle  Dynamics  Division 


ABSTRACT 


The  blade  loads  analysis  calculates  rotor  blade  flapwise, 
chordwi.se,  and  torsional  deflections  and  loads,  together  with 
rotor  performance,  control  system  forces,  and  vibratory  hub 
loads.  Articulated  and  hingeless  rotors  with  2  to  9  blades 
and  large  twist  may  be  analyzed.  The  blades  may  be  of  arbi¬ 
trary  planform,  twist,  and  radial  variation  in  airfoil  section. 

The  analysis  considers  coupled  f lapwise-chordwise-torsion  de¬ 
flections  of  the  rotor  blades.  Boundary  conditions  for  either 
articulated  or  hingeless  rotors  are  applied  and  the  solution 
is  obtained  by  expanding  the  variables  in  a  10-harmonic  Fourier 
series . 

Airload  calculations  include  the  effects  of  airfoil  section 
geometry,  compressibility,  stall,  3-dimensional  flow,  unsteady 
aerodynamics,  and  nonuniform  inflow.  The  unsteady  aerodynamic 
loads  are  calculated  by  modifying  the  static  loads  resulting 
from  the  airfoil  tables  to  include  Theodorsen’s  shed  wake  func¬ 
tion,  dynamic  stall  effects  based  on  oscillating  airfoil  data, 
and  yawed  flow  across  the  blade. 

The  nonuniform  inflow  calculations  are  based  on  a  tip  and  root 
vortex  trailed  from  each  blade.  The  vortex  wake  is  assumed  to 
be  rigid  and  to  drift  relative  to  the  hub  with  a  constant  re¬ 
sultant  velocity  composed  of  thrust-induced  uniform  downwash 
and  the  aircraft  airspeed. 

The  upgraded  aeroelastic  stability  analyses  handle  all  the 
classical  mechanisms  of  instability  such  as  whirl  flutter,  air 
and  ground  resonance,  and  aeromecnanical  instability.  Addi¬ 
tional  phenomena  dealt  with  include  individual  blade  flutter 
resulting  from  pitch-lag-flap  coupling  which  may  occur  when 
rotor  blades  are  deflected  under  load;  limit-cycle  wing-rotor 
flutter  caused  by  a  similar  mechanism;  and  the  effect  on  aero¬ 
elastic  behavior  of  rotor  tilt  to  the  free  stream. 

In  the  area  of  rotor  loads,  correlation  of  prop/rotor  and 
helicopter  rotor  loads  test  data  has  been  made  with  the  fully 
coupled  rotor  loads  analysis.  Predictions  are  made  for  alter¬ 
nating  blade  loads  and  steady  in-plane  hub  moments  for  prop/ 
rotors  and  for  alternating  blade  loads,  pitch  link  loads,  and 
airloads  for  helicopters. 

In  the  area  of  aeroelastic  stability,  phenomena  beyond  the 
scope  of  analysis  with  previous  stability  methodology  are  in¬ 
vestigated  and  stability  trends  associated  with  certain  sensi¬ 
tive  parameters  are  presented.  These  include  the  influence  of 
blade  frequencies,  blade  deflections  and  thrust,  and  advance 
ratio.  Parametric  behavior  of  propeller-rotor  static  deriva¬ 
tives  is  presented.  An  unexpected  discovery  of  this  study  was 


iii 


the  sensitivity  of  the  rotor  derivatives  to  blade  lead-lag 
motion  when  the  frequency  of  this  mode  is  near  1  per  rev. 
This  result  was  subsequently  confirmed  by  test  and  adds  sub¬ 
stantially  to  our  understanding  of  the  influence  of  large 
propeller-rotors  on  static  divergence  and  flying  qualities. 

The  impact  of  the  new  prediction  capabilities  and  possible 
refinement  or  development  of  design  criteria  through  use  of 
the  analysis  are  discussed. 
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Part  I.  New  Methodology  for  Aeroelastic 

Stability  Analysis 


A.  K.  Amos 
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INTRODUCTION  ' 


SCOPE  OF  ANALYSIS 


The  new  methodology  has  been  developed  to  provide  significant 
improvements  in  areas  where  existing  tools  have  not  been 
particularly  successful  in  the  prediction  or  correlation  of 
aeroelastic  instability  behavior  (1).  It  also  provides 
greater  capability  and  flexibility  in  the  representation  of 
rotor,  airframe/empennage,  and  landing  gear  structural  inter¬ 
actions  and  is  therefore  applicable  to  entire  aircraft  sys¬ 
tems.  These  include  configurations  such  as  tilt-rotor  and 
tilt-wing  V/STOL  designs  and  single-rotor  as  well  as  tandem- 
rotor  helicopters.  Flight  conditions  of  such  systems  ranging 
from  hover  through  transition  to  cruise  are  accommodated, 
along  with  ground  effects  involving  landing  gear  dynamics. 
System  features  included  in  the  analysis  are: 

o  Large-diameter  rotors  with  flexible  blades 

o  Airframe  structure  with  flexible  aerodynamic 
lifting  surfaces 

o  Landing  gear  with  linear  viscoelastic  properties. 

Starting  with  a  description  of  the  system  and  flight  condi¬ 
tions  and  using  basic  aeroelastic  properties  of  the  system 
and  environment  as  supplied,  the  analysis  establishes  a  number 
of  possible  states  of  perturbation  about  steady  flight  condi¬ 
tions  and  the  inherent  stability  of  each. 

FEATURES  OF  ANALYSIS 

The  analysis  provides  for  large  aeroelastic  steady  deflections 
of  the  rotor  blades,  which  may  include  mechanical  precone  and 
static  droop.  These  deflections  are  one  source  of  kinematic 
and  dynamic  couplings  between  blade  bending  and  blade  tor¬ 
sional  vibrations.  The  effects  of  these  couplings  on  the  sta¬ 
bility  of  the  system  have  been  discussed  in  Reference  1.  Their 
inclusion  in  the  analysis  represents  a  significant  advance  in 
the  state  of  the  art. 

Lumped  mass  and  inertia  representation  and  modal  synthesis  of 
the  system  perturbations  from  steady  state  form  the  basis  of 
the  analysis.  Fully  coupled  pitch-flap-lag  blade  modes  are 
provided  for.  The  aerodynamics  of  the  rotor  blades  and  air¬ 
frame  lifting  surfaces  are  based  on  two-dimensional  quasi¬ 
steady  theories,  but  with  provision  for  reducing  the  blade 
aerodynamics  to  a  steady  (or  static)  formulation.  The  air  in¬ 
flow  through  the  rotor  disc  is  treated  very  generally,  provid¬ 
ing  for  any  inclination  from  axial  to  purely  tangential.  This 
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eliminates  one  of  the  most  serious  application  restrictions 
inherent  in  existing  analyses. 

LIMITATIONS  OF  ANALYSIS 

The  consideration  of  steady  blade  deflections  is  limited  to 
collective  effects  only;  i.e.,  cyclic  blade  deflections  are 
assumed  negligible.  Cyclic  perturbations  are  considered, 
however.  The  blade  aerodynamic  treatment  does  not  account 
for  stall  effects,  although  provision  is  made  for  the  use  of 
experimentally  determined  lift  curve  slope  values. 
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SYSTEM  DEFINITION 


STRUCTURAL  IDEALIZATION 


The  entire  system  is  partitioned  structurally  into  three  sub¬ 
systems:  rotor,  airframe,  and  landing  gear.  Each  subsystem 

is  made  up  of  a  number  of  available  components  which  serve  as 
building  blocks  for  defining  particular  problems  during  actual 
application  of  the  analysis. 

The  landing  gear  subsystem  has  four  possible  oleo-tire  assembly 
units.  Each  is  idealized  as  a  rigid  mass  supported  by  a  number 
of  linear  springs  and  dashpots  as  shown  in  Figure  1.  The 
springs  and  dashpots  represent  the  oleo  and  tire  viscoelastic¬ 
ity.  The  point  P  lies  in  and  has  compatible  displacements  with 
the  airframe  subsystem;  it  is,  however,  isolated  rotationally 
from  the  airframe.  The  mass  has  an  independent  translational 
freedom  parallel  to  the  Z  axis  (Figure  1) .  In  the  other  two 
coordinate  directions  it  follows  rigidly  the  translations  of 
point  P. 

Two  general  configuration  rotors  are  provided  for  in  the  rotor 
subsystem.  They  can  be  particularized  to  tandem-rotor  heli¬ 
copter,  tilt-rotor,  tilt-wing,  and  conventional  propeller- 
driven  aircraft  configurations;  or  to  a  single-rotor  helicopter 
system  by  retaining  only  one  of  the  two.  Each  rotor  consists 
of  three  or  more  blades  attached  to  a  central  hub.  The  hub  is 
mounted  to  a  rotating  shaft  which  is  encased  in  a  nacelle  or 
pylon  structure  as  illustrated  in  Figure  2.  The  hub  and  shaft 
combination  is  isolated  from  the  nacelle/pylon  in  rotation 
about  the  shaft  axis.  Additionally,  the  hub  is  free  to  teeter 
in  pitch  and  yaw  relative  to  the  nacelle/pylon. 

The  blades  are  mounted  to  the  hub  with  variable  offsets  from 
the  hub  center  of  rotation  as  illustrated  in  Figure  2.  Each 
blade  is  flexible  in  both  bending  and  torsion  and  is  idealized 
as  a  slender  elastic  beam.  Its  steady-state  configuration  is 
characterized  by  possibly  large  deflections  in  flap,  lead-lag, 
and  pitch,  including  geometric  twist  and  collective  pitch 
setting.  The  hub  freedoms  of  motion  are  made  up  of  (1)  all 
translations  of  the  nacelle/pylon  structure  at  its  interface 
with  the  hub;  (2)  nacelle/pylon  rotations  other  than  that  about 
the  shaft  axis;  and  (3)  teetering  and  shaft  axis  rotational 
freedoms  of  the  hub.  The  blade  kinematics  are  separated  into 
rigid  components  arising  from  the  hub  motions  and  elastic 
deformations  relative  to  the  hub. 

The  following  components  constitute  the  airframe  subsystem: 
two  wings,  two  tails,  one  fin,  and  two  nacelle/pylon  struc¬ 
tures.  The  wings,  tails,  and  fin  are  each  considered  as  a 
two-dimensional  lifting  surface.  Each  is  idealized  as  a 
finite-element  aeroelastic  system. 
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Figure  1.  Typical  Landing  Gear  Representation 


i. 


Figure  2.  Schematic  Diagram  of  Rotor  Hub  Showing  Blade 
Mount  Eccentricities 


MATHEMATICAL  MODELING 


The  mathematical  model  of  the  system  envisions  the  aircraft  in 
its  steady-state  configuration  as  fixed  in  space  with  the  air 
streaming  past  it.  A  set  of  coordinates  fixed  to  the  aircraft 
in  the  steady  state  is  therefore  also  fixed  in  space. 

Tne  structural  idealization  of  the  landing  gears  also  defines 
the  mathematical  model  of  the  subsystem.  In  the  airframe  and 
rotor  subsystems,  the  mathematical  model  is  derived  from  a 
finite-element  discretization  of  the  individual  components. 

The  inertial  and  aerodynamic  properties  of  the  elements  are 
lumped  at  representative  points  or  sections  and  the  elastic 
properties  are  averaged  between  adjacent  elements.  As  an 
alternative  to  considering  the  actual  elastic  properties,  the 
model  postulates  the  existence  of  characteristic  functions  or 
modes  which  reflect  the  inertial  and  elastic  properties  of  the 
finite-element  system.  The  kinematic  freedoms  of  the  finite 
elements  are  therefore  represented  by  modal  synthesis,  em¬ 
ploying  modal  freedoms  as  generalised  coordinates. 
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COORDINATE  SYSTEMS  AND  TRANSFORMATIONS 


FIXED-COORDINATE  SYSTEM 


A  system  of  orthogonal  coordinates  is  established  in  space 
such  that  the  X  axis  is  coincident  with  the  steady-state 
longitudinal  axis  of  the  aircraft.  The  origin  is  located  at 
some  convenient  location  along  this  axis.  The  Z  axis  is  such 
that  it  would  represent  a  vertically  downward  direction,  if 
the  aircraft  were  in  level  flight.  Figure  3  illustrates  the 
fixed-coordinate  system. 

SPATIAL  COORDINATES  IN  AIRFRAME  SUBSYSTEM 


Each  finite  element  of  the  airframe  components  is  associated 
with  a  reference  point  whose  spatial  coordinates  in  the  fixed- 
coordinate  system  under  steady  state  are  represented  by  the 
vector  {XFR}0.  A  local  coordinate  system  is  defined  with 
origin  at  the  reference  point  under  steady-state  conditions. 
Transformation  from  the  local  to  the  fixed-coordinate  systems 
is  represented  by  the  matrix  [TAF]  . 

Figure  4  illustrates  the  local  axes  for  typical  elements  of 
the  various  airframe  components.  Except  for  the  fin  element, 
each  local  system  shown  represents  a  rotation  of  the  fixed 


system  through  a  pitch  angle 
transformation  is  thus  given 

about  the 
generally 

Y 

by 

axis. 

The  coordinate 

-  cos  P, 

0, 

sin 

P  - 

tTAF]  - 

0  , 

1, 

0 

(1) 

.-sin  P, 

0, 

cos 

P  - 

where  P  =  Pw  for  wing  element 

=  Pps  for  fuselage  element 
=  P-p  for  tail  element 
=  0  for  fin  element 

Under  perturbations,  the  reference  point  is  displaced  and 
rotated  relative  to  its  steady-state  location.  These  are  re¬ 
presented  in  the  local  coordinate  system  by  the  displacement 
vector  {dR}  and  the  rotation  vector  {0R}. 

The  spatial  coordinate  vector  of  the  reference  point  in  a 
perturbed  configuration  is  thus  given  in  the  fixed-coordinate 
system  by 


{XFr1  -  {XFR}0  +  [TAF]{dR}. 


(2) 
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Figure  4.  Local  Coordinate  Axes  for  Typical  Segment  of 
Various  Airframe  Components 


Consider  now  a  point  A  in  the  element  whose  coordinates  rela¬ 
tive  to  the  reference  point  under  steady-state  conditions  is 
given  in  the  local  coordinate  system  by  the  vector  {AAK  Its 
perturbed  spatial  coordinate  vector  in  the  fixed  system  is 


where 


{XFa^  ~  {xFR^o  +  tTAF^  +  [t6r]{^A^)  (3) 


1  , 

"8rz  ' 

6ry 

RZ  ' 

1  , 

"9rx 

RY ' 

8rx 

1 

See  Reference  2  for  derivation  of  transformation. 


Modal  synthesis  of  the  displacement  and  rotation  vectors  {dR] 
and  {0R}  is  expressed  by 

{dR}  =  l  {*df}j  ^Aj 

j 

{  6r}  =  l  Urf>  .  <JAj  (5) 

j 

where  { } j  is  a  vector  of  the  displacement  components  and 
{(j)rf}j  is  a  vector  of  the  rotational  components  of  the  jth 
airframe  mode  specified  at  the  reference  point. 

Variations  of  the  spatial  coordinates  with  the  modal  freedoms 
can  now  be  established  as  below: 


[TAF^  ^df^j 

[TAFHHdf  J  j  +  tT0R],j{AA})  (6) 


where  [T0RJ,j  -  ^  [T0R]  - 

0 

t 

_4>rf  j  (3)  , 

<t>rf  j  (2) 

<J>rf  j 

(3)  , 

0 

"#rf  j (1) 

_*rf  j 

(2)  , 

4>rf  j  ^  1)  • 

0 

i  _  3{XFR} 
lXFR]'j  “  IqJT- 

,  ■,  3  ^  XFA^ 

{XFA}'j  =  ~9qA j 


(7) 
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SPATIAL  COORDINATES  IN  ROTOR  SUBSYSTEM 


Figure  5  shows  the  rotor  shaft  axis  in  relation  to  the  fixed 
coordinate  system  under  steady-state  conditions.  The  hub 
location  on  this  axis  is  used  to  establish  the  origin  of  a 
steady-state  shaft  coordinate  system  as  shown  in  Figure  6. 
The  coordinate  transformation  from  this  system  to  the  fixed 
system  is  represented  by  [Tp]  given  by 


[Tn]  = 


-sin  PQ,  0  ,  cos  PQ" 
0,1,  0 
-cos  PQ,  0  ,  -sin  P0 


(8) 


In  a  perturbed  configuration,  the  hub  is  displaced  and  rotated 
relative  to  its  steady-state  location.  These  are  represented 
by  a  displacement  vector  {dH}  and  a  rotation  vector  { 6 pj } 
referenced  to  the  steady-state  shaft  coordinate  system. 


A  perturbed  shaft  coordinate  system  is  now  defined  with  its 
origin  at  the  perturbed  hub  location  and  representing  a  rota¬ 
tion  of  the  steady-state  shaft  system  through  the  rotation 
vector  { ©h) •  The  transformation  between  the  perturbed  and 
steady-state  shaft  coordinate  systems  is 


1  , 

_eHZ  ' 

0HY 

HZ  ' 

1  , 

_0HX 

HY' 

8HX' 

1 

(9) 


In  the  perturbed  configuration,  the  hub  is  considered  teetered 
so  that  the  disc  plane  is  not  normal  to  the  shaft  axis  at  the 
hub. 

(3)  (3)  (3) 

A  nonrotating  disc  coordinate  system  (X,  Y,  Z  )  is  defined 

(3) 

with  its  origin  at  the  perturbed  hub  location  and  the  Z  di¬ 
rection  normal  to  the  disc  plane.  It  represents  a  rotation  of 
the  perturbed  shaft  coordinate  system  through  a  teeter  angle 

(4)  (4) 

with  components  $TEX  ab°ut  X  and  a^out  Y  •  The 

transformation 


[Tte] 


'cos  $^ey;  3/2  (sin  4>tey^  (sin  $TEx)  >  sin  $TEY 
1/2 (sin  4>tey)  (sin  $TEX) !  cos  $TEX;  -sin  $TEX 
—sin  'f'T'EY i  sin  ^^EX'  (cos  “^teY^  (cos  ^TEX^ 


(10) 
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ROTOR  SHAFT  AXIS 


Figure  5.  Rotor  and  Aircraft  Coordinate  Axes 


Figure  6.  Shaft  and  Disc  Coordinate  Systems 


relates  the  nonrotating  disc  system  to  the  perturbed  shaft 
coordinate  system.  A  rotating  disc  coordinate  system  is 
established  by  rotating  the  nonrotating  system  through  the 
azimuth  angle  ip  as  shown  in  Figure  7.  Transformation  from 
rotating  to  nonrotating  system  is  given  by 


'V  - 


sin  \p  cos  ip  0 

-cos  sin  ip  0 

0  0  1 


(11) 


The  displacement  vector  (XBA)  referred  to  the  rotating  disc 
system  locates  a  predetermined  reference  point  on  the  chord 
axis  of  a  representative  blade  segment  section.  The  plane  of 
the  blade  section  is  used  to  define  a  normal  blade  section 
coordinate  system  as  shown  in  Figure  8.  The  transformation 
from  this  to  the  rotating  disc  system  is  given  by 


[Tsl 


(l~l/2o£)  /  ©L  / 

-eL  ,l-l/2(02  +  e2),  _9[, 

-V2Vf  '  eF  '  (1-1/2eF) 


(12) 


A  local  blade  coordinate  system  (XB,  YB,  ZB)  with  XB  parallel 
to  the  chord  axis  is  next  defined  with  the  origin  at  the  ref¬ 
erence  point.  This  represents  a  rotation  of  natural  blade 
section  coordinate  system  by  the  blade  section  twist  angle  eT 
(1) 

about  Y  .  Hence  the  transformation  matrix 


[Tq] 


cos  0T,  0  sin  6rp 

0,1  0 

-sin  0ip,  0  cos  0ip 


(13) 


Consider  now  a  point  A  in  the  blade  section  with  coordinates 
in  the  blade  coordinate  system  represented  by  {AA}.  In  the 
normal  blade  section  system  it  has  coordinates  represented  by 


(1) 

(XA)  =  [Tq  ]  UA}  .  (14) 
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In  the  other  coordinate  systems  point  A  has  the  coordinates 
shown  below: 


Rotating  Disc  System: 

(2) 

{XA}  =  {XBA}  +  tTs] [T0 ] { XA}  (15) 

Nonrotating  Disc  System: 

<V  -  iV  “V1  +  iVITeHV)  <16) 

Perturbed  Shaft  System: 

(4) 

«A>  =  tTTE] [T„] ({XBA}  +  [Ts](T9]{Xa})  (17) 

Steady-State  Shaft,  System: 

(5) 

(XA)  =  (dB)  +  [Tw]  [Tipg]  [T^]  ({XBA)  +  [Ts]  [T0]  {AA})  (18) 

Fixed  System: 

^XFA^  =  +  ^Tw^TTE^T^^  ^XBA^  +  [Ts^Te^AA^ 

(19) 

The  blade  coordinates  are  made  up  of  steady-state  and  pertur¬ 
bation  components.  The  latter  are  obtained  from  modal 
synthesis.  Thus  the  blade  twist  angle  is  given  by 

eT  = (Y  +  I  eRj  ^Bj ) r  ,  (20) 

j 

where  y  is  the  net  steady-state  twist  angle  locally; 

0Rj  is  the  blade  twist  component  of  the  jth  blade  mode.  The 
components  of  {XBA>  are  given  by 

XBAX  =  (?X  +  l  nxj'  3Bj)  r 

j 

XBAY  =  r  “  Sy  "  I  nyj  ^B j  “  1  ^jk  ^Bj  ^Bk  (21) 

j  jk  ' 

XBAZ  =  ?z  +  l  nzj  qBj 

j 

T  =  +1  for  clockwise  rotation  of  the  blades  as 
viewed  from  the  shaft  side  and  -1  for 
counterclockwise  rotation 


where  ;x  and  ?z  are  the  steady-state  blade  deflections  in  lag 
and  flap;  is  the  radial  shortening  due  to  the  above  deflec¬ 
tions;  r  is  che  radial  distance  of  the  segment  reference 
section; _nxj  and  nzj  are  modal  components  in  lag  and  flap,  and 
rtyj  and  rij^  are  modal  contributions  to  the  blade  radial 
shortening. 

The  angles  9l  ana  0p  are  slope  angles  associated  with  the  lag 
and  flap  displacement  components; 

1  * e • '  eL  =  XBAX , r  =  (5x,r  +  l  nx,rj  qBj)r 


eF  ~  xBAZ,r  “  ?z,r  +  l  nz,rj  qBj 


Using  equations  20  and  22  in  equations  12  and  13  yields 


IT  ]  = 

6  o 


cos  y  Or  sin  y~ 


-  r  sin  y  0  cos  y 


(~-l [T0])  = 

9c3Bi  0  o 


l  ,  o,  reRi 

o,i,  o 


-  r6Ri'  °-  0 


tTj  = 
b  o 


-!/2n;Xjr  52,r! 


[l-l/2(4>r+^,r^ 


1//2r^x,r^z , 


(l-l/2c^r) 


[Ts])  =  - 


n  . ;  r n  -l/2r(c  n  .  +?  n 

x,r  x,ri  x,ri  x,r  z,n  z,r  x,rx 


i  ;  -  U  n  .  +  e  n  . ) ;  -n 

x,n  x,r  x,n  z,r  z,ri  z,ri 


■  1/2  r  ( ?  n  .  +  c  n  . )  ;  n  .  ;  -  x>  n 

x,r  z,n  z,r  x,-ri  z,ri  z,r  z,n 
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—a*  ,335S*f s3is«*SiSX*te*=*s-to 


From  equation  21 


<xba> 


z  J 


r rn 

'xi 


(9qBi  {Xba>)o  \  "V  ) 


Zl- 


( 


.  .  { Xq a})  —  j  -n  . 

9qBi  9qBj  °  \  3 


'  0 
-n -k 


“A 


J 


The  azimuth  angle  ij;  for  an  individual  blade  is 

*  =  +  rnt  + 

where  q,.,  is  the  rotor  hub  rotational  freedom  about  the 
axis.  Hence: 


[TJ  = 


sin  (if>0  +  Tfit)  ;  cos(^0+rfit);  0" 

-cos  ( ^ 0  +  rfit)  ;  sin  (if»0  +  rfi  t)  ;  0 

o  0 


tJ 


at  [I>])  =  ™ 

0 


cosU'0  +  rfit);  -sin  ( (j/0  +  Tfit)  ; 
sin(ij>  +ffit);  cos(^0+rnt); 


3t< 


[Tx]J  =  -Q2 


sin  (ij;0  +  Tftt)  ;  cos  U'0  +  Tfit)  ; 


0 


-cos  (^0  +  Tftt)  ;  sin  (ij>0  +  Tftt)  ;  0 

0  0:0 


 J*W'*)*m  -li 

-ii  fi  -  -**i2EfeE5wi^j«;- 


(24) 

(25) 
shaft 
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(9q,  3t  llV*  ~ 
^4*  o 


-ra 


^  ^2 


[T,.  ] )  =  -a2 


sin(4»0  +rat)  ;  cos  (4i0  +  rat)  ; 

-cos  (i|/0  +  rat)  ;  sin  (4>0  +  rat)  ; 

o  0; 

cos  (4/0  +  rat) ;  -sin(4>  +  rat)  ; 
sin  (4j0  +  rat) ;  cos  (ip  +  rat)  ; 


0 

0 

0. 

0 

0 

0- 


a<3t h 


[T  J) 


cos  (4j0  +  rat)  ;  -sin  (ip0  +  rat) ;  0 

sin  (ip0  +  T  at)  ;  cos  (i|/0  +  rat)  ;  0 

0  0  ;  0 


tT*n0  ■ 

Ip 


-sin(4>0  +  rat) 
cos  (4*0  +  rat) 
o 


-cos (\p0  +  rat) ;  o 
-sin  (4>q  +  r  at)  ;  0 


0 


0 


(26) 


The  hub  teeter  angles  4>TEX  an<3  4>TEY  have  steady-state  and  per¬ 
turbation  components; 

i.e.,  4>Tex  =  ^TE  +  ^TEX 

4>tey  =  eTE  +  ^TEY  (27) 

Substituting  in  equation  10, 

cos  SfE'  1/2  sin  0ij£  sin  4'TE'  sin  ®TE 
^tte^0  =  !/2  sin  ©te  sin  ^te !  cos  4>te'  -sin  4»te| 
_-sin  6rpg;  sin  4>te '  cos  ®TE  cos  4'TE 


a<?TEX 


1  ^TE^  ^  -  ~  tTijigJ 


0  0  0 

0  0-1 

0  10 
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[T,pp ]  )  —  [Trpg ] 


9qTEY  TE  o 


0 

0 

-1 


0 

0 

0 


1 

0 

0 


9qTEX 

3qTEX 

L  TE J  ' 

o 

3  2 

[ttej ) 

O 

9qTEX 

9qTEY 

3  2 

[tte])  = 

3qTEY 

9qTEY 

ltej 


lte 


lte- 


0 

0 

0 

0 

1/2 
_  0 

'-1 

0 

0 


0 

-1 

0 

1/2 

0 

0 

0 

0 

0 


0 

0 

-1_ 

o" 

0 

0_ 

o’ 

0 

-1J 


(28) 


The  shaft  perturbation  displacements  and  rotations  are  repre¬ 
sented  by: 


{V  =  l  4df} .  qAj 

j  3 

{ 6H}  =  l  ^rf }  .  qAj 

j  3 

Therefore,  from  equation  9, 


‘V  * 

[I] 

o 

0 

<f>rfi  (3)  ? 

4>rfi  (2) 

(— -  [T  ]  )  = 

9qAi  W 

<t>rf  i  C  3 )  ; 

0 

-(J’rfi  U) 

-'t’rfi (2)  ; 

♦  rfi (1)  ; 

0 

(29) 
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=■*  *>  ->  jJWS'j  ~J  tt  ISC?*/*  *. 


Also , 

{dH} 

=  {O} 

o 

(t—  (dH>) 
aqAi  H  o 

=  Hdf\ 

l 

(30) 

SPATIAL 

COORDINATES 

IN  LANDING  GEAR  SUBSYSTEM 

The  spatial  coordinates  of  the  reference  point  P  are 

{Xpp}  =  {Xpp}  +  'l  q^j' 

°  3  3 

The  mass  point  spatial  coordinates  are  given  by 

0 


f  $df  (x)] 

{XFH}  =  l  V ^dQ(y)J  qA j  +|  J  jqL 


(31) 


32) 


3 


{ 4>df }  j  is  a  vector  of  displacements  at  P 
frame  mode. 

qL  is  the  individual  vertical  freedom  of 


due  to  the  j th  air- 


the  mass. 
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ANALYTICAL  METHODS 

LAGRANGIAN  FORMULATION  OF  STABILITY  EQUATIONS 

Lagrange's  equations  of  motion  for  nonconservative  systems 
ma"  be  written  as 

d  (9T_)  _  3T_  +  3U_  +  3R_  _  3We  _  Q 
dt  3Uj_  3Uj_  3u^  3Uj_  3u  ' 


where  T  =  kinetic  energy, 

U  =  strain  energy, 

R  =  dissipation  function  of  the  system. 

We  =  virtual  work  of  external  forces, 

Uj_  =  ith  component  of  generalized  coordinates  of  the 
system. 

Application  of  the  above  to  a  dynamic  system  leads  to  equili¬ 
brium  equations  which  can  be  represented  generally  by 

Fi  (Uj ;  Uj  t  Uj ,  t)  =  0  (34’ 

and  are  generally  nonlinear  in  the  generalized  coordinates. 

The  coordinates  for  an  aircraft  system  can  be  separated  into 
steady-state  and  perturbation  components.  The  latter  may  be 
considered  as  a  linear  combination  of  known  characteristic 
functions  with  their  undetermined  coefficients  constituting 
the  perturbation  freedoms? 


u .  --  u 
1  oi 


+  I  *ik 


where  uQj_  =  steady-state  component  of  ith  coordinate, 

<j>ik  =  kth  characteristic  function  in  the  ith  coordinate, 
qk  =  kth  perturbation  freedom. 

Hence  equation  34  becomes 

Fi  (^o j  '  ^oj  '  ^oj  '  ?  t)  —  0  .  (36) 
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ifiiffiKpfc'sS  .■-■.Aakfet-  git  ^fcaiSs<3£ifeWwii^&6ft ^m*^*^*^* 


A  Taylor* s  series  expansion  leads  to 


3F .  3F-  3F. 

(F.)  +  J  [(— i)  q,  +  (-i)  q.  +  (-^)  qU 


3q-i 
Do  Jo 


32F. 


3q 


J  o 


3q 


J  o 


32F; 


2p  . 


3ZF 


n  n  0  F  -I  0  A  s  0*1  • 

+ 1/2  1  "^ir*  q3qk  +  ‘tttt’  q^k  +  q’qk 

jk  3qj3qk  o  3qj3qk  o  3qj  3qk  o 


+  ( 


32F.  ..  3  2F  •  ..  3  2F  •  .  . 

— 4-)  g^qk  +  (-— — )  q^qk  +  q^v 

8c3j3<3k  o  J  9qj3qk  o  3  ^qj^qk  o  3 


+  (■ 


3  2F  • 


-) 


3  2F 


32F, 


qjqk 


+  <—■-:->  ijqk  +  (■— -V>  9-j^k 

3qj5qk  o  a<5 j  8<Jk  o 


3qj3qk  o 

+  (higher  order  derivative  terms)  =  0  .  (37) 


The  zero  subscripts  in  the  above  relation  indicate  evaluation 
of  the  quantities  involved  for  steady-state  equilibrium  con¬ 
ditions.  These  conditions  are  governed  by  the  equilibrium 
equations 


(u 


°1 


u 


oj '  uo j 


t)  =  0 


(38) 


The  remaining  terms  of  equation  37  must  then  vanish  to  ensure 
equilibrium  of  the  system  in  any  perturbed  configuration  — 
this  is  the  condition  for  neutral  stability  of  the  steady 
state.  Thus,  when  attention  is  focused  only  on  small  per¬ 
turbations,  the  stability  equations  are  given  by 

I  [  OFi/Sq-j^qj  +  OFi/SqjJ^j  +  ( 3F.j/3q  j )  Qq  j  ]  =  0  .  (39) 

j 


In  matrix  notation  the  above  is  represented  by 


where 


[M] {q}  +  [C] {q}  +  [K] {q }  =  Q 


M  ( i ,  j ) 

3Fi 

1  =  — n — 

3q  j 

C (i,  j] 

3Fi 

(40) 


25 


3F-: 

K  (i  f  D  )  =  ~ 

9qj 

{q}  =  vector  of  all  system  freedoms. 

KINETIC  ENERGY  TERMS 

The  kinetic  energy  of  the  system  is  expressed  generally  by 


T  =  1/2  If  pm  {Xf)T{Xf}  da  , 


(41) 


where 


Xp  is  a  spatial  coordinate  vector  of  an  element  of 
mass , 

pm  is  the  mass  density, 
dc  is  an  element  of  volume. 


The  integration  extends  over  the  volume  of  a  finite  element 
and  the  summation  is  over  the  totality  of  elements  in  the 
system. 

Time  dependency  of  the  spatial  coordinates  is  both  implicit 
through  the  perturbation  freedoms  (q^)  and  explicit  through 
the  azimuth  angle  increment  (Qt)  of  the  rotor  blades? 


i.e. ,  {XF}  =  {Xp(q£,t) }  .  (42) 

Hence,  {Xp}  =  ^  {XF}  +  £ -~-{X p}  q.. 

j  ^ 


(XF  -qi 

,  qi#  t)  } 

— (XF) 

+  I  «lq7 

j  3 

{Xp}  qj  +  .  { Xp }  qj 

3q  j 

a2 

Bt?{XF 

S2 

}  +  T  (2— 

L  3q 

— {xf}  +  fi7<V 

j 

D  D 

+  I  I 

j  k 

3  2 

3qjaqk(XF 

1  qjqk 

{Xp  (qi 

'  qi'  qi' 

t)  } 

(43) 


The  following  relations  follow  from  the  above  results: 


-  'V  =  fr  (XF> 


3q. 


3q. 


<3  ,  *  ,  3 2  • 

8q j  “  3q j  8t  {XF}  +  1  aqTTq^  {XF*  <3j 


(44) 


From  equation  41, 


8T  v  ,  8  *  T  * 

•  ~  l  J  Pm  7T~  ( XF  ^  da 

°<Sj  3qj 


-  I  /  Pm  I(gq.(xF^  (^{Xp})  +  I  (|— -{xp}) 

1  i  ^  1 


(8q,  {XF})<3k]  ^CT 

K 


d  ,3T 
dt  39j 


.  v  r8  /9T  1  *  .  3  ,8T  . 

8t  3q j  8qk  9qj  3qk  3qj  k 


I  /  pm  t ( 


+  ^((l5^ttXF))  (%{XF»9k  +  2(f57(XF>>T<f!^(XFi)c|k 

k 


+  (■ 


3qi^qJ{XF})T(lt{XF})  +  (|^-{XF> ) 1  (f^~{XF> )  9k} 


T  ,8 


+  ^<(|q-(XF})T(|l^7{XF})+  ( 


k  1 


^qk9qi 


8q j  8qi 


{XF>)T<lqjtxF))> 


V^l1 


da 
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Alt..'  , 


3T 

3qj 


I  /  Pm  <|57<V>T{V  d<1 


a  2  T  S { Xp  }  2 

0  r>,  M1,  F  \  ,  VT  ,3 


=  i!  pm  +  n  ( 

3 


3t 


{X„}) 


3q  .  9q,  F 
kl  3  k 


9  r  ■*,  TV  •  •  r  r/3  r  1  \  ^  /  3 


(377~{xf}>  <3,.  ?,  +  I  {(■ 


*1  ‘  4  lV3i~3t{XF})  (3^{XF}) 


+  ( 


T  ,  3 


3q.3q,{XF})  (3t{XF})}  ^k^  do 

D  K 


d  3T 
Hence,  ^r(~ 


3T 


3q^  3q. 


I  /  pm  L(gq^Xp})  (3t.2{XF})  +  H  (|^— {Xp})T 

kl  3 


9  2 


^g^-{Xp})  qkqi  +  l  {(^-{Xp})T(^-{Xp})  qk 

V  «' 


+  2^3q.^XF^  ^3qv3t{XF^  ^k^  da 


The  right-hand  side  of  the  above  is  now  expanded  in  a  Taylor' 
series  about  the  steady-state  configuration  and  truncated  at 
the  linear  terms  in  the  system  freedoms  and  their  time  deriva 
tives  to  give: 

d  |ST  |  _  3T_ 

dt  3<q j  3q j 

S 

3  T  3 

—  I  ^  I/p m  ^  9q  ^  Xp ^ )  ("g^j  t Xp } )  da}  q^. 

■»_  j  o  k  o 


+  2l  { /pm 

k 


{xF}) 


T 


at 


(XF}) 


da}  qk 


0 


0 


0 


+  l  <J Pm  [( 

k  j  k 


a2  T 

o  r  v  i  \  /  o 


m  1  l3q'lq"{XF})  (7EJ{XF}) 


+  <Wr(xF»T'l^trfxF))  1  <J»>  q, 

3  0  k  0 


+  h m  (f— (Xf))T(^(Xf})  do] 
qj  0  0 


(45) 


The  kinetic  energy  contributions  to  the  stability  equations 
are  therefore: 

S 


MKE ( 1 ' 3  ' 

1  ~  II  Pm  ( 3q£  {XF^0  (3qj  {XF})0  do 

‘“KE  (i  '  3  • 

S  3  T  ^2 

1  pm  <3q.  (XFn  < 3q j st  <XF»0 

KKE  ( i ' 3 ' 

s 

1  -  l(sqi5qj  (XFn0  (3t2  IXf))0 

+  (3  (XF>)  (  t2  {XF})  ]  da 

(46) 

The  last  integral  in  equation  45  is  the  kinetic  energy  con¬ 
tribution  to  the  steady-state  equilibrium  equations. 

EXTERNAL  LOAD  TERMS 


The  external  loads  acting  on  any  finite  element  of  the  system 
are  represented  generally  by  a  force  vector  and  a  moment 
vector.  The  static  equivalents  of  these  at  the  origin  of  the 
fixed-coordinate  system  are  denoted  by  {Ff}  and  {Mp}.  Under  a 
virtual  displacement  comprising  perturbations  from  steady 
state,  the  finite  element  undergoes  conf iguralioi.  changes. 

The  actual  displacements  and  rotations  of  the  element  have 
kinematical  equivalents  at  the  origin  of  the  fixed-coordinate 
system  represented  by  6{Xp}  and  6{<j>p).  These  are  defined  as 
the  displacements  and  rotations  which  must  be  applied  at  the 
origin  to  produce  the  true  deflections  and  rotations  of  the 
element  if  it  were  rigidly  linked  to  the  origin  in  the  steady 
state . 
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The  virtual  work  of  the  external  forces  on  the  element  are 
then  obtained  as: 

Wp  =  ( 6  {X„}T{F  }  +  (6{(}>})T{  }  .  (47) 


But  5 { Xp } 

=  I 

‘4 

(XF»0 

qj 

j 

<5  ( <P  F } 

-  I 

(-i- 
®q  j 

uFno 

qj  ? 

(48) 

j 

IK 

j 

CS 

T  T 

therefore , 

W 

e 

T"~ — {  Xp  } 

3q_j  F 

)  {Fp}  +  O/aq-Wp})  {Mp}]  q. 

o  o  J 

and 


3W 

jf  =  +  <33-Uf))T(mf}  . 

1  DO  DO 


(49) 


If  it  is  assumed  that  the  forces  and  moments  are  functions  of 
the  system  freedoms  and  their  first  two  time  derivatives, 


i  .  e .  , 


{ Fp }  =  (Fp  (qi,  qi,  t)  } 

{Mp}  =  Mp  (q±,  q±l  q£/  t) }  , 

then  a  Taylor's  series  expansion  of  9We/9qj  yields 


(50) 


3We  =  { x„) )  T{Fr,}  „  +  (-1_{^})T{M„} 


3q j  '^“F"0  l‘F  0  3<3j  F"0  l“FJ° 


T,  a 


T  ,  9 


+  I  ^  9a  ^la- {FF^  +  <  9a  {<^F^  4a  {MF^ 


8qi 

k  3 


o  9qk 


3qj 


9q, 

o  k 


T,  9 


T  ,  9 


+  lU—  {Xp})  (  7  ( Fp } )  +  ( - UF1)  (— {Mf})]  qk 

L  Srr  .  t  I  9g_.  ^  3rt.  *  K 


k  9qD 


3qk 


o  9(3k 
T-  9 


+  l  [  (-^-{XF})T(-IU{FF})  +  (— UF>>*  (-^-{Mf})  ]  qk 

o  3c3k  o  9t4  o  9(3k 


k  9<4 


+  higher  order  derivative  terms, 


(51) 
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The  external  load  contributions  to  the  stability  equations 
are  obtained  from  the  above  as 

S 

MEX  (  if  j  )  =  I  t  {  Xp  }  )  (■■■::■  ■{  Fp  }  )  +  (— — {<|>p})  (— r. — {Mp})  ] 

3qi  o  3q j  o  3qi  o  3qj  o 

S 

cEX(i,j)  =  IK— {xF} ) T (— |_{Fp} )  +  (-i_Up})T(-i-{MF})  ] 

Sqp  o  3q  j  o  3q^  o  3q^  o 


'Ex(i-i)  -  H<a|7«F»0T(^ 


■{Fpj )  ^  +  (^UpnJf^Mp))^ 


(52) 


The  external  loads  considered  in  the  present  analysis  comprise 
aerodynamic  forces  and  moments  on  the  rotor  blades  and  on  the 
airframe  components  (except  fuselage) .  The  blade  aerodynamics 
are  considered  below  in  illustration  of  the  external  load 
analysis . 


The  aerodynamic  forces  and  moments  on  a  blade  segment  consist 
of  a  lift  force  AL,  drag  force  AD,  and  pitching  moment  AM  as 
shown  in  Figure  9. 


They  are  expressed  by 

AL  =  l/2pU2 (2b) Ar .  CL 
AD  =  l/2pU2 (2b) Ar .  CD 

AM  =  1/2PU2 (2b) 2Ar.  CM  (53) 

The  coefficients  CL,  CD,  and  CM  are  functions  of  the  airspeed 
(U) ,  the  angle  of  attack  (a)  and  its  time  derivatives  (a,  a), 
the  segment  velocity  and  acceleration  (h,  h)  normal  to  the 
airstream.  The  airspeed  and  angle  of  attack  depend  on  the 
velocity  components  U<p  and  Up  shown  in  Figure  9  as  follows: 

U  =  (Up2  +  UT2) 1//2 

a  =  0T  -  tan-1  Up/UT  (54) 

Hence  the  aerodynamic  independent  variables  may  be  taken  as: 

Ut ,  Up ,  b(ji ,  Up  ,  Uip ,  Up ,  0ip ,  0,p ,  9<p ,  h,  h  •  (55) 


These  in  turn  are  functions  of  the  system  freedoms. 


„L_ 


The  velocity  components  UT  and  Up  are  derived  from  a  velocity 
vector  {Vg}; 


where  is  normal  to  both  U<j>  and  Up. 


The  vector  is  made  up  of  components  from  the  blade  rotational 
speed  ft,  the  aircraft  speed  with  components  Vx»  Vy,  and  Vg  in 
fixed-coordinate  system,  and  induced  velocities  from  various 
sources  represented  in  the  rotating  dine  coordinate  system  by 


(57) 


It  is  easily  shown  that 


/.  ( 2 ) 

{Vg}  =  -  ft [Tg]  [ Z ] { XA } 

-  [R]  [Ts]T([Ti|_,]T[TTE)T[Tp]{V}  -  {V} 


where  {V}  = 


[R]  = 


0 

0 

1 


=  +1  for  clockwise  rotating  blades 
=  -1  for  counterclockwise  rotating  blades 


[Z]  = 


0 

1 

0 


-1 

0 

0 


0 

0 

0 


(58) 


as  viewed 
from  the 
shaft  side 


This  establishes  the  relationship  between  the  airspeed  compo¬ 
nents  and  the  system  freedoms.  Equation  20  does ^ the  same  for 
the  blade  twist  0T.  The  blade  segment  velocity  h  is  given  by 


(2) 

h  =  -{L3}T[T^JT[Ts]T(d/dt{XA})  , 


(59) 
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where 


!l3}  »  (0 


* 


cos  p ,  0,  sin  p 

0  ,  1,  0 

-sin  <p ,  0,  cos  <t> 


p  -  tan" 


_i 


A  force  vector  { FgA)  and  a  moment  vector  {M0A}  are  defined  in 
the  normal  blade  section  coordinate  system  from  AL,  AD,  and 
AM  as 


<FBA}  =  0 


{MRa}  =  AM 


These  are  resolved  into  the  fixed-coordinate  system  yielding 

iFF£'  -  1TfHFba} 

{ MpA )  =  [  Tp  ]  { Mg A  }  (  6  < 

where  [Tp]  =  [Tp]  [Tw]  ITTE]  [T^]  [Tg]  .  (6: 

Their  statical  equivalents  at  the  origin  of  the  fixed- 
coordinate  system  are 

{Fp}  =  {FFA} 

{Mp>  =  (Mfa)  +  [Sp] T{ FFA}  (6- 

where  [Sp]  =  [{Sp^},  (Sp2)/  ( Sp3 } ] 

{ SF1 }  =  [ X ] { XFA } 

{SF2>  =  EY] {XFA} 


{ SF3 }  =  [Z]{XFA} 
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i  ,'S,  v,..^ 


..  -v  £vC-  >5 


and  [Xj  = 


[Y]  » 


0 

0 

0 

0 

0 

-1 


0  0 

0  -1 

1  0 


L 


0 

0 

0 


1 

0 

0 


Therefore , 


(JL-IFj,})  =  (JL.[TF]HFBA) 

3<3j  o  3c3j  o 


+  ‘V  (4t{fba!) 

0  ^3 


But  ^  !FBA) 


(9^(FBA,>  '  ^  +  (^(FBA» 


d  Uip 


+  (-^{Fba})  • 


auT 


l3 
9U>j> 

4 

3UT 

3qj 


3Ui 


5Up 

3qj 

dU 


+  ( — S — {  fba}  )  •  — 
3q^  9U„  BA  3q 


3  9Up 

+  <T7T-^BA>)  * 


3q- 


+  (ae^{FBA}) 


39T  (3/f  })  .  9h_ 
9^~  W  BAj;  9q. 


+  ("T{fBA}) 

ah  BA 


3h 

sq. 


This  is  expressed  in  matrix  notation  by  defining  the  3  X 
matrix  [KF]  and  the  nine-component  vector  { A j }  as 


[KF]  =  [  (__£ — { F  BA)  )  ;  ^?~-‘FbA^  ;  *7X“*FBA}  *  ;  ^  A  {tBA}' 


3Un 


aui 


BAJ 


3  Ur 


9  Ur 


(4-1  F*a»;  (”{FBA})  ;  {7?{FBA}); 


3Ut 


auT 


BA 


BA' 


9h 


<^{fba})J 


3h 


(65) 

(66) 


(67) 
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(68) 


A.}T  =  [3U<r-  SUP.  ™T.  ££p.  9Ut.  9Up<  30t  3h  3h 

j  3qjf  3qj'  3c2j'  9c3 j '  sij"'  sq^'  3q~'  3g“;  3q~]  * 


Hence'  JET  {FBA}  =  tKF ) { A j } 
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From  equation  64, 
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Following  equation  69,  3/3qj{MBA)  can  be  written 


as 


3q~  {MBA}  = 


(72) 


where  [KM]  =  [<^{MBA>);  (£•{*%*}),  (4-{MBA});  <4-{MBA»  ; 
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Hence,  (“-{Mp}) 
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The  derivatives  of  {Fp}  and  {Mp}  with  respect  to  and  qj  are 
similarly  obtained  as  J  J 
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STRAIN  ENERGY  AND  DISSIPATION  TERMS 


The  strain  energy  is  a  function  of  the  elasticity  and  the 
state  of  strain  in  the  structural  members  of  the  system.  The 
latter  depends  on  the  system  freedoms  but  not  on  their  time 
derivatives ; 

i  .  e . ,  U  =  U  (qj_ ;  t)  (77) 
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The  strain  energy  term  in  Lagrange's  equations  is  therefore 
represented  by  the  Taylor's  series 


3U _ 3U_  n  ,  3 2U  , 

3qj-  3q  j  0  +  1  (3qj3c3k  o  qk 

k 


+  V2  II 
k  1 


( — — - ) 

9qj3qk3ql  Q 


^k  *1 


+  higher  order  terms.  (78) 

The  only  contribution  to  the  stability  equations  is  in  the 
stiffness  matrix; 


l  .e, 
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SE 


(i»j)  =  ( 
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9qi3q.j 


) 


'J  0 


(79) 


The  dissipation  function  is  assumed  to  be  a  function  of  the 
first  time  derivatives  of  the  system  freedoms, 

R  =  R(qi,t)  .  (80) 


As  such  it  contributes  only  to  the  damping  coefficients  in  the 
stability  equations; 


2  n 

i.e.,  CDISS(i,j)  =  (■■■?■  ft.-)  •  (81) 

3qi 3q j  o 

The  coefficient  matrices  of  the  stability  equations  are  seen 
to  comprise  constituents  from  kinetic  energy,  strain  energy, 
dissipation,  and  external  load  considerations. 


[Mj  = 

[M]Ke 

+ 

Mex 

[C]  = 

t^KE 

+ 

^diss  + 

[K]  = 

+ 

[K]gE  +  tKlEX 

Each  of  these  constituents  is  made  up  of  contributions  from 
the  three  basic  subsystems  as  discussed  previously  under 
kinetic  energy  and  external  loading.  Separating  the  strain 
energy  constituent  into  the  subsystem  contributions, 
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SE 


^KR^SE 
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[KLG]SE 
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The  strain  energy  of  a  single  landing  gear  (see  Figure  1)  is 
given  by 

I AK1 


ULG  =  1/2  {Xfm} 


ak2/ 


ak4 

T 


+  1/2({XFM}  -  ( Xpp } ) 


r 


AK- 


({XFMJ  -  ( Xpp } )  .  (84) 


°  2  T 
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where  £  indicates  summation  over  all  landing  gears  and 
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LGK1J  = 

ak2 

AK3  +  AK4 

tGK2J  = 


o 


AKr 


(86) 


The  airframe  component  is  obtained  from  the  modal  frequencies 
( a>A ( i )  )  and  inertias  (M&(i))  using  the  diagonal  stiffness 
matrix  property  associated  with  normal  modes; 


i.e.,  Ikaf^SE  ^wA(i)  7MA(i)  ^ 


(87) 


In  the  rotor  subsystem,  the  modal  stiffness  is  made  up  of  the 
kinetic  energy  and  strain  energy  constituents: 


[Kr]  KE  +  ^KR^SE 


(88) 
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The  kinetic  energy  part  is  of  the  form 

[KR^KE  =  [kR^ 

The  modes  corresponding  to  a  rotating  blade  (in  vacuo)  at  an 
angular  speed  of  AreF  are  assumed  to  give  rise  to  a  diagonal 
stiffness  matrix  applicable  to  that  speed; 

i.e.,  [kr1ref  =  ['“b  ^  2mb  ^  J 

=  ^2REF  tkR^  +  ^Kr] SE 

Hence,  t^R^SE  =  ~  ^2REF^kR^  (89) 

where  oiB(i)  is  the  ith  rotating  mode  frequency  of  the  blade 
measured  at  fiREF 

Mr (i)  is  the  ith  modal  mass. 

The  constituent  [CJrisS  i-s  separable  into  landing  gear,  air¬ 
frame,  and  rotor  contributions  similar  to  [K] SE  *  It  can  be 
shown  by  consideration  of  the  dissipation  functions  of  the 
landing  gears  that 
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where 
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Structural  damping  in  the  airframe  and  rotor  subsystems  is 
specified  as  fractions  of  the  critical  viscous  damping  coeffi- 


cients 

of  the  various 

modes . 

Hence , 

^AF^DISS  =  [ 

2uA(i)aJA(i)MA(i)] 

tCR]DISS  =  [ 

2pg  ( i)  uig  (i) Mg  (i)  J 

where  pa(i)  =  fraction  of  critical  viscous  damping  in  ith 

airframe  mode 

uB(i)  =  fraction  of  critical  viscous  damping  in  ith 
blade  mode. 

RESULTANT  EQUATIONS  AND  SOLUTION 


The  foregoing  analytical  developments  establish  the  numerical 
and  algebraic  manipulations  involved  in  assembling  the  coeffi¬ 
cients  of  the  stability  equations.  Two  basic  processes  are 
involved.  The  first  is  the  formation  of  a  number  of  trans¬ 
formation  matrices  and  vectors  as  defined  in  coordinate 
systems  and  transformations,  using  basic  input  data  which 
defines  the  system  and  operating  conditions.  The  second 
process  is  execution  of  the  matrix  algebra  indicated  in  this 
and  the  previous  section. 

The  transformations  involving  azimuth  angle  (i/j)  contain 
trigonometric  functions  of  the  time-varying  quantity  (fit)  as 
indicated  in  equation  26.  The  coefficients  of  the  stability 
equations  therefore  contain  these  harmonic  terms  and  are 
arrived  at  by  modifying  standard  matrix  operations  to  include 
harmonic  analyses  concurrently. 

The  solution  of  systems  of  equations  with  harmonic  coeffi¬ 
cients  has  been  the  subject  of  various  publications  in  the 
past  (see,  for  example,  Reference  3).  These  solutions  have 
largely  been  based  on  Floquet's  theorem.  The  present  develop¬ 
ment  uses  a  quasi-normal  mode  approach  which  replaces  individ¬ 
ual  blade  modal  freedoms  with  rotor  freedoms  (equivalent  to 
tip-path-plan^  modal  freedoms  in  flapping  blades) .  The  latter 
are  further  represented  by  a  truncated  general  Fourier  series 
giving  rise  to  the  so-called  collective  (zero-order  harmonic) 
and  cyclic  (first-  and  higher-order  harmonics)  freedoms.  Only 
the  cir=t  harmonic  is  considered  in  the  present  treatment. 
References  4  and  5  contain  more  detailed  discussion  and 
definition  of  quasi-normal  modes. 
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Part  II.  Aeroelastic  Stability  Analysis 
Correlation  and  Trend  Studies 

H.  R.  Alexander 


INTRODUCTION 


AEROELASTIC  STABILITY  CORRELATION 

Hingeless  rotors,  and  more  particularly  the  bingeless  tilt- 
rotor  V/STOL  configuration,  present  a  number  of  technical 
challenges  which  are  not  found  in  an  acute  form  in  conventional 
articulated  rotor  helicopters  or  aircraft  powered  by  relatively 
rigid  propellers  with  high  disc  loadings. 

Three  problem  areas  have  been  discussed  in  detail  in  Volume  I 
and  include  such  phenomena  as  whirl  flutter  and  divergence, 
ground  and  air  resonance,  and  aeromechanical  instability. 

'T’hese  classical  stability  problems  have  the  common  feature 
that  they  may  be  analyzed  and  predicted  in  terms  of  airframe 
degrees  of  freedom  and  a  simple  blade  representation.  Thus, 
whirl  flutter  tends  to  be  associated  with  flapping  motion  of 
the  blade  and  pitching  or  yawing  of  the  hub;  resonance  phenom¬ 
ena  such  as  ground  resonance  and  aeromechanical  instabilities 
may  be  analyzed  in  terms  of  linear  hub  motion  in  the  rotor 
plane  coupling  with  blade  lead-lag  degrees  of  freedom.  How¬ 
ever,  experience  has  shown  that  additional  instabilities  may 
occur  when  the  blades  adopt  deflection  shapes  under  load, 
introducing  strong  coupling  between  pitch  or  torsion  and  lag 
or  flap.  Furthermore,  such  coupling  mechanism  may  alter  the 
blade  behavior  to  the  point  where  a  wing-rotor  instability 
occurs  even  though  the  blade  by  itself  is  stable. 

CORRELATION  WITH  1/10-SCALE  TILT-ROTOR  DYNAMIC  MODEL 

Incidents  involving  this  type  of  coupling  have  generally  been 
of  limited  amplitude.  Correlation  with  two  such  incidents  is 
presented  in  the  next  section.  These  occurred  on  the  Model 
160  1/10-scale  dynamically  similar  tilt-rotor  model  (Figure 
10)  . 

ONERA  13-FOOT-DIAMETER  PROPELLER 

The  blade  itself  may  become  unstable  when  subjected  to  thrust 
and  torque  conditions  causing  deflections  which,  in  turn, 
couple  the  flap,  lag,  and  torsion  degrees  of  freedom.  A  cor¬ 
relation  has  been  made  with  an  incident  of  this  type  which 
occurred  during  a  propeller  performance  test  at  ONERA  (France), 
Figure  11.  The  blade  experienced  a  severe  instability  at  high 
advance  ratio,  which  is  predictable  in  terms  of  blade  proper¬ 
ties  alone.  At  very  low  advance  ratios  a  rather  mild  limit- 
cycle  instability  was  experienced  and  in  this  case  the  blade 
docs  not  go  unstable  by  itself,  but  when  the  hub  degrees  of 
freedom  are  introduced  the  instability  is  obtained. 

MODEL  160  TILT-ROTOR  1/22-SCALE  DYNAMICALLY  SIMILAR  MODEL 

In  Volume  I  it  was  indicated  that  tilting  the  rotor  in  wind 
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tunnel  tests  of  the  Model  160  1/22-scale  dynamically  similar 
tilt-rotor  (Figure  12)  had  produced  significant  and  apparently 
contradictory  shifts  in  the  stability  boundaries.  The  method¬ 
ology  previously  available  did  not  permit  correlation,  but 
application  of  the  computer  program  described  in  Part  I  gives 
the  trends  experienced  in  test. 

ROTOR  DERIVATIVE  TREND  STUDIES  AND  CORRELATION 


Parametric  studies  of  rotor  hub  forces  and  moment  derivatives 
were  performed  as  part  of  the  subject  contract  to  establish 
the  effect  of  selected  parameters  such  as  blade  frequency  and 
Lock  number  on  the  derivatives  and  therefore  on  the  static 
divergence  of  rotor-wing  systems.  These  initial  studies  were 
made  with  a  simple  mathematical  model  (Figure  13)  and  indicated 
unexpected  and  interesting  trends  in  the  rotor  hub  forces  and 
moment  derivatives  when  the  lead-lag  frequency  of  the  blade 
was  parametrically  varied  over  a  range  including  1  per  rev. 
Specifically,  the  normal  force  derivative  was  observed  to 

reach  a  minimum  at  1  per  rev  and  the  pitching  moment  derivative 
Cma  reversed  sign.  In  a  subsequent  test  of  a  dynamically 

similar  1/9-scale  tilt/stowed-rotor  conversion  model  (Figure 
14) ,  the  rotor  derivatives  were  carefully  measured  over  a  range 
of  rpm  in  which  the  lead-lag  mode  frequency  progressed  from 
less  than  1  per  rev  at  900  rpm  to  values  significantly  greater 
than  1  per  rev  as  the  rotor  was  feathered.  The  measured  values 
demonstrated  the  predicted  behavior  trend  and  subsequent  cor¬ 
relation  calculations  using  the  model  properties  and  a  method 
based  on  the  stability  technology  outlined  in  Part  I  gave  ex¬ 
cellent  agreement  with  test. 

STABILITY  TREND  STUDIES 


In  selecting  which  parameters  to  investigate  we  have  been 
guided  by  two  considerations:  (a)  to  use  the  capability  provided 
by  the  new  stability  methodology,  (b)  to  explore  those  param¬ 
eters  which  we  would  intuitively  expect  to  be  critical  or  are 
known  from  test  to  be  those  to  which  instabilities  are  sensi¬ 
tive.  Then  we  have  concentrated  on  instability  mechanisms  in¬ 
volving  blade  pitch,  lag,  and  flap  degrees  of  freedom  because 
the  classical  mechanisms  are  handled  adequately  by  the  previ¬ 
ously  available  technology.  In  these  areas  the  effect  of 
thrust,  precone,  and  drag  offset  have  been  investigated,  along 
with  the  variation  in  blade  torsional  frequency  and  lead-lag 
frequency.  As  points  of  departure  in  these  studies  we  have 
used  existing  scaled  model  parameters  or  preliminary  blade 
design  parameters.  Both  individual  blade  flutter  and  wing- 
rotor  limit-cycle  flutter  have  been  studied.  The  results  ob¬ 
tained  suggest  that  further  extensive  studies  might  be  usefully 
accomplished . 


ROTOR  DERIVATIVE  TREND  STUDIES 


As  mentioned  earlier,  rotor  derivative  correlation  was  pre¬ 
ceded  by  an  investigation  of  the  influence  of  blade  parameters 
such  as  flapping  and  lead-lag  natural  frequencies  and  Lock 
number.  Some  of  these  results  were  included  in  Volume  I  and 
are  repeated  here,  along  with  the  additional  results  for  the 
influence  of  lead-lag  natural  frequency. 

ROTOR  LOADS 


Part  IV  of  this  volume  presents  the  results  of  the  rotor  loads 
correlation  using  the  coupled  loads  analysis  described  in  Part 
III.  The  correlation  is  divided  into  two  main  sections: 
prop/rotor  loads  and  helicopter  loads.  Test  data  for  the 
prop/rotor  correlation  were  obtained  from  wind  tunnel  tests  of 
three  prop/rotor  models.  Photographs  of  the  three  prop/rotor 
wind  tunnel  models  are  shown  in  Figures  14,  15,  and  16.  These 
models  offered  a  wide  range  of  blade  dynamic  characteristics 
as  well  as  a  wide  range  of  flight  conditions  and  reliable  test 
data.  Alternating  blade  root  moments  as  well  as  steady  rotor 
in-plane  moments  were  correlated  with  the  loads  analysis  and 
the  results  are  presented  in  Part  IV.  For  the  helicopter  loads 
correlation,  flight  test  data  from  three  tandem-rotor  helicop¬ 
ters  were  used:  the  CH-47C,  the  CH-47  with  advanced-geometry 
blades,  and  the  Model  347  with  CH-47C  blades.  Test  data  from 
these  aircraft  offered  correlation  in  predominately  unstalled 
and  predominately  stalled  conditions.  Alternating  flap  and 
chord  bending  and  alternating  pitch  link  loads  data  were  used 
for  correlation.  In  addition  to  blade  loads  correlation,  air¬ 
loads  test  data  obtained  from  a  wind  tunnel  test  of  a  full- 
scale  CH-34  rotor  were  correlated  with  airloads  predicted  by 
the  coupled  loads  analysis.  Results  of  the  helicopter  corre¬ 
lation  are  presented  in  Part  IV. 
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Figure  10.  Model  160  1/10-Scale  V/STOL  Tilt-Rotor 
Dynamically  Similar  Model 


,L., 


Figure  11.  NASA  Ames  and  AARL  13-Foot  Fiberglass 
Performance  Model  Rotor 
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Figure  12.  Model  160  1/22-Scale  Dynamically  Similar 
V/STOL  Tilt-Rotor  Model 
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Figure  13.  Mathematical  Model  for  Rotor  Derivative 
Parametric  Trend  Studies 


Figure  14 . 


>1  213  1/9-Scale  Dynamically  Similar 
?0L  Tilt/Stowed-Rotor  Conversion  Model 
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AEROELASTIC  STABILITY:  CORRELATION  AND  TREND  STUDIES 


CORRELATION 

Correlation  With  Model  160  1/10-Scale  Tilt-Rotor 
Dynamically  Similar  Full-Span  Model  Instability- 

The  Model  160  1/10-scale  dynamically  similar  full-span  model 
was  found  on  test  to  have  unstable  regions  as  indicated  in- 
Figure  17.  The  instability,  which  was  limit-cycle  in  type, 
occurred  when  the  collective  pitch  settings  were  reduced  to 
the  negative-thrust  region  in  the  q  and  tilt  conditions  being 
tested.  The  unstable  mechanism  involved  flap,  lag,  and  torsion 
motion  of  the  blades,  and  the  wing  and  rigid-body  pitch  degrees 
of  freedom  were  also  observed  to  participate  actively.  In  the 
zero-  and  negative-thrust  regions  the  blades  are  bent  signifi¬ 
cantly  away  from  the  precone  condition  (5  degrees)  so  that 
strong  inertial  coupling  exists  between  lead-lag  blade  motion 
and  blade  torsion  and  pitch.  The  C-39  program  described  in 
Part  I  shows  successful  correlation  as  indicated  by  the  cal¬ 
culated  flutter  boundaries. 

Correlation  With  the  Model  160  1/10-Scale  Tilt-Rotor 
Dynamically  Similar  Semi-Spain'  Model 

This  model  also  was  unstable  in  a  limit-cycle  mode  involving 
wing  bending  and  torsion,  and  blade  lead-lag,  flap,  and  torsion 
degrees  of  freedom  (Figure  18) .  Although  the  model  frequency 
parameters  are  rather  different  from  those  of  the  later  full- 
span  version,  the  behavior  is  similar.  The  analysis  predicts 
the  onset  of  instability  at  slightly  greater  negative-thrust 
conditions  than  those  actually  experienced.  However,  in  view 
of  the  excellent  correlation  with  the  full-span  version,  it 
seems  logical  to  attribute  this  to  less  precise  knowledge  of 
the  earlier  model  parameters. 

13-Foot-Diameter  Fiberglass  Blade  Instability  Correlation 

During  the  aerodynamic  performance  test  of  a  13-foot-diameter 
prop-rotor  conducted  under  NASA  Contract  NAS2-5025  at  the  ONERA 
facility  at  Modane,  France,  two  instabilities  were  encountered. 
A  violent  instability  occurred  at  a  tunnel  speed  of  643  feet 
per  second  at  980  rpm  and  a  low-growth-rate  instability  was  ex¬ 
perienced  at  60  feet  per  second  and  1,040  rpm.  Correlation 
with  each  of  these  instabilities  is  shown  in  Figure  19.  These 
incidents  were  clearly  sensitive  to  the  thrust  and  blade  de¬ 
flections  combination.  The  blades  were  not  scaled  to  any 
realistic  full-scale  design  and  the  vibratory  modes  have  all 
relatively  high  per-rev  values.  The  modal  deflections  are 
highly  coupled  and  three  predominantly  flexural  modes  occur  at 
lower  frequencies  than  the  blade  torsion  mode.  Thus,  a  minimum 
of  four  blade  modes  is  required  in  any  attempt  at  correlation. 
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In  the  high-speed  case  the  calculated  flutter  boundary  is 
plotted  as  a  function  of  collective  pitch  and  rpm.  At  the 
collective  pitch  of  the  incident,  55  degrees,  the  boundary  is 
at  900  rpm  or  80  rpm  before  the  speed  at  which  the  instability 
was  reported;  that  is  to  say,  the  calculation  is  conservative 
but  the  accuracy  is  within  10  percent.  The  low-speed  incident 
could  not  be  predicted  in  terms  of  blade  properties  alone,  but 
when  the  stand  modal  properties  were  included  in  the  calcula¬ 
tion  it  seems  that  a  very  mild  instability  is  predicted  near 
1,040  rpm,  the  rotor  angular  velocity  at  which  the  incident 
occurred. 

Model  160  1/22-Scale  Dyanmically  Similar  V/STOL  Tilt-Rotor 
xnodel  Correlation  Study 

The  Model  160  1/22-scale  model  consisted  of  a  nacelle  and  wing 
spar  of  variable  stiffness  and  a  3-bladed,  articulated,  flap¬ 
ping  rotor.  Correlation  with  analysis  was  reasonably  good  and 
on  the  conservative  side  for  untilted  axial-flow  cases,  i.e., 
classical  whirl  flutter.  When  the  nacelle  was  tilted  relative 
to  the  wing  on  a  wing  spar  relatively  stiff  in  torsion,  a 
significant  reduction  in  the  flutter  speed  was  experienced  as 
shown  in  Figure  20.  However,  with  a  low  torsionally  stiff 
wing  spar,  an  increase  in  flutter  speed  was  recorded  as  shown 
in  Figure  21  when  the  nacelle  and  rotor  were  tilted.  The 
analysis  available  for  correlation  did  not  address  the  problem 
of  nacelle  tilt  and  nonaxial  flow  through  the  rotor  and  the 
current  study  is  the  first  attempt  to  correlate  with  this 
data.  In  Figure  20  the  calculated  boundaries  corresponding  to 
the  test  condition  are  shown.  It  is  seen  that  the  trend  of  a 
significant  reduction  in  flutter  speed  with  4.5  degrees  tilt 
of  the  nacelle  is  also  in  the  calculated  values;  however,  both 
boundaries  are  conservative  in  relation  to  the  test  data. 

Figure  21  shows  calculated  boundary  as  a  function  of  tilt 
angle  for  the  low  torsionally  stiff  wing  spar.  This  is  seen 
to  follow  the  test  data  closely,  although  the  calculation 
becomes  slightly  unconservative  at  the  higher  tilt  cases. 

Correlation  With  Static  Rotor  Derivatives 

Although  the  C-39  computer  program  does  not  compute  directly 
the  steady  hub  force  and  moment  derivatives,  it  evaluates  and 
prints  out  the  partiax  derivative  of  forces  in  each  of  the 
rotor  variables  and  6  orthogonal  linear  and  angular  hub  dis¬ 
placements.  The  derivatives  are  with  respect  to  unit  ampli¬ 
tude,  rate,  or  acceleration  in  each  of  the  hub  displacement 
rotor  variables.  Proceeding  from  this  information,  it  is 
relatively  straightforward  to  compute  the  rotor  steady  response 
and  hence  rotor  hub  forces  and  moments  with  respect  to  such 
quantities  as  incremental  shaft  angle.  This  provides  a  con¬ 
venient  check  on  the  rotor  technology  in  cases  where  no  insta¬ 
bility  boundaries  are  encountered  in  the  test  region.  In  such 
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cases  the  rotor  hub  derivatives  may  still  be  measured.  In  a 
flexible  blade  the  dynamic  response  of  the  blade  has  to  be  com¬ 
puted  accurately  to  provide  the  correct  hub  forces  and  moments. 
This  correlation  with  static  rotor  derivatives  goes  a  long  way 
toward  substantiating  the  accuracy  of  stability  predictions 
when  the  onset  of  instability  does  not  occur  in  the  test  range. 

Correlation  of  this  sort  is  shown  in  Figure  22.  The  measured 
derivatives  were  obtained  from  the  Model  213  1/9-scale  conver¬ 
sion  model,  Figure  14.  The  calculated  values  show  excellent 
agreement  with  the  test  data.  In  this  rotor  the  first  natural 
frequency  of  the  blade  crosses  1  per  rev  at  500  rpm,  so  that 
the  test  data  confirms  the  unexpected  sensitivity  in  this 
region  to  blade  lead-lag  frequency  mentioned  earlier.  At  85 
feet  per  second  excellent  correlation  is  obtained  when  the 
first  and  second  blade  flexural  modes  are  included  (Figure  22) . 
At  higher  advance  ratios  (141  feet  per  second)  the  influence 
of  blade  torsion  has  to  be  included  to  achieve  a  similar  level 
of  correlation  (Fig’u  a  23). 

TREND  STUDIES 

The  preceding  sections  have  demonstrated  a  substantially 
reliable  capability  to  correlate  with  nonclassical  rotor  in¬ 
stability  behavior.  Selected  parametric  variations  are  now 
developed  to  gain  information  on  the  sensitivity  to  parameters 
expected  to  be  critical. 

Effect  of  Advance  Ratio  and  Thrust  at  Constant  RPM  on 
Individual  Blade  Stability 

In  Figure  24,  calculated  pitch-lag-flap  flutter  boundaries  for 
an  individual  blade  are  shown.  These  show  that,  for  the  par¬ 
ticular  design  chosen  for  investigation,  instability  would 
occur  at  zero  airspeed  and  zero  thrust.  As  advance  ratio  in¬ 
creases  the  amount  of  negative  thrust  that  can  be  tolerated 
also  increases,  but  the  amount  of  positive  thrust  that  can  be 
pulled  without  an  instability  is  reduced  at  the  same  time. 

Whether  this  would  prove  to  be  a  problem  in  practice  would  de¬ 
pend  on  the  thrust-speed  schedule  and  also  on  how  this  behavior 
is  modified  when  the  rotor-to-wing  dynamics  are  included. 

Sensitivity  of  Blade  Stability  to  Lead-Lag  Natural  Frequency 

The  j.ead- lag  frequency  of  the  blade  mentioned  in  the  preceding 
paragraph  was  varied  over  a  range  from  0.4  per  rev  to  1.4  per 
rev  (baseline  value  0.73).  Two  regions  of  flutter  are  identi¬ 
fied  as  shown  in  Figure  25.  Negative  thrust  is  destabilizing 
when  the  lead-lag  frequency  is  above  1  per  rev  and  in  the 
vicinity  of  the  flap  frequency  of  1.22  per  rev.  Positive 
thrust  is  destabilizing  for  lead-lag  frequencies  below  1  per 
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rev,  but  below  0.6  per  rev  large  values  of  negative  thrust 
produce-  static  divergence. 

Sensitivity  of  Individual  Blade  Flutter  to  Blade  Torsional 
Natural  Frequency 

The  torsional  frequency  of  the  blade  discussed  in  the  preceding 
two  paragraphs  was  varied  over  the  range  2.0  to  6.5  per  rev. 

It  is  seen  in  Figure  26  that  in  the  positive-thrust  regime  the 
stability  is  relatively  insensitive  to  torsional  frequency. 

This  is  true  down  to  a  value  of  2.2  per  rev  where  a  pitch-flap 
flutter  boundary  occurs.  In  the  negative-thrust  regime  this 
merges  into  a  pitch-flap  divergence  boundary.  This  suggests 
that  individual  blade  flutter  is  relatively  insensitive  to 
variations  in  blade  torsional  frequency  over  practical  ranges 
of  torsional  stiffness. 

Sensitivity  of  Wing-Rotor  Pitch-Lag-Flap  Instabilities  to 
Precone  Angle  and  Thrust 

For  this  study  data  pertaining  to  the  Model  160  1/10-scale 
dynamically  similar  model  was  used  at  an  advance  ratio  of  0.13. 
It  was  expected  that  significant  changes  would  be  seen  in  the 
amounts  of  positive  and  negative  thrust  required  to  produce 
instability  when  the  precone  angle  is  changed.  This  was  found 
not  to  be  the  case  as  shown  in  Figure  27.  Over  a  wide  range 
of  precone  values  there  is  remarkably  little  change  in  the 
thrust  level  required  for  the  onset  of  the  wing-rotor  insta¬ 
bility.  This  indicates  that  the  blade  deflection  away  from 
the  precone  position,  and  hence  the  flap-pitch  and  lag-pitch 
inertial  coupling,  is  much  less  important  than  the  magnitude 
of  the  steady  aerodynamic  forces  (thrust  and  drag)  on  the  rotor 
blade . 

Sensitivity  of  Wing-Rotor  Pitch-Lag-Flap  Flutter  to  Drag 
Offset  and  Thrust 


The  Model  160  1/10-scale  dynamically  similar  model  was  selected 
for  study.  Since  the  drag  offset  affects  the  blade  elastic 
deflections,  it  was  anticipated  that  a  significant  influence 
on  pitch-lag-flap  flutter  would  be  demonstrated.  However,  a 
result  very  similar  to  that  for  precone  was  evaluated  as  shown 
in  Figure  28.  We  must  conclude,  as  in  the  case  of  precone, 
that  the  inertial  coupling  variations  introduced  by  offset  are 
much  less  significant  than  the  steady  blade  force  modal 
couplings . 
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Figure  17.  Model  160  1/10-Scale  Tilt-Rotor  Dynamically 
Similar  Full-Span  Model  Correlation 
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Figure  18. 
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Figure  19.  Blade  Flutter  Correlation  13'  Diameter  Rotor 
Performance  Test  (ONERA  1968) 
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Figure  20.  Model  160  1/??  o ,  7 

Tilt-Rotor  ModelCCo?relationaTlY  ?imilar  V/STOL 
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Figure  21.  Model  160  1/22-Scale  Dynamically  Similar  V/STOL 
Tilt-Rotor  Model  Correlation  Torsionally  Weak 
Wing  Spar 
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Figure  22.  Model  213  1/9-Scale  Conversion  Model  -  85  Ft/'Sec 
Derivative  Variation  with  RPM 
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Figure  23.  M213  1/9-Scale  Conversion  Model  -  141  Ft/Sec 

Derivative  Variation  With  RPM 


*>**»JaZ£***~  , 


i 


=  0.7  3 


STABLE 


UNSTABLE 


-.16 


UNSTABLE 


Figure  24.  Calculated  Pitch-Lag-Flap  Flutter  Boundarie 
for  a  Preliminary  Design  Blade 


Figure  25.  Blade  Flutter  Sensitivity  to  Changes  in 
Lead-Lag  Frequency  Ratio  and  Thrust 
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Figure  26.  Individual  Blade  Flutter  Sensitivity  to 
Torsional  Natural  Frequency 
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Figure  27.  Sensitivity  of  Wing/Rotor  Pitch-Lag-Flap 
Flutter  to  Precone  Angle  and  Thrust 

69 


rigure  28.  Sensitivity  of  Wing/Rotor  Pitch-Lag-Flap 
Flutter  to  Drag  Offset  and  Thrust 
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ROTOR  DERIVATIVE  TREND  STUDIES 


ROTOR  DERIVATIVE  SENSITIVITY  TO  ADVANCE  RATIO,  FLAP 
FREQUENCY,  AND  LOCK  NUMBER  ~~ 

The  static  derivatives  of  the  rotor  provide  a  convenient  means 
of  determining  the  accuracy  of  aeroelastic  blade  theory.  The 
parametric  behavior  of  the  derivatives  is  important  in  itself 
since  the  static  divergence  characteristics  of  V/STOL  tilt- 
rotors  are  governed  by  the  rotor  normal  force  and  moment 
derivatives.  This  is  discussed  at  some  length  in  Volume  I  and 
need  not  be  repeated  here.  In  view  of  this  importance  of  the 
rotor  static  derivatives,  studies  have  been  accomplished  in¬ 
vestigating  the  effects  of  blade  advance  ratio,  flap  natural 
frequencies,  lag  natural  frequency,  and  Lock  number.  With  the 
exception  of  lead-lag  natural  frequency,  these  are  discussed 
in  Volume  I.  Figures  29  and  30  showing  the  effects  of  advance 
ratio,  blade  flapping  frequency,  and  Lock  number  are  included 
in  this  volume  for  ease  of  reference.  As  is  noted  in  Figure 
29,  the  derivatives  are  extremely  sensitive  to  blade  flapping 
frequency,  particularly  in  the  1.2-  to  1.4-per-rev  range 
typical  of  propeller/rotor  designs  for  V/STOL. 

In  Figure  30  some  results  are  presented  for  Lock  number  vari¬ 
ation.  These  show  a  significant  sensitivity  to  Lock  number, 
particularly  in  the  pitching  moment  derivative  Cma  and  side 
force  derivative  Cya ,  where  a  reversal  of  trend  is  also  noted 
as  the  blade  frequency  increases  from  1.2  to  1.6  per  rev. 

Overall,  the  results  indicate  the  need  for  care  in  the  evalua¬ 
tion  of  blade  parameters  to  be  used  in  static  and  dynamic 
calculations  involving  the  blade. 

ROTOR  DERIVATIVE  SENSITIVITY  TO  LEAD-LAG  NATURAL  FREQUENCY 

In  the  derivative  studies  in  the  preceding  paragraphs  the 
blade  assumed  lead-lag  frequency  was  sufficiently  high  that 
the  freedom  did  not  affect  the  results.  In  Figures  31  and  32 
the  parametric  behavior  with  lead-lag  natural  frequency  is 
presented.  As  in  the  earlier  work,  the  mathematical  model  of 
Figure  13  was  used.  A  dramatic  variation  of  all  derivatives 
is  noted  when  the  lead-lag  frequency  is  in  the  vicinity  of  one 
per  rev.  Particularly  notable  is  the  change  in  sign  of  the 
pitching-moment  and  side-force  derivatives.  These  variations 
are  of  profound  importance  to  the  static  and  dynamic  aero¬ 
elastic  properties  of  rotor-powered  aircraft  whose  designs 
incorporate  blade  in-plane  frequencies  below  or  near  one  per 
rev.  The  discovery  of  this  effect  is  a  significant  new  con¬ 
tribution  to  our  body  of  knowledge  of  rotor  system  behavior. 

An  intuitive  explanation  of  this  sensitivity  to  lead-lag  fre¬ 
quency  may  be  given.  The  phase  of  the  lead-lag  motion  to 
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one-per-rev  aerodynamic  forcing  caused  by  the  shaft  angle  a 
tends  to  increase  rapidly  from  low  value  to  180  degrees  as  the 
blade  passes  through  resonance.  Then  the  coupling  between 
lead- lag  motion  and  flapping  motion  undergoes  a  change  of  sign 
as  the  blade  lead-lag  frequency  traverses  one  per  rev.  In 
addition,  the  lead-lag  amplitude  will  be  magnified  in  this 
region  so  that  the  cross-couplings  overcome  the  natural  ten¬ 
dency  of  the  blade  to  flap  back  (i.e.,  Cmcx  positive)  with 
increasing  shaft  angle. 

As  noted  earlier,  these  trends  have  been  fully  confirmed  by 
subsequent  test  data  and  correlation. 
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Variation  of  Rotor  Derivatives  With  F Lap  Frequency.  Lag  Frequency  High 


Part  III.  New  Methodology  for  Aeroelastic 
Prop/Rotor  Loads  Analysis 


F.  J.  Tarzanin 
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INTRODUCTION 


An  ideally  designed  prop/rotor  aircraft  has  all  its  struc¬ 
tural  load  limits  beyond  the  power-limited  flight  envelope. 
When  structural  loads  limit  the  flight  envelope,  the  aircraft 
is  inefficient  and  costly,  since  it  contains  an  oversized 
powerplant  which  cannot  be  fully  used.  In  addition,  problems 
of  limited  flight  envelope,  limited  component  life,  safety, 
frequent  overhauls,  and  high  maintenance  costs  reduce  the 
effectiveness  of  the  aircraft.  When  the  structure  has  been 
greatly  oversized  and  many  structural  limits  are  significantly 
beyond  the  flight  capability,  a  penalty  is  paid  through  in¬ 
creased  aircraft  weight  and  cost.  Only  when  a  design  includes 
structural  limits,  power  available,  and  flight  envelope  that 
are  matched  can  the  best  aircraft  be  obtained. 

To  obtain  the  ideal  design,  rotor  loads  for  any  proposed 
design  must  be  accurately  known  throughout  the  flight  envel¬ 
ope  for  both  maneuvers  and  level  flight.  To  obtain  these 
loads  requires  a  mixture  of  test  and  analysis.  Analysis  is 
the  most  efficient  approach,  since  parameter  variations, 
flight  envelope  investigations,  and  tradeoff  studies  can  be 
performed  with  a  minimum  of  time  and  cost.  However,  wind 
tunnel  tests  must  also  be  performed  to  check  the  analytical 
results,  establish  confidence,  and  account  for  the  unknowns 
of  actual  hardware  without  the  scatter  of  flight  data.  Some 
flight  data  for  a  similar  aircraft  is  needed  to  determine 
actual  load  scatter  in  the  real  environment. 

The  rotor  loads  that  can  limit  the  aircraft  flight  envelope 
include  blade  loads,  control  loads,  and  hub  loads.  Hub  loads 
directly  cause  fuselage  vibration,  which  can  become  so  severe 
that  specifications  are  exceeded  or  pilot  effectiveness  is 
reduced.  Blade  loads  include  flapwise  and  chordwise  blade 
bending  moments  and  blade  torsional  moments  which  are  used  to 
determine  blade  life.  Control  loads  determine  pitch  link, 
swashplate,  and  control  actuator  life  and  can  affect  aircraft 
controllability.  All  these  loads  are  required  at  critical 
points  of  the  flight  envelope,  including  level  flight,  maneu¬ 
vers,  and  emergency  conditions,  to  effectively  design  the 
rotor . 

However,  rotor  analysis  will  have  a  critical  role  in  the  se¬ 
lection  and  design  of  advanced  prop/rotor  aircraft.  The 
Vertol  Division  of  The  Boeing  Company  has  an  operational 
steady-state  rotor  loads  analysis  that  has  shown  good  agree¬ 
ment  with  measured  helicopter  flight  data.  However,  the 
Vertol  rotor  program  is  intended  to  analyze  helicopter  rotors 
and  must  be  further  developed  to  correctly  analyze  propeller- 
type  rotors . 


For  highly  twisted  blades,  the  principal  axes  in  the  chord- 
wise  and  flapwise  directions  change  their  orientation  relative 
to  the  rotor  shaft  at  each  spanwise  location.  The  helicopter 
analysis  does  not  consider  this  principal  axis  rotation.  To 
accurately  calculate  the  blade  deflections,  including  prin¬ 
cipal  axis  rotation,  requires  that  either  all  calculations  be 
performed  in  a  local-axis  system  coincident  with  the  local 
principal  axis  orientation,  or  in  a  fixed-axis  system  with 
all  the  cross-inertia  terms  included.  In  either  case,  the 
analysis  is  best  performed  if  the  chordwise  and  flapwise  de¬ 
flections  are  coupled.  In  addition,  since  the  chordwise  mass 
center  is  not  coincident  with  the  shear  center,  coupling  with 
torsion  will  also  result.  Therefore,  to  properly  analyze 
rotors  with  large  twist  (propellers) ,  a  coupled  chordwise- 
f lapwise-torsion  analysis  is  required. 

The  present  helicopter  rotor  analysis  assumes  that  the  pitch 
axis,  vertical  neutral  axis,  and  shear  center  are  coincident. 
This  assumption  was  fairly  accurate  when  applied  to  long, 
slender  rotor  blades  with  spar-type  construction.  However, 
recent  developments  in  composite  materials,  new  construction 
techniques,  and  the  larger  chord- to-radius  ratio  of  propellers 
have  made  the  difference  between  the  pitch  axis  and  shear 
center  more  important.  Control  input  (cyclic  and  collective) , 
as  well  as  flexible  control  input  resulting  from  control 
loads,  are  input  about  the  pitch  axis  while  elastic  twist 
resulting  from  blade  torsional  loads  pitches  about  the  shear 
center.  When  the  difference  between  the  pitch  axis  and  shear 
center  is  significant,  both  axis  systems  must  be  taken  into 
account  when  calculating  the  rotor  loads. 

In  light  of  the  need  for  improved  knowledge  of  prop/rotor 
behavior,  contract  number  RFPF  33615-71-R-1075 ,  "Airload 
Predictions  for  the  Avoidance  of  Dynamic  Problems  on  Rotor/ 
Prop  Powered  Military  V/STOL  Aircraft,"  was  awarded.  A  por¬ 
tion  of  this  contract  provided  for  the  continued  development 
of  the  Vertol  rotor  loads  program  (C-60)  to  be  expanded  to 
allow  analysis  of  propellers.  These  improvements  included 
the  effect  of  large  twist,  shear  center  chordwise  location, 
vertical  neutral  axis  chordwise  location,  and  fully  coupled 
blade  deflections.  This  section  outlines  the  prop/rotor 
methodology  and  provides  user's  instructions  for  the  computer 
programs  developed  under  this  contract. 
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HISTORY  OF  DEVELOPMENT 


This  program  has  undergone  continual  modification  and  improve¬ 
ment  since  its  inception.  The  original  analysis,  which  used 
10  mass  and  elastic  bays,  4  harmonics  of  loads  and  deflections, 
uniform  inflow,  and  linear  aerodynamics,  was  developed  by  the 
Vertol  Division  Dynamics  Group  for  their  study  of  helicopter 
rotor  hub  vibratory  forces  under  Navy  Bureau  of  Weapons  Con¬ 
tract  NOa(s)  60-6112c.  The  program  development  and  parameter 
study  results  are  given  in  Reference  6.  Following  the  comple¬ 
tion  of  this  contract,  the  analysis  was  expanded  to  include 
15  mass  and  elastic  bays,  10  harmonics  of  loads  and  deflec¬ 
tions,  and  nonlinear  aerodynamics,  including  stall  and  com¬ 
pressibility  effects.  As  programmed  on  the  IBM  650,  the  actual 
computing  time  of  the  expanded  analysis  was  approximately  2 
hours  and  required  nand  assembly  of  subprogram  inputs.  With 
the  replacement  of  the  IBM  Model  650  computer  by  the  faster 
computers  and  subsequent  reprogramming  into  a  continuous 
chain,  the  unit  running  time  was  reduced  to  approximately  10 
minutes . 

Taking  advantage  of  the  faster  computer,  the  Research  and 
Development  Group  extended  the  rotor  analysis  to  include 
special  input  formats,  aerodynamic  performance  calculations , 
and  an  internal  nonuniform  downwash  calculation.  Since  down- 
wash  is  a  major  source  of  higher  harmonic  rotor  and  hub  loads, 
the  development  of  an  internal  nonuniform  downwash  computation 
significantly  improved  rotor  load  prediction  capability  (7) . 
These  improvements  are  reported  in  the  program  description  and 
users  instructions  (8) . 

In  subsequent  years,  further  improvements  have  been  m£  ■»  in  the 
thrust  routine,  pitch  and  elastic  matrices,  flap  boundary  con¬ 
ditions,  and  downwash  representation.  Special  problems,  such 
as  mechanical  flap-pitch  coupling  (63),  linear  hub  motion,  and 
flap  dampers  can  now  be  analyzed.  The  engineer  has  the  option 
to  perform  either  a  flexible  or  rigid  blade  analysis.  Either 
linear  or  nonlinear  aerodynamics  may  be  employed  and  additional 
options  are  available  to  represent  downwash.  The  program's 
aerodynamic  theory  was  expanded  to  include  an  approximation 
for  the  unsteady  shed  wake  of  an  oscillating  airfoil  (described 
in  Reference  9)  which  significantly  improved  pitch  link  load 
prediction  capability  (10)  .  Any  airfoil  section  may  be  used 
in  the  load  calculations  by  inputting  two-dimensional  static 
airfoil  data  tables.  Only  representative  stall  and  compressi¬ 
bility  effects  were  accounted  for  previously.  In  addition, 
blades  with  variable  airfoil  sections  that  change  thickness 
ratio  or  camber  with  blade  span  can  be  accounted  for  by  de¬ 
fining  up  to  three  airfoil  sections  along  the  blade  and 
linearly  interpolating  between  them.  A  description  of  the 
program  with  all  the  above  changes  has  been  documented  in 
Reference  11.  A  comparison  of  predicted  loads  with  test  is 
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presented  in  Reference  12  and  the  effect  of  linear  hub  motion 
is  given  in  Reference  13. 

Since  1967  the  theory  has  been  further  improved  to  consider 
blade  tab  and  trailing-edge  bends,  an  increased  number  of  mass 
stations  (from  15  to  20) ,  an  increased  number  of  aero  stations 
(from  10  to  15),  viscous  damping  in  the  pitch  links,  angular 
hub  motion  in  the  blade  load  calculations  (14) ,  hub  loads  for 
rotors  with  2  to  9  blades,  and  a  3-point  or  4-point  swashplate 
actuator  system  for  the  lower  control  load  calculations.  In 
addition,  there  were  two  major  improvements:  1  -  the  develop¬ 
ment  of  a  nonuniform  downwash  theory  that  includes  compatibil¬ 
ity  between  vortex  strength,  downwash,  and  the  airload 
distribution  (15),  and  2  -  the  inclusion  of  a  newly  developed 
nonlinear,  unsteady  aerodynamic  theory  (16).  The  downwash 
modifications  improved  blade  loads  (17)  and  probably  hub 
loads.  The  new  unsteady  aerodynamics  provide  a  significant 
improvement  in  rotor  performance  prediction  (18)  and  can  now 
predict  the  sharp  increase  in  control  loads  that  occur  at  the 
high  C<p /a  or  high  airspeeds  associated  with  blade  stall  (16). 


PROGRAM  FLOW  DIAGRAM 


A  flow  diagram  of  this  analysis  is  shown  iu  Figure  33.  The 
solutions  for  the  nonlinear  aerodynamic  loads  and  the  coupled 
flap-lag-pitch  blade  response  are  performed  in  series  and  up 
to  10  iterations  between  the  airloads  and  blade  response  are 
used  to  obtain  the  final  solution.  An  iterative  solution  is 
used  to  account  for  the  nonlinear  coupling  between  the  blade 
deflections  and  airloads  that  result  from  airfoil  stall  and 
compressibility.  Iteration  techniques  are  also  used  to  obtain 
compatibility  between  the  airloads,  downwash,  and  vortex 
strength  and  to  obtain  a  match  with  a  specified  rotor  thrust. 

A  brief  outline  of  the  computer  procedure  is  given  below, 
followed  by  a  description  of  the  methodology  used  in  this 
analysis . 

The  prop/rotor  loads  program  is  started  by  calculating  initial 
deflections  and  defining  boundary  conditions  from  input  (i.e., 
collective  and  cyclic  pitch  and  root  flap  deflections) .  These 
inputs  are  either  known  (as  in  the  case  of  a  model  test  or 
where  these  quantities  were  measured  in  flight)  or  are  ob¬ 
tained  from  an  aerodynamic  trim  analysis.  The  trim  analysis 
calculates  the  rotor  trim  (i.e.,  aircraft  angle  of  attack, 
thrust,  collective  pitch,  cyclic  pitch,  and  blade  root  flap 
angle)  by  considering  aircraft  gross  weight,  center  of  gravity, 
fuselage  drag,  rigid  blade  properties,  quasi-linear  static 
airfoil  characteristics,  nonuniform  downwash  correction 
factors,  forward  speed,  and  rotor  speed  to  determine  the  air¬ 
loads  required  to  maintain  the  free-flight  aircraft  in 
equilibrium  for  a  steady-state  flight  condition. 

Next,  the  rotor-induced  velocities  are  calculated  to  provide 
a  downwash  field  for  each  blade  (see  Figure  34)  .  First  uni¬ 
form  downwash  is  determined  either  from  input  or  a  simple 
calculation.  If  only  uniform  downwash  is  required,  the  pro¬ 
gram  exits  from  the  downwash  routine  and  proceeds  to  the 
airload  routine;  if  nonuniform  downwash  is  required,  a  complex 
iteration  loop  is  initiated.  The  downwash  field  resulting 
from  this  routine  is  used  throughout  the  program  with  no 
updates  or  modifications. 

For  clarity,  when  calculating  tandem-rotor  downwash,  the  two 
rotors  are  designated  as  the  prime  rotor  and  the  secondary 
rotor.  If  the  program  is  calculating  aft  rotor  loads  of  a 
tandem  helicopter,  the  aft  rotor  is  the  prime  rotor  and  the 
forward  rotor  is  the  secondary  rotor.  In  the  case  of  a 
single-rotor  aircraft,  the  main  rotor  is  always  the  prime 
rotor . 

For  a  tandem  helicopter,  the  program  first  calculates  a  com¬ 
patible  self-induced  downwash  for  the  secondary  rotor.  Then, 
the  resulting  vortex  field  of  this  rotor  is  used  to  calculate 
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Figure  33.  Flow  Chart  for  Rotor  Loads  Program 


Figure  34.  Schematic  Diagram  of  Nonuniform  Downwash  Routine 
for  Normal  Conditions 
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the  rotor  interference  downwash  on  the  prime  rotor.  The  prime 
rotor  self-induced  downwash  is  then  calculated  in  the  already- 
defined  rotor  interference  downwash  field.  The  final  prime 
rotor  downwash  field  is  then  the  sum  of  the  rotor  interference 
downwash  plus  the  self-induced  downwash  of  the  prime  rotor. 

For  a  single-rotor  aircraft,  only  the  self-induced  downwash  is 
calculated. 

Calculating  the  self-induced  downwash  requires  an  iteration 
between  airloads,  trailed  vortex  strength,  and  downwash.  This 
iteration  is  performed  up  to  10  times  to  insure  compatibility 
among  these  three  quantities.  As  Figure  16  shows,  each  time 
the  partial  airload  routine  is  entered  an  iteration  between  the 
thrust  routine  and  airloads  may  be  performed.  This  thrust  rou¬ 
tine  iteration  is  performed  when  agreement  with  the  thrust  de¬ 
fined  by  the  trim  program  is  desired.  The  thrust  routine  should 
be  used  when  analyzing  a  specified  flight  condition  where  rotor 
lift  and  propulsive  force  are  defined.  To  obtain  this  constant 
thrust,  every  time  the  airload  routine  is  entered  an  iteration 
process  is  performed  on  the  airloads  by  altering  the  steady  root 
collective  pitch  until  a  collective  angle  is  found  which  corre¬ 
sponds  to  the  required  thrust.  If  no  thrust  match  is  required, 
the  thrust  routine  iteration  is  not  performed. 

A  thrust  match  is  needed  whenever  the  analysis  of  a  specific 
aircraft  flight  condition  is  required.  For  any  steady  state 
flight  condition,  the  rotor  thrust  and  propulsive  force  must 
overcome  the  aircraft  weight  and  drag  (i.e.,  the  equilibrium 
of  the  aircraft  must  be  maintained) .  To  analyze  a  steady  state 
flight  condition,  a  trim  analysis  calculates  the  rotor  thrust 
required  for  equilibrium,  assuming  rigid  blades  and  uniform 
downwash.  In  general,  the  value  of  collective  defined  by  the 
trim  program  will  not  produce  the  same  thrust  with  nonuniform 
downwash  and  a  flexible  blade.  Therefore,  a  thrust  match  must 
be  performed  to  guarantee  that  the  required  thrust  is  main¬ 
tained.  A  t-hrust  match  is  not  performed  when  the  value  of 
collective  must  be  held  constant.  This  is  often  the  case  in  a 
wind  tunnel  test  where  the  collective  pitch  is  held  constant 
during  the  test  and  the  rotor  thrust  is  allowed  to  vary  by 
changing  rotor  shaft  tilt.  The  model,  however,  is  maintained 
in  equilibrium  by  the  forces  exerted  by  the  test  stand;  there 
is  no  guarantee  that  an  aircraft  in  free  flight  would  remain 
in  equilibrium  under  those  conditions. 

After  the  downwash  field  is  calculated,  the  rigid  blade  deflec¬ 
tions  (initial  deflections)  are  used  to  calculate  the  coriolis 
forces.  Then,  the  nonuniform  downwash,  nonlinear  aerodynamic 
coefficients,  and  rigid  blade  deflections  are  combined  to  cal¬ 
culate  the  airloads.  Here  again,  a  thrust  routine  iteration  is 
provided  to  guarantee  a  thrust  match  if  desired.  Following  the 
airload  routine,  the  airloads  and  coriolis  loads  are  harmonical¬ 
ly  analyzed  and  used  as  forcing  functions  to  calculate  the 
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coupled  flap-lag-pitch  dynamic  response  of  the  blade.  Since  the 
forcing  airloads  and  coriolis  loads  are  nonlinear  functions  of 
the  forced  response  (due  to  stall,  compressibility,  and  blade 
shortening) ,  an  iteration  between  aerodynamic  loads  and  the 
blc.  e  dynamic  response  is  required  to  provide  feedback.  To  per¬ 
form  the  iteration,  the  most  recently  calculated  coupled  flap- 
lag-pitch  deflections  are  substituted  back  into  the  coriolis 
and  aerodynamic  analyses,  the  forcing  functions  are  evaluated 
again,  and  the  coupled  blade  response  analysis  repeated.  The 
iteration  is  continued  until  a  specified  number  of  iterations 
(usually  10)  has  been  completed.  The  number  of  iterations 
specified  should  be  sufficient  to  insure  a  converged  solution. 

Next,  blade  radial  forces  are  calculated  by  considering  blade 
shortening  resulting  from  flap  and  lag  deflections.  Pitch  link 
loads  are  calculated  by  determining  the  blade  system  pitching 
moment  and  all  loads  are  resolved  from  the  undeflected  blade 
system  into  the  deflected  blade  system.  Finally,  fixed  and 
rotating  system  hub  and  lower  control  loads  are  evaluated  by 
combining  the  root  shears  and  moments  with  the  system  geometry. 

The  program  flow  diagram  shows  that  the  aeroelastic  rotor 
analysis  is  basically  an  aerodynamic  analysis  coupled  by 
iteration  to  a  dynamic  analysis.  A  discussion  of  the  dynamic 
and  aerodynamic  methodologies  follows. 
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DYNAMIC  ANALYSIS 


BLADE  IDEALIZATION 


The  rotor  blade  is  essentially  a  long,  flexible  beam  with  a 
nonsymmetrical  cross  section,  nonlinear  mass  and  elastic 
properties,  and  five  degrees  of  elastic  freedom.  (The  five 
degrees  of  freedom  are  flapwise  and  chordwise  linear  trans¬ 
lations,  torsion,  flap,  and  lag  rotations.  The  elastic  radial 
deflection  can  be  neglected.)  To  calculate  the  blade's  re¬ 
sponse  by  a  classically  rigorous  method  is  certainly  a  waste 
cf  time  and  money  since  the  blade's  physical  properties  and 
aerodynamic  forcing  loads  are  only  approximations.  In  order 
to  obtain  an  accurate  dynamic  response  at  reasonable  cost,  an 
approximate  blade  analysis  will  be  performed.  Therefore,  a 
simple  blade  idealization  that  lends  itself  to  a  direct- 
solution  method  is  required. 

Figure  35  compares  the  actual  prop/rotor  blade  with  the  blade 
idealization.  The  actual  blade  has  large,  nonlinear  twist, 
variable  shear  center  and  vertical  neutral  axis,  nonlinear 
mass  and  elastic  properties,  variable  planform,  and  variable 
cross  section.  The  idealized  blade  is  obtained  by  cutting  the 
actual  blade,  into  sections,  each  section  idealized  to  have 
uniform  mass  and  geometric  properties.  Then,  the  distribution 
of  physical  properties  for  the  idealized  blade  approximates 
the  actual  blade  distribution  with  a  stepwise  variation  of  the 
blade  properties  (see  Figure  36) . 

The  blade  section  boundaries  (see  Figure  35)  define  the  blade 
idealization.  Each  blade  section  acts  dynamically  as  a  rigid 
body  and  each  is  elastically  connected  to  the  adjacent  sec¬ 
tions  through  equivalent  stiffness.  The  blade  sections  also 
act  as  idealized  aerodynamic  sections,  each  having  constant 
planform  and  cross  section  and  a  constant  distribution  of 
velocity,  angle  of  attack,  and  airloads  along  its  length. 

Since  each  blade  section  is  rigid  with  a  uniform  mass  distri¬ 
bution,  the  blade  section  can  be  replaced  by  pitch  inertia 
about  the  pitch  axis,  flap  and  lag  inertia  about  the  section 
midpoint,  and  a  concentrated  mass  at  the  section  midpoint 
some  distance  from  the  pitch  axis  (at  the  section  center  of 
gravity) .  The  pitch  axis  is  defined  as  a  straight  line  about 
which  an  inelastic  blade  would  pitch  for  an  applied  control 
input. 

An  elastic  bay  is  defined  as  a  weightless  beam  of  constant 
stiffness  between  the  blade  masses.  The  constant  stiffnesses 
are  defined  as  the  equivalent  stiffness  needed  to  match  the 
section  deflection  calculated  from  the  actual  blade  elastic 
properties.  Torsional  elastic  twist  is  calculated  about  the 
shear  center  and  chordwise  bending  is  calculated  about  the 
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vertical  neutral  axis.  Flap  bending  can  be  calculated  any¬ 
where  since  blade  thickness  is  neglected. 

Therefore,  the  blade  is  idealized  as  an  equivalent  system 
composed  of  20  discrete  mass  points  possessing  flap,  chord, 
and  torsional  moments  of  inertia,  connected  by  20  weightless 
uniformly  elastic  beams.  Associated  with  each  mass  is  a 
rigid  airfoil  segment.  The  airloads  generated  by  each  airfoil 
segment  are  applied  on  the  j;itch  axis  at  the  middle  of  the 
segment. 

DYNAMIC  BLADE  RESPONSE 


Background 

The  aim  of  the  prop/rotor  loads  analysis  is  to  predict  the 
blade  and  hub  loads  for  aircraft  in  steady  flight.  The  prob¬ 
lem  of  predicting  transient  maneuver  loads  must  wait  until  the 
steady-state  loads  have  been  obtained. 

The  steady-state  blade  response  may  be  calculated  by  two 
methods.  One  approach  uses  a  timewise  extrapolation,  which 
uses  the  loads  and  deflections  at  a  given  time  to  extrapolate 
to  loads  and  deflections  at  a  future  time.  This  approach  is 
mandatory  for  transient  load  and  deflection  calculations  but 
is  extremely  time-consuming  for  steady-state  calculations, 
since  the  program  must  start  from  arbitrary  initial  conditions 
and  continues  until  the  response  is  periodic. 

The  method  used  for  this  prop/rotor  program  is  quicker  and 
more  direct.  For  this  approach  the  solution  is  assumed  to  be 
steady-state  and  hence  periodic.  This  allows  the  blade  re¬ 
sponse  to  be  expanded  in  a  series  and  the  solution  is  obtained 
by  finding  the  series  coefficients.  Since  a  rotor  has  a  nat¬ 
ural  period  of  one  rotor  cycle  and  since  virtually  all  steady- 
state  rotor  responses  are  multiples  of  this  natural  period, 
the  blade  response  will  be  represented  by  a  Fourier  series 
based  upon  multiple  frequencies  of  the  rotor  speed.  Each 
multiple  of  the  rotor  speed  is  defined  as  a  separate  harmonic. 
A  deflection  that  occurs  at  the  period  of  the  rotor  frequency 
is  a  first-harmonic  deflection.  A  load  that  experiences  three 
cycles  in  the  period  of  a  rotor  cycle  is  a  third-harmonic 
load . 

A  possible  exception  to  the  Fourier  series  solution  based  on 
rotor  speed  is  the  subharmonic  blade  response  which  can  arise 
from  coupling  and  beating  of  the  rotor  with  other  aircraft 
components,  such  as  a  wing  or  a  feedback  stability  system. 
However,  for  this  analysis  only  the  rotor  will  be  considered 
and  subharmonic  responses  will  be  neglected. 
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This  method  of  solution  will  be  further  simplified  by  neglect¬ 
ing  interharmonic  coupling  of  the  blade  response.  Simply 
stated,  this  assumption  means  that  moments  resulting  from 
terms  like  second-harmonic  deflection  times  third-harmonic 
load  are  neglected.  In  fact,  only  steady  deflections  times 
steady  loads,  steady  deflections  times  harmonic  loads,  and 
harmonic  deflections  times  steady  loads  will  be  considered. 
This  assumption  allows  each  harmonic  to  be  solved  separately, 
requiring  solution  of  the  load  and  deflection  coefficients 
for  one  harmonic  at  a  time.  Therefore,  first  the  steady 
solution  is  found,  then  the  first  harmonic,  second  harmonic, 
etc.  The  final  solution  is  obtained  as  a  continuous  function 
by  reconstituting  the  series  at  24  discrete  positions  around 
the  azimuth. 

To  develop  the  equations  for  the  dynamic  analysis,  a  mass 
coordinate  system  was  chosen  and  the  following  assumptions 
were  made : 

1.  The  blade  idealization  described  above  adequately 
represents  the  blade. 

2.  The  airloads  and  coriolis  loads  are  assumed  to  be  known, 
continuous,  and  periodic. 

3.  The  solution  is  steady-state  and  expandable  in  a  Fourier 
series . 

4.  All  angular  deflections  are  small. 

5.  Steady  centrifugal  untwist  can  be  calculated  once  and 
superimposed  on  the  blade  built-in  twist.  Live  centrif¬ 
ugal  untwist  is  accounted  for  by  increasing  the  effective 
GJ. 

6.  Only  linear  terms  will  be  considered  in  the  dynamic 
response.  Large  nonlinear  deflections  such  as  radial 
deflection  of  the  mass  due  to  flap  and  lag  and  incre¬ 
mental  flap,  lag,  and  pitch  deflections  of  the  shear 
center  are  calculated  separately. 

The  blade  is  analyzed  by  transforming  the  lumped  mass  and 
elastic  elements  of  the  idealized  blade  into  a  sequence  of 
transfer  matrix  products  by  means  of  the  associated-matrix 
method.  This  method  replaces  each  blade  element  with  an 
equivalent  transfer  matrix  that  transfers  the  dynamic  system 
variables — shear,  moment,  deflection,  and  slope — inboard 
across  the  element. 

Since  the  blade  response  is  expanded  in  a  Fourier  series  and 
interharmonic  coupling  is  neglected,  the  blade  response  for 
each  harmonic  is  independent  of  the  others  and  can  be 
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calculated  separately.  Since  this  program  considers  a  steady 
plus  10  harmonics  of  blade  response,  11  dynamic  matrix  equa¬ 
tions  must  be  solved,  one  for  the  steady  and  each  harmonic. 
Each  dynamic  equation  considers  cosine  and  sine  components  of 
coupled  flap-lag-pitch  blade  deflections  and  loads.  The  un¬ 
known  loads  for  each  harmonic  are  shear  normal  to  the  chord, 
shear  in  the  chord  direction,  pitching  moment,  flap  bending, 
and  chord  bending.  The  unknown  deflections  for  each  harmonic 
are  flapwise  and  chordwise  translation,  flapwise  and  chord- 
wise  rotation,  and  pitch  rotation.  Therefore,  the  transfer 
matrix  represents  the  transfer  of  harmonic  coefficients  for 
a  single  harmonic  from  one  point  on  the  blade  to  the  next. 

Transfer  Matrices 

For  this  method,  the  blade  variables  (the  unknown  loads  and 
deflections)  at  each  discrete  point  along  the  blade  are  listed 
in  a  column  matrix  called  a  state  vector.  If  q  is  a  state 
vector  representation  of  the  problem  variables,  [M]  is  the 
transfer  matrix  across  discrete  mass,  [F]  is  the  transfer 
matrix  across  a  discrete  load,  and  [E]  is  the  transfer  matrix 
across  an  elastic  beam,  then  the  matrix  transfer  equations 
can  be  written. 

For  a  discrete  mass  at  station  n, 

HO- 

{qMn  {^}n 


{q}  is  the  vector  of  loads  and  deflections  outboard  of  the 
mass  and  (q'}n  is  the  vector  of  loads  and  deflections  inboard 
of  the  mass.  The  matrix  transfer  equation  is 

{ q ' ) n  =  M n{q}n  • 

For  a  discrete  load  (forcing  function)  at  station  n, 

_ ,  t 

F 

n 

{ q }  { q  1  } 
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{q')n  is  the  vector  of  loads  and  deflections  outboard  of  the 
load  and  {q>n  is  the  vector  of  loads  and  deflections  inboard 
of  the  load.  The  matrix  transfer  equation  is 

{q>n  =  tF]n{q’ } 

To  obtain  the  discrete  loads  for  the  transfer  matrix,  the 
coriolis  force,  aerodynamic  pitching  moment,  lift,  and  drag 
are  harmonically  analyzed  and  the  appropriate  harmonic  compo¬ 
nents  are  inserted  into  the  transfer  matrix  to  represent  a 
concentrated  external  load.  To  help  the  program  converge, 
linear  aerodynamic  damping  and  coupling  terms  are  provided  in 
the  load  transfer  matrix  to  approximate  anticipated  airload 
changes  resulting  from  elastic  deflections.  Substituting 
{q'}n  into  the  above  equation, 

{q>n  =  [F]n[M]n{q} 

For  the  elastic  beam  element  of  length  in, 


{<3Vl 


<q) 


n 


{q}n  is  the  vector  of  loads  and  deflections  outboard  of  the 
elastic  beam  and  (qln+1  is  the  vector  of  loads  and  deflections 
inboard  of  the  elastic  beam.  Then: 

{q}n+l  =  fE]n{^}n  ' 

or,  substituting  for  {q}R  gives: 

^q^n+l  =  fE]j1[F3n[M]n{q}ri  , 

which  is  a  transfer  equation  relating  a  vector  of  the  loads 
and  deflections  acting  at  one  mass  station  to  the  loads  and 
deflections  at  the  next  most  inboard  mass  station. 

Therefore,  multiplying  the  state  vector  of  a  given  harmonic 
by  the  transfer  matrix  gives  the  harmonic  state  vector  inboard 
of  the  blade  element.  Similarly,  for  the  next  inboard  station 
(n+i)  the  transfer  equation  becomes: 

n+1  ^  n+1  ^  n+l^^n+l  • 

[EW[F!n+llM1n+liEln!Fl„tM1„<‘J,n  • 


(9>n+2  = 
or  = 
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.oniTlwtiriFfn’-rt' 


As  shown  here,  the  state  vectors  can  be  premultiplied  by  a 
successive  string  of  transfer  matrices,  each  representing  the 
change  in  loads  and  deflections  across  the  mass,  airload,  and 
elastic  sections  it  represents.  A  boundary  matrix  (or 
matrices)  is  added  between  the  most  inboard  mass  and  the  blade 
hub  to  account  for  multiple  load  paths  (such  as  pitch  links 
and  lag  dampers)  and  hub  geometry  (such  as  hinge  configuration 
and  delta  three) .  In  this  way  the  root  state  vector  can  be 
equated  to  the  tip  state  vector  through  a  string  of  transfer 
matrices  and  the  elastic  blade  idealization  is  then  reduced 
to  a  sequence  of  transfer  matrix  products  (see  Figure  37)  . 
Collapsing  the  string  of  transfer  matrices  (i.e.,  multiplying 
the  matrices  out) ,  a  single  collapsed  transfer  matrix  is  ob¬ 
tained  and  the  blade  equation  for  each  harmonic  becomes : 

(S}root  =  [T .M.  ]  { S  )<pip  • 

where  {S}R00t  =  <VZ(KC)'  MY(KC)'  B (KC) '  ^(KC)'  VY(KC)' 


MZ (KC)  ' 

5  (KC)  ' 

Y  (KC)  '  mX(KC)  ' 

0f (KC) ' 

VZ (KS)  ' 

MY  (KS) 

'  6 (KS) '  Z (KS) ' 

VY(KS) ' 

MZ (KS)  ' 

5  (KS)  ' 

Y(KS) '  MX(KS) ' 

0f (KS) 1 

{ S }ipip  =  a  column  of  the  same  variables  at  (S}R0OT 
but  defined  at  the  blade  tip. 

[T.M.]  =  K  R  E  (20)  A  T19  E(19)B  M19  E(19)A  T1B  E(18)B 

m1B  e  (18)  A  t17  e  ( 17 )  B  m17  e(17)A - E(4)A 

t3  e ( 3 ) B  b18  e (18) A  E (3 ) A  t2  e(2)B  f2  m2 
E  (2)  A  T1  e(1)B  f1  m1 

=  a  21  x  21  matrix  obtained  from  multiplying 
the  indicated  transfer  matrices 

Vg  =  shear  normal  to  the  chord 

My  =  flap  bending  moment 

8  =  flap  angle 

Z  =  translation  normal  to  the  chord 

Vy  =  shear  in  the  chord  direction 
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Mz  =  chord  bending  moment 
£  =  lag  angle 

Y  =  translation  in  the  chord  direction 
=  torsion  moment  in  the  disc  plane 
6f  =  flexible  torsion  rotation 

The  subscript  (KC)  indicates  the  harmonic  cosine  component 
and  (KS)  indicates  the  sine  component. 

M  =  mass  matrix 

F  =  aerodynamic  matrix,  a  discrete  load  matrix 
E  =  elastic  matrix 
T  =  twist  matrix 
K  =  control  system  matrix 

R  =  is  composed  of  a  sequence  of  additional 
matrices  which  vary  with  the  boundary 
configuration  and  hinge  sequence  (for 
articulated  rotors) . 


The  detailed  definitions  of  each  of  the  transfer  matrices  are 
defined  in  the  supplemental  sections  of  this  report. 

Solution  of  Dynamic  Equations 

Since  there  is  a  dynamic  equation  for  each  harmonic,  the  fre¬ 
quency  of  vibration  is  known.  Therefore,  the  velocity  and 
acceleration  terms  can  be  written  as: 


X  =  d/dt(X^c  cos  +  XKS  sin  Ky) 

=  -wK  Xj^-,  sin  Kip  wK  X^g  cos  K<p 

X  =  d’2/dt2  (Xv_  cos  +  X„_  sin  Kip ) 

KC  Kb 

=  ~co2K2  X-^q  cos  KiJ;  -  oj2K?  XKg  sin  Ki|i 

where  u  is  the  rotational  frequency  of  the  rotor  and  K  is  the 
harmonic  number. 


Since  the  time  derivatives  are  removed  from  the  equation  and 
replaced  by  parameters  dependent  on  w,  the  dynamic  equation 
above  reduces  to  a  simple  set  of  simultaneous  algebraic 
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equations.  To  solve  these  equations,  boundary  conditions 
must  be  applied. 


All 

the 

loads  at  the 
tip  becomes: 

blade 

tip  are  zero,  so 

the  state 

vector  at 

{S}TIP  = 

rS 

o 

■'w 

o 

'•i 

6 (KC) '  Z (KC)  '  °' 

O 

vT^f 

n 

Y(KC)  '  °' 

6f (KC) 

t  If  0/  0/  3  ^ )  * 

Z (KS)  '  0  ' 

0 '  £ (KS)  ' 

Y(KS) '  °'  ef (KS) } 


Multiplying  the  zeros  through  the  collapsed  transfer  matrix 
reduces  the  equation  to 

{S}ROOT  =  {S)TIp 

21x1  21x11  llxl  v  ' 

where  {S}TIP  =  {3(KC),  Z (KC) ,  5 (KC)  ,  Y(KC),  6f (Kcj  ,  1,  3 (KC) 

Z (KS)  ’  ?(KS)'  ^  (KS)  '  ®  f  (KS ) ^ 

Next,  apply  the  10  known  root  boundary  conditions  and  select 
the  resulting  10  equations  that  contain  the  known  root  end 
conditions  and  the  10  unknown  tip  deflections.  Invert  the 
resulting  reduced  collapsed  transfer  matrix  and  solve  the 
equations  for  the  coupled  flap-lag-pitch  tip  deflections.  By 
progressively  applying  the  transfer  matrices  to  the  tip  de¬ 
flections,  the  deflections  and  loads  at  each  blade  element  are 
determined.  Repeating  the  procedure  for  each  harmonic  and 
then  superimposing  the  separate  harmonic  responses,  the  total 
coupled  flap-lag-pitch  coupled  loads  and  deflections  are 
found . 

An  example  of  this  procedure  will  be  illustrated  below  for  a 
hingeless  rotor  boundary  condition.  The  root  end  conditions 
for  a  hingeless  rotor  are  zero  flapwise  and  chordwise  deflec¬ 
tions,  defined  flap  and  lag  angles  for  the  zero  harmonic 
(precone  and  prelag)  and  zero  for  all  other  harmonics,  and 
pitch  deflections  defined  by  the  control  input.  Therefore, 

3,  Z,  £,  Y,  and  9  are  known.  Selecting  the  10  equations  foi 
the  known  root  end  conditions  (i.e.,  3kc  equation,  B^s  equa¬ 
tion,  ZK£  equation,  etc.)  from  the  21  equations  of  matrix 
equation  92  gives 

{ sB }  =  (tTm.J  {S}TIP 

10x1  10x11  llxl  (94) 
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where  {Sg}  =  { 0 (KC) /  Z(KC)/  C (KC) ,  5(KC).  6f(KC) 
3(KS)'  Z(KS)'  C(KS)'  Y(KS)'  8 f (KS) 

=  { 0 ,  0  ,  0 ,  0 ,  6g  (kc)  , 

0,  0 ,  0/  0/  ®B(KS)^ 


eB(KC).  anc*  0B(KS)  are  t^ie  known_control  inputs.  Since  (S)tip 
contains  a  unity  element  and  [T7M.]  is  not  square,  equation 
93  can  be  partitioned  and  written  as: 


Rearranging  the  equation  and  solving  for  the  tip  unknowns 
gives : 


Form  the  complete  tip  state  vector  and  premultiply  it  by 
e(2)A  Ti  E(i)g  F]_  Mi  to  obtain  the  state  vector  before  the 
second  mass.  Repeat  this  process  to  obtain  the  state  vector 
before  each  mass  station.  The  state  vectors  for  each  mass 
station  and  each  harmonic,  along  with  the  harmonically  recon¬ 
stituted  loads  and  deflections,  are  printed  out. 
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AERODYNAMIC  ANALYSIS 


AIRLOADS 


This  program's  airload  calculations  include  the  effects  of 
airfoil  section  geometry,  compressibility,  stall,  three- 
dimensional  flow,  unsteady  aerodynamics,  and  nonuniform 
inflow.  Static  airfoil  tables  are  used  to  account  for  com¬ 
pressibility,  static  stall,  and  airfoil  shape.  The  unsteady 
aerodynamic  loads  are  calculated  by  modifying  the  static  loads 
resulting  from  the  airfoil  tables  to  include  Theodorsen's  shed 
wake  function,  dynamic  stall  effects  based  on  oscillating  air¬ 
foil  data,  and  yawed  flow  across  the  blade. 

In  the  development  of  this  theory  it  is  assumed  that  the 
separate  effects  of  stall,  compressibility,  dynamic  stall 
hysteresis,  and  yawed  flow  can  be  superimposed  to  give  a 
theory  that  adequately  defines  the  unsteady,  three-dimensional 
environment  of  helicopter  rotors.  Though  this  assumption  has 
not  been  substantiated  by  specific  tests,  it  has  been  used  to 
demonstrate  a  significant  improvement  in  rotor  performance 
prediction  in  the  rotor  stall  region  (18)  and  control  load 
prediction  (16) . 

DEVELOPMENT  OF  r  IE  THEORY 

Theodorsen  in  Reference  19  developed  an  unsteady  aerodynamic 
theory  for  airfoils  pitching  and  plunging  about  a  zero  mean 
angle  of  attack.  The  theory  assumed  linear  incompressible 
aerodynamics  and  accounted  for  a  sinusoidal  shed  wake  trail¬ 
ing  from  the  airfoil.  To  update  this  theory  to  account  for 
the  helicopter  aerodynamic  environment,  compressibility  and 
blade  stall  were  included  and  the  restrictions  of  small 
angles  and  zero  mean  angle  of  attack  were  removed. 

In  the  lift  equation  this  was  accomplished  by  replacing  the 
linear  aerodynamic  expression  not  by  the  static  aerodynamic 
lift  coefficient  Cl  as  obtained  from  a  table  involving  the 
angle  of  attack  and  Mach  number.  For  the  aerodynamic  pitch¬ 
ing  moment,  the  static  aerodynamic  moment  about  the  quarter 
chord  (as  given  by  the  coef ticient)  was  added  to  the  cal¬ 
culation.  The  Cm  term  was  added  since  the  moment  equation  in 
Reference  19  does  not  account  for  a  static  moment  about  the 
quarter  chord. 

To  use  the  shed  wake  representation  of  Reference  19,  the  dif¬ 
ferences  between  the  fixed-wing  aerodynamic  environment  and 
the  helicopter  environment  must  be  reconciled  by  making  three 
assumptions.  First,  the  curvature  of  the  helicopter's  shed 
wake  pattern  is  neglected  and  a  pattern  parallel  to  the  blade 
is  assumed  (see  Figure  38) .  This  assumption  is  reasonable, 
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•PARALLEL  WAKE 


ACTUAL  SHED  WAKE 
PATTERN 


PARALLEL  WAKE 
ASSUMPTION 


•LINEAR  RELATION  BETWEEN  VORTEX  STRENGTH  AND 
ANGLE  OP  ATTACK 


SHED 

VORTEX 

STRENGTH 


•REDUCED  FREQUENCY  USES  THE  INSTANEOUS 
RELATIVE  VELOCITY  AND  ASSUMES  ONLY  1/REV 
BLADE  MOTION  IS  PRESENT 


Figure  38.  Shed  Wake  Assumptions 
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since  the  most  important  vortices  are  close  to  the  blade 
where  the  effects  of  curvature  are  the  smallest. 

Second,  a  linear  relationship  between  angle  or  attack  and 
shed  wake  is  assumed.  This  assumption  is  valid  in  the  sub¬ 
stall  region  but  incorrect  in  the  stall  region.  However,  the 
error  introduced  is  minimal  since  the  dominant  characteris¬ 
tics  of  the  airloads  in  the  stall  region  are  determined  by 
stall  hysteresis.  The  third  requirement  is  to  assume  an  in¬ 
stantaneous  value  of  reduced  frequency  based  on  the  one-per- 
revolution  flap  and  pitching  motions  and  the  instantaneous 
relative  wind  velocity. 

The  expression  for  this  pseudo-reduced  frequency  is 

K  =  CR/2U  ,  (95) 

where  U  is  the  total,  instantaneous,  two-dimensional  velocity 
relative  to  the  airfoil  and  Q  is  the  one-per-revolution  rota¬ 
tional  frequency. 

To  include  yawed  flow  (i.e.,  radial  flow)  in  the  theory,  an 
empirical  approximation  to  account  for  three-dimensional  flow 
on  the  lift  coefficient  was  developed  in  Reference  20.  The 
equation  used  to  define  the  three-dimensional  Cl  is  given  as: 

Cl ( 3D)  =  Cl (2D) / cos  A  ,  (96) 

where  A  is  the  yaw  angle. 

The  equation  has  the  restriction  that  dCL(3D)A*a  maY  not  ex" 
ceed  the  maximum  value  of  dCL(2D)/da*  This  means  that  the 
effect  of  yawed  flow  is  to  delay  stall  to  increasingly  larger 
values  of  Cl  with  increasingly  larger  yawed  angles.  A  compar¬ 
ison  of  equation  95  with  test  data  obtained  from  Reference  21 
is  shown  in  Figure  39. 

Therefore,  Theodorsen's  theory  has  been  extended  to  provide 
an  unsteady  aerodynamic  theory  that  accounts  for  compressibil¬ 
ity  and  stall  through  static  Cl*  Cq,  and  Cj^  tables;  includes 
unsteady  shed  wake  induced  angle  of  attack  variations  by  using 
the  Reference  19  F  and  G  functions;  and  provides  for  three- 
dimensional  flow  with  an  empirical  equation.  However,  this 
theory  does  not  account  for  the  unsteady  stall  effects  shown 
in  References  22,  23,  and  24  and  illustrated  in  Figure  40. 

To  develop  a  theory  to  include  the  unsteady  stall  effects  ob¬ 
served  above,  analytical  investigations,  tests  of  oscillating 
airfoils,  reviews  of  existing  flight  test  data,  and  extensive 
wind  tunnel  testing  of  model  rotors  were  carried  out  over  a 
3-year  period.  Evaluation  of  the  oscillating  airfoil  test 
results  has  led  to  the  development  of  an  empirical 
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CL ( 3D)  CL ( 2D) /cos  A 

A  =  YAW  ANGLE 


Figure  39.  Engineering  Approximation  Accounting  for  Yawed 
Flow 
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Figure  40.  Typical  Oscillating  Airfoil  Test  Data  for  Vertol 
23010-1.58  Airfoil 
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relationship  between  the  static  stall  angle  and  the  dynamic 
stall  angle  for  lift  and  pitching  moment.  The  development  of 
this  relationship  is  shown  in  Reference  25  and  the  resulting 
equation  is  given  below: 

aD.S.  "  aS.S.  =  YtCa/2U]1/2  •  (97) 

This  relationship  is  used  as  a  bridge  between  static  and 
dynamic  stall  by  defining  the  dynamic  stall  angle  of  attack 
as  a  function  of  the  static  stall  angle  of  attack.  Therefore, 
with  this  relationship  the  dynamic  stall  can  be  derived  from 
two-dimensional  static  airfoil  data.  Since  dynamic  and  static 
stall  can  be  related,  a  pseudo-dynamic  angle  of  attack  can  be 
defined  from  the  quasi-static  angle  of  attack  in  the  same 
manner.  The  dynamic  angle  of  attack  (cxdyn)  can  then  be  ex¬ 
pressed  as  a  combination  of  the  quasi- static  blade  element 
angle  of  attack  (aBE)  and  its  time  derivative  daBE/dt  as: 

aDYN  =  “BE  ~  Y  (C  (dctgg/dt)  /2V]  ,  (98) 

where  the  quasi-static  angle  of  attack  is  defined  in  the  con¬ 
ventional  manner  as 


BE 


=  9  + 


tan-1  UD/U n 


(99) 


and  the  parameter  y  is  derived  from  the  oscillating  airfoil 
test  points  chat  exhibit  dynamic  stall  and  is  experimentally 
dependent  only  on  Mach  number  and  airfoil  shape  (see  Reference 
25) .  The  dynamic  angle  of  attack  can  now  be  used  to  calculate 
a  dynamic  lift  coefficient. 

The  static  lift  coefficient,  CT  ,  can  be  defined  as 

CL  =  (cl/aBE)  “BE  '  (100) 


where  Cl/«be  i s  the  quasi-static  lift  slope.  The  dynamic  lift 
coefficient  can  be  defined  using  a  similar  equation  as 


CL (DYN)  ~  (CL 


(ATaDYN)/“DYN 


“BE 


(  101) 


where  the  quasi-static  lift  slope  is  replaced  by  a  dynamic 
lift  slope  composed  of  the  pseudo-dynamic  angle  of  attack  kdyn 
and  a  value  of  CE  obtained  from  the  static  tables  using  the 
dynamic  angle.  With  this  method  CE  will  not  show  the  effects 
of  stall  until  cxdyn  reaches  stall,  even  though  agE  may  be 
significantly  beyond  the  static  stall  angle. 
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In  any  region  of  the  two-dimensional  static  airfoil  data  where 
the  dynamic  lift  slope  equals  the  quasi-static  lift  slope, 
i  .e . , 


Cl/aBE  "  CL  ^ataDYF^aDYN 


(102) 


this  technique  for  generating  a  lift  hysteresis  is  completely 
removed.  This  occurs  since  the  region  where  equation  99  holds 
is  the  linear  lift  portion  of  the  CL-versus-a  curve.  An  ex¬ 
ample  of  a  dynamic  Cl  calculation  using  this  technique  is 
shown  in  Figure  41  for  an  airfoil  at  a  mean  angle  of  attack 
of  14  degrees  oscillating  with  a  5-degree  amplitude. 

A  different  approach  is  used  to  calculate  the  dynamic  Cm* 
Quasi-statically ,  Cj^  is  obtained  by  looking  up  its  value  from 
a  static  CM-versus-a  curve  using  the  quasi-static  angle  of 
attack  cxbe-  Therefore,  dynamic  Cm  will  be  obtained  by  looking 
up  its  value  using  the  dynamic  angle  of  attack  a£>YN  •  Note 
that  the  static  stall  and  dynamic  stall  angles  for  Cl  and  Cm 
are  not  the  same,  and  hence  the  pseudo-dynamic  angles  also 
differ.  The  subscript  L  and  M  will  be  used  henceforth  with 
the  dynamic  angle  of  attack  to  indicate  whether  the  dynamic 
angle  of  attack  was  obtained  from  the  lift  or  moment  data. 

With  this  approach,  the  expression  for  dynamic  Cm  is 

CM(DYN)  =  CM (ataDYN (M) )  *  (103) 


The  effect  of  this  equation  is  shown  in  Figure  42  for  an  air¬ 
foil  with  a  mean  angle  of  attack  of  14  degrees  and  an  ampli¬ 
tude  of  oscillation  of  5  degrees. 


Superimposing  the  separate  effects  of  compressibility,  linear 
shed  wake,  yawed  flow,  and  dynamic  stall  hysteresis  gives  the 
following  equations  for  dynamic  CL  and  CM: 


-  C 


CL (DYN)  ~ 


CLO  +  “  a 


(at  aDYN (L) J 


LO 


DYN 


cc^  A 


equ 


^Lfat  a  ) 

lat  dyn  (l)  , 

+  -  ^DEL  f  (H,  0) 

daDYN  (L)  DMj 


,  (104) 
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Cl(DYN) 


CL  (AT  a  DYN) 

a  DYN  “B.E. 


Figure  41.  The  Dynamic  CL  Calculation  in  Stall  Is  Dependent 
on  the  Dynamic  Stall  Lelay 


CM(DYN)  “  CM[AT  a  DYN (M) ] 


Figure  42.  The  Dynamic  in  Stall  Is  Dependent  on  the 
Dynamic  Stall  Delay 
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CM(DYN) 


M(at  ctDyN(M)  ) 
+  CQ (PA  -  1/4) 
+  Cq (PA  -  1/4) 
dCT 

(at  “DYH1 
d 


CL (DYN)  cos  “BE 
CD(DYN)  sin  “BE  “  -1/2  C0 


“DEL  cos  “BE  +  F(H/ 


e)  . 

(105) 


The  last  terms  in  the  and  CM  equation  can  be  neglected 
since  they  are  generally  small. 

The  equivalent  angle  of  attack  aequ  is  based  upon  the  mechan¬ 
ical  angle  of  attack,  flapping  motion,  rotor  inflow  from  shaft 
tilt,  and  downwash.  All  these  terms  give  a  relative  velocity 
that  is  nearly  uniform  across  the  chord.  This  angle  is  a 
modification  of  cxbe  resulting  from  shed  wake  and  can  be  con¬ 
sidered  as  the  overall  angle  of  attack  of  the  airfoil.  The 
equation  is 


“equ  *  0o  +  IF  -  Cc<3/4  -  PA)  G  -f]  0y 
+  tan-1  [ { H0  +  FHV 

+  [C  (3/4  -  PA) F  +  —  V]  (0  +  R  B  ) 
o  w  V 

+  C0(3/4  -  PA)  R  Bo)  /UT]  ,  (106) 

• 

where  H  represents  the  relative  velocity  normal  to  the  chord, 

0  is  the  mechanical  pitch  angle,  and  B  is  the  flap  angle.  The 
delta  angle  of  attack,  a^L'  results  directly  from  the  pitch¬ 
ing  velocity  which  gives  an  angle  of  attack  that  is  indepen¬ 
dent  of  shed  wake  and  varies  directly  as  the  chord, 


“DEL 


1/4  CD 


(0  +  RB) i 

U 
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This  term  is  treated  separately  since  its  impact  on  pitching 
moment  is  much  larger  than  the  other  lift  terms.  It  is 
arbitrarily  treated  as  a  perturbation  on  the  angle  of  attack 
and  uses  the  dynamic  lift  curve  slope  to  determine  its  con¬ 
tribution  to  CL. 

The  ability  of  this  theory  to  predict  CL  and  CM  has  been  eval¬ 
uated  by  comparing  calculated  values  with  test  data  obtained 
by  oscillating  an  airfoil  through  a  prescribed  pitch  motion. 


The  results  of  this  comparison  for  test  conditions  below 
stall,  at  stall,  and  beyond  stall  are  shown  in  Figure  43. 

METHOD  OF  CALCULATION 


Before  airloads  can  be  calculated,  the  blade  deflections  and 
downwash  field  must  be  known.  The  total  velocity  on  each 
aerodynamic  bay  of  the  rotor  blade  is  calculated  for  15-degree 
azimuth  positions  by  summing  the  downwash,  rotational  veloc¬ 
ity,  aircraft  forward  speed,  and  blade  motions.  The  pitch 
angle,  tangential  velocity,  and  axial  velocity  are  combined 
to  calculate  the  blade  element  (ag^)  angle  of  attack  (see 
equation  98) .  Next,  the  dynamic  angle  of  attack  for  lift  and 
pitching  moment  is  calculated  (apyw)  from  equation  9?  by  using 
the  stall  delay  parameter  y  and  the  pseudo-reduced  frequency. 
The  value  of  y  is  found  from  a  linear  interpolation  of  input 
Y  values  for  Mach  numbers  of  0,  0.2,  0.4,  0.6,  0.8,  and  1.0 
for  both  lift  and  pitching  moment. 

The  static  values  of  Cg  (at  agyN) /  CM  (at  aDYN )  >  and  CD  for 
equations  100  and  102  and  the  drag  calculations  are  found  from 
a  double  table.  The  tables  define  curves  of  the  aerodynamic 
coefficients  Cl ,  Cg,  Cm  versus  airfoil  angle  of  attack  from 
two-d j inensional  static  airfoil  tests. 

Compressibility  is  accounted  for  by  providing  separate  aero¬ 
dynamic  curves  for  a  range  of  Mach  numbers  from  0  to  1.  The 
tables  consist  of  punched-card  decks  which  must  be  input  with 
the  job  and  equations  which  are  contained  in  the  program.  A 
description  of  the  card  decks  and  equations  is  included  in  the 
user's  instruction  portion  of  this  document. 

Using  the  dynamic  angles  of  attack  (ctgyisj)  f°r  both  lift  and 
pitching  moments  and  local  Mach  number,  a  linear  double  inter¬ 
polation  with  Mach  number  and  angle  of  attack  is  used  to  find 
the  static  aerodynamic  coefficients.  Next,  the  equivalent 
angle  of  attack  (aeqU)  is  calculated  from  equation  105,  based 
upon  the  blade  deflections  and  the  Theodorsen  F  and  G  func¬ 
tions.  The  values  of  F  and  G  are  calculated  in  the  program; 
the  method  is  defined  in  detail  in  the  supplemental  document. 

Substituting  the  static  aerodynamic  coefficients  into  equa¬ 
tions  103  and  104  along  with  blade  motions  and  the  equivalent 
angle  of  attack  defines  the  dynamic  Cl  (Cl(dyn))  and  the 
dynamic  (cm(DYN) ) •  The  dynamic  Cq  is  set  equal  to  the 
value  of  Cg  obtained  from  the  static  tables. 

The  lift,  drag,  and  pitching  moment  are  then  calculated  in  the 
conventional  lifting-line  approach  using  the  dynamic  aerody¬ 
namic  coefficients;  therefore. 
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DYNAMIC  CL 


DYNAMIC  CM 


Figure  43.  Comparison  of  Test  and  Theoretical  Dynamic  Cl  and 
Cm  for  a  Condition  Below  Stall,  at  Stall,  and 
Beyond  Stall 
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L  =  Lift  =  l/2pV2C£.Cjj  (dyN) 

D  =  Drag  =  l/2pV2C£Cp (dyn) 

M  =  Pitching  Moment  =  l/2pV2C2£CM (DYn) 

The  wind  axis  airloads  (L,  D,  and  M)  are  then  converted  into 
blade  system  loads  for  every  aero  blade  station  at  15-degree 
increments.  The  dynamic  airloads  for  each  station  are  har¬ 
monically  analyzed  into  10  harmonics  which  then  become  part  of 
the  harmonic  forcing  loads  applied  in  the  response  analysis. 

NONUNIFORM  DOWNWASH 


Background 

The  initial  effort  in  the  representation  of  nonuniform  down- 
wash  in  the  rotor  loads  program  was  made  in  1960-61  with  an 
attempt  to  use  the  theory  of  Reference  26.  Later  the  work  in 
References  27  and  28  was  employed  in  digital  computer  pro¬ 
grams.  This  effort  was  minimal  and  these  programs  were  not 
readily  adaptable  for  integration  with  the  existing  rotor 
analysis.  Later,  effort  was  extended  to  use  the  tabulated 
induced  velocities  prepared  in  Reference  29;  but  about  this 
period  of  time,  Davenport  completed  his  computer  analysis 
which  was  reported  in  Reference  30.  This  program  can  compute 
nonuniform  induced  velocities  for  an  n-bladed  rotor  with  pre¬ 
scribed  loading,  either  to  yield  rotor  self-induced  or  rotor- 
rotor  interference  velocities  and  was,  in  fact,  the  first 
source  from  which  the  downwash  for  a  tandem  could  be  computed. 

This  analysis  was  suited  for  incorporation  as  a  subroutine  for 
the  rotor  loads  analysis,  but  it  had  two  basic  deficiencies 
in  that: 

1.  With  10  trailing  vortices  from  each  blade,  the  computer 
running  time  was  greatly  increased. 

2.  To  compute  the  necessary  airioads  for  input  to  this 
routine  required  that  it  be  made  a  part  of  the  existing 
iterative  loop  between  the  airload  and  the  blade  dynamic 
response  routines.  The  nonlinearities  which  already 
existed  in  the  program  from  nonlinear  aerodynamics  consid¬ 
erations  were  therefore  compounded  and  resulted  in  a  highly 
unstable  rotor  analysis  whi would  be  extremely  time- 
consuming  . 

To  reduce  the  analysis  computation  time,  the  downwash  analysis 
was  modified  to  reduce  the  trailed  vortices  on  each  blade  from 
10  to  2  and  to  employ  a  simplified  expression  for  the  vortex 
strength.  The  analysis  considered  one  vortex  trailed  at  the 
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blade  tip  and  the  other  from  the  root-end  blade  cutout.  The 
strength  of  the  trailed  vortices  was  calculated  by  considering 
the  blade  to  have  a  constant  bound  circulation  along  its  span 
from  the  blade  cutout  to  the  tip  of  sufficient  strength  to 
maintain  the  rotor  thrust.  Assuming  a  constant  bound  circula¬ 
tion  along  the  blade  span  eliminated  the  need  to  include  the 
downwash  calculations  in  the  airload-blade  response  iteration, 
since  the  trailed  vortex  can  be  calculated  from  the  assumed 
bound  circulation  without  knowing  the  actual  lift  distribution. 
Correlation  of  the  loads  program  with  test  data  (Reference  7) 
verified  the  acceptability  of  this  approach. 

h'ith  the  advent  of  faster  computers,  it  became  feasible  to 
use  the  actual  blade  lift  distribution  to  calculate  the 
trailed  vortex  strength.  Studies  showed  the  principal  factor 
perturbing  the  blade  lift  distribution  near  the  tip  and  cutout 
was  not  the  flexible  blade  deflections  but  the  downwash  it¬ 
self.  Figure  44  shows  that  the  downwash  alters  the  aerody¬ 
namic  angle  of  attack.  The  altered  angle  of  attack  changes 
the  lift  distribution  which  modifies  the  trailed  vortex 
strength.  Therefore,  as  a  first  approximation,  the  effect  of 
downwash  on  the  lift  will  be  included  and  the  blade  flexibil¬ 
ity  neglected  in  rhe  calculation  of  the  trailed  vortex 
strength . 

Therefore,  a  downwash  distribution  that  is  directly  a  function 
of  the  blade  tip  and  cutout  lift  distribution  is  obtained  by 
iterating  among  lift,  downwash,  and  vortex  strength  until 
there  is  compatibility  among  the  three.  Comparisons  of  the 
currently  used  compatibility  method  with  the  earlier  constant 
circulation  method  show  a  significant  improvement  in  load  pre¬ 
diction  ability  with  the  current  method  (see  References  15 
and  17)  . 

Method  of  Calculation 

The  nonuniform  downwash  is  calculated  by  considering  the  blade 
lift  distribution  at  15-degree  increments  around  the  rotor 
azimuth.  At  each  increment,  a  vortex  is  assumed  to  trail  from 
the  blade  cutout  and  blade  tip  circulation  discontinuities  of 
strength  determined  by  the  lift  of  the  adjacent  aerodynamic 
bay  (see  Figure  45) .  The  lift  and  hence  the  vortex  strength 
are  assumed  to  vary  linearly  from  azimuth  position  to  azimuth 
position.  Since  the  vortices  are  trailed  from  all  blades  on 
both  rotors  (for  a  tandem  helicopter) ,  a  complex  pattern  of 
vortex  spirals  is  formed  in  space.  By  summing  the  effects  of 
all  the  vortex  segments  on  a  given  blade  element,  the  nonuni¬ 
form  downwash  for  each  blade  element  is  determined. 

Next,  the  blades  are  rotated  15  degrees  and  the  previously 
trailed  vortices  drift  relative  to  the  hub  with  a  velocity 
composed  of  the  aircraft  forward  speed  and  the  uniform 
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NO  DOWNWASH 


WITH  DOWNWASH 


THE  DOWNWASH  FROM  THE  TRAILED  VORTICES  REDUCES  THE  ANGLE  OF 
ATTACK.  THE  REDUCTION  IN  ANGLE  OF  ATTACK  REDUCES  THE  LIFT  (LIFT 
DEFICIENCY)  AND  HENCE  REDUCES  THE  STRENGTH  OF  THE  TRAILED  VORTEX 
SINCE  THE  TRAILED  VORTEX  STRENGTH  IS  PROPORTIONAL  TO  THE  BLADE 
SECTION  LIFT. 


Figure  44.  Effect  of  Downwash 
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downwash.  New  trailed  vortices  are  calculated  for  this 
azimuth  position  and  the  downwash  of  all  the  vortex  segments 
is  again  summed  on  each  blade  element.  Then  the  blades  are 
incremented  another  15  degrees  to  the  next  azimuth  position. 

When  the  whole  rotor  azimuth  has  been  covered  by  this  stepwise 
procedure,  a  nonuniform  downwash  field  has  been  defined  by  a 
downwash  value  at  each  blade  airfoil  section  all  around  the 
azimuth  at  15-degree  increments.  A  description  of  the  down- 
wash  calculations  follows. 


Theory 

Reference  30  presents  a  detailed  description  of  the  analytical 
method  which  is  based  on  the  assumption  of  a  rigid  wake.  The 
essential  steps  involved  in  the  computation  of  rotor-induced 
velocities  are  summarized  below. 


Considering  the  origin  as  the  center  of  rotation  of  the  rotor 
trailing  the  vortex,  the  position  of  the  vortex  filament  is 
established  from  the  vector  diagram  in  Figure  46. 

It  is  shown  that: 


Qx  =  r  cos  8(6)  cos  0  +  yx(tjj~e) 
Qy  =  -r  cos  8(0)  sin  0  +  Vy  ( ^— 6 ) 
Qz  =  -r  sin  8(0)  +  yz(^-0) 


where  r 
R 

Q 

fi 


ip ,  9 


^x'  yy' 


radial  distance  along  blade 
rotor  radius 

vortex  filament  position  vector  with 
components  Qx,  Qy,  and  Qz 

rotor  angular  velocity 

rotor  blade  flapping  angle 

rotor  blade  azimuth  angles 

ratio  of  components  of  vortex  system 
drift  velocity  to  EUQ 


Computing  the  Q  vector  components  for  a  series  of  points  spaced 
at  equal  increments  in  time  (i.e.,  equal  increments  in  <p)  ,  the 
geometry  of  the  complete  vortex  filament  is  defined.  Between 
successive  computed  points  the  filament  is  represented  by 
straight  line  segments.  A  linear  variation  of  the  circulation 
strength  along  the  filament  is  assumed. 
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Figure  46.  Vortex  Vector  Diagram 


The  segment  between  successive  points  A  and  B  along  the  fila¬ 
ment  is  given  by 


Then,  let  P  be  the  vector  from  the  origin  to  the  point  at 
which  it  is  desired  to  compute  the  velocity.  The  vectors  be¬ 
tween  the  points  A  and  B  on  the  vortex  filament  and  the  point 
at  which  the  downwash  is  to  be  obtained  are 

SA  =  Qa  “  p 

SB  =  Qb  ”  p 

For  defining  circulation  along  the  filament,  the  positive  di¬ 
rection  is  a  vector  pointing  from  A  toward  B  (i.e,  downstream 
along  the  filament)  using  the  right-hand  rule  for  the  sense  o 
rotating.  Positive  circulation  would  therefore  induce  a  ve¬ 
locity  coming  out  of  the  paper  as  shown  below. 


u»-~ 


A  unit  vector  in  the  direction  of  the  velocity  induced  by  the 
segment  is  given  by 

SA  x  SB 


Using  the  Biot-Savart  law,  the  magnitude  of  the  induced  veloc¬ 
ity  can  be  written 


V 


B 

/ 


A 


r  p 

4tt  R3 


ds 


Integrating  the  above  expression  for  a  linear  variation  of  the 
circulation , 


V  = 


4  tt  P  d 
rB 


-  £d  +  i2  +  P2 

✓  0  2  +  p  2  -  . 

v  u  +  d) 2  +  P2 


4tt  P  d 


^ U  +  d) 2  +  P2  + 


£d  +  £2  4-  P2 
)/£ 2  +  P2 


Velocity  components  induced  at  P  by  the  filament  A  B  are  those 
of  the  vector. 


V  =  V  U 

The  circulation  value  at  the  corners  A  and  B  of  each  filament 
is  defined  from  the  spanwise  lift  distribution  at  the  ip  values 
corresponding  to  the  corners.  Bound  circulation  at  any  radial 
station  on  the  blade  for  a  particular  y  can  be  found  by  the 
Kutta-Joukowski  law: 

Lift  =  (density) (speed) (circulation)  =  pVI 
From  the  above,  the  nondimensional  bound  circulation  is 

CL  r 

r  =  C/R  —  (B  +  y  cos  &s  sin  ip)  , 

where  the  term  in  parentheses  is  an  approximation  to  the  blade 
element  velocity  and 
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CL  =  blade  element  lift  coefficient 

C  =  blade  chord  length 

y  =  ratio  of  helicopter  flight  speed  to  fiR 

as  =  shaft  angle  of  attack 

r  =  radial  distance  along  blade 

R  =  rotor  radius 

ip  =  rotor  blade  azimuth  angle 

Using  the  above  expression,  the  strength  of  each  trailed 
vortex  is  established  from  the  lift  distribution  along  the 
blade  at  each  \p  position.  The  velocity  at  point  P  is  then 
obtained  by  successive  additions  of  the  contribution  from 
each  segment  of  each  vortex,  going  back  along  the  filament 
spirals  far  enough  to  include  all  substantial  contributions. 


Part  IV.  Prop/Rotor  Loads 
Analysis  Correlation 

R.  B.  Taylor 
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INTRODUCTION 


Extensive  loads  correlation  was  performed  with  the  C--70  analy¬ 
sis  to  determine  its  prediction  capability  for  a  wide  range  of 
rotors  and  test  conditions.  For  correlation  of  the  analysis 
with  highly  twisted  prop/rotors,  model  test  data  from  3  prop/ 
rotor  wind  tunnel  tests  was  used.  The  test  data  range  in¬ 
cluded  hover,  transition,  and  cruise  for  the  tilt-wing  and 
tilt-rotor  aircraft,  in  addition  to  windmilling  cruise  for  the 
stowed-rotor  aircraft.  Blade  flap  and  chord  bending,  in  addi¬ 
tion  to  blade  torsion,  were  of  primary  concern  fcr  this  part 
of  the  correlation.  In  addition,  steady  rotor  forces  and 
momercs  were  also  correlated  with  the  prop/rotor  test  data. 
That  portion  of  the  analysis  which  predicts  loads  for  low- 
twisted  rotors  for  which  flapwise  and  lagwise  motions  are 
largely  uncoupled  (C-60)  was  used  for  helicopter  loads  corre¬ 
lation.  For  correlation  of  the  analysis  with  low- twisted 
articulated  rotors,  flight  test  data  for  the  CH-47  and  Model 
347  was  used  in  addition  to  CH-34  full-scale  wind  tunnel  test 
data.  Correlation  of  pitch  link  loads,  blade  flap,  and  chord 
bending  was  of  primary  concern  for  this  part  of  the  correla¬ 
tion.  Also,  airloads  were  correlated  with  CH-34  airloads  test 
data. 
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PROP/ROTOR  LOADS  CORRELATION 


INTRODUCTION 


Blade  loads  and  rotor  loads  test  data  from  three  hingeless 
prop/rotor  wind  tunnel  models  was  used  for  correlation  with 
the  C-70  analysis  to  determine  its  prediction  capability  for 
highly  twisted  prop/rotors.  The  wind  tunnel  models  selected 
provided  a  wide  range  of  flight  conditions  and  also  a  large 
variation  of  blade  dynamic  properties.  Two  of  the  models  were 
stiff  in  plane  (first  chordwise  bending  frequency  greater  than 
operating  rotor  speed) ,  and  the  other  model  was  soft  in  plane 
(first  chordwise  bending  frequency  less  than  operating  rotor 
speed) .  Due  to  the  high  twist,  all  three  model  rotors  had 
highly  coupled  flapwise  and  lagwise  motions  and  therefore  pro¬ 
vided  a  good  test  of  the  flap-chord  coupling  capabilities  of 
the  C-70  analysis. 

Table  I  is  a  brief  summary  of  the  prop/rotor  models  used  for 
correlation.  The  160  performance  model  is  a  1/10-scale,  full- 
span  powered  wind  tunnel  model  of  the  Model  160  tilt-rotor 
aircraft,  a  current  Boeing-Vertol  tilt-rotor  design.  The  wind 
tunnel  model  has  stif f-in-plane  blades;  and  since  the  model  is 
intended  to  obtain  primarily  performance  data,  the  blades  are 
geometrically  scaled  only.  The  1/3-scale  LIT  model  is  a  pow¬ 
ered,  4-bladed  rotor  mounted  on  a  wing.  The  blades  are 
Froude-scaled  and  are  stiff  in  plane.  The  1/9-scale  conver¬ 
sion  model  is  a  half-span,  unpowered  windmilling  model  of  the 
Model  213  stowed-rotor  aircraft  which  is  a  current  Boeing- 
Vertol  stowed-rotor  design.  The  4-bJaded  rotor  and  the  wing 
are  Froude-scaled  and  the  blades  are  soft  in  plane. 

The  test  data  obtained  from  these  models  covered  the  3  main 
flight  regimes  of  the  tilt-rotor  and  tilt-wing  aircraft:  hover, 
transition,  and  cruise.  In  addition,  as  shown  in  Table  II, 
data  was  obtained  for  correlation  in  windmilling  cruise,  an 
important  flight  regime  of  the  stowed-rotor  aircraft. 

Table  III  is  a  summary  of  the  prop/rotor  model  test  data  used 
for  correlation.  As  shown,  blade  loads  and  steady  rotor  hub 
moments  for  all  three  models  were  correlated  with  the  analysis. 

CORRELATION  RESULTS 

Hover 


Test  data  from  two  prop/rotor  models  was  used  for  correlation 
with  the  analysis  in  hover:  the  1/3-scale  LIT  rotor  and  the 
Model  160  performance  model.  Both  models  were  tested  out  of 
ground  effect  to  reduce  interference  effects;  cyclic  applica¬ 
tion  with  the  rotor  disc  plane  in  a  horizontal  position  pro¬ 
duced  the  blade  loads  and  hub  moments  for  the  hover 
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TABLE  II 

SUMMARY  OF  TEST  CONDITIONS  OF  PROP/ROTOR 
MODELS  USED  FOR  CORRELATION 


Model 

Hover 

Transition 

Powered 

Cruise 

Windmilling 

Cruise 

160  Performance 

X 

- 

- 

- 

1/3- Scale  LIT 

X 

X 

X 

- 

1/9-Scale  Conversion 

- 

- 

- 

X 

TABLE  III 

SUMMARY  OF  PROP/ROTOR  MODEL  TEST  DATA  USED  FOR  CORRELATION 


Model 

Blade 

Flap 

Bending 

Blade 

Chord 

Bending 

blade 

Torsion 

Steady  Rotor 
In-Plane 
Moments 

160  Performance 

X 

X 

V 

4\ 

X 

1/3-Scale  LIT 

X 

X 

X 

X 

1/9-Scale  Convers: 

Lon  X 

X 

X 

X 

c  • 

i 

!  i 


fs 

i 


correlation.  The  correlation  is  presented  as  a  comparison  of 
predicted  and  measured  alternating  blade  moments,  first-  and 
second-harmonic  blade  moments,  and  blade  moment  waveforms. 

Wind  tunnel  testing  experience  with  prop/rotors  has  shown  that 
alternating  bending  moments  are  nearly  all  1-  and  2-per-rev, 
with  the  exception  of  blade  torsion  near  stall  flutter  and 
chord  bending  in  conditions  where  pitch-lag  coupling  is 
evident.  For  a  condition  such  as  this,  the  waveforms  clearly 
show  the  higher  harmonic  content.  Since  the  prediction  of 
steady  rotor  moments  due  to  hover  cyclic  is  closely  related  to 
and  as  important  as  the  prediction  of  blade  bending  moments, 
correlation  is  also  performed  with  steady  rotor  in-plane 
moments . 

Figures  47  through  56  are  the  hover  cyclic  correlation  results 
for  the  1/3-scale  LIT  rotor.  As  mentioned  in  the  introduction , 
this  4-bladed  rotor  was  mounted  on  a  wing  and  the  blades  are 
stiff  in  plane  with  40  degrees  of  twist.  Correlation  with 
alternating  blade  root  flap  bending  (one  half  peak-to-peak)  in 
Figure  47  is  very  good,  with  the  analysis  slightly  over¬ 
predicting  flap  bending  sensitivity  to  cyclic.  Correlation  of 
the  first  two  flap  bending  harmonics,  which  are  the  prime  har¬ 
monics  of  the  alternating  moment,  is  shown  in  Figure  48.  The 
analysis  correlates  very  well  with  the  first-harmonic  flap 
bending  sensitivity  to  cyclic  but  overpredicts  second-harmonic 
flap  bending.  The  flap  bending  waveform  in  Figure  49  shows 
that  the  predicted  first  harmonic  agrees  well  in  phase  and 
magnitude  with  the  measured  data,  but  inclusion  of  the  predicted 
second  harmonic  shifts  the  total  waveform  peak  to  150  degrees 
azimuth.  The  alternating  blade  root  chord  bending  correlation 
is  shown  in  Figure  50,  and  the  prediction  again  agrees  very 
well  with  the  test  data.  Correlation  of  the  first  two  chord 
bending  harmonics  in  Figure  51  shows  the  analysis  overpredicts 
first-harmonic  chord  bending  by  13  percent  and  agrees  very 
well  with  the  measured  low  value  of  2-per-rev  chord  bending. 
Figure  52  shows  the  predicted  chord  bending  waveform  has  a 
large  3-per-rev  contribution;  this  can  be  traced  to  the  over¬ 
prediction  of  2-per-rev  flap  bending  combined  with  the  large 
1-per-rev  flap  bending,  since  1-  and  2-per-rev  flapping  pro¬ 
duce  3-per-rev  and  4-per-rev  lag  motion.  Note  that  the  pre¬ 
dicted  f irst-harmonj c  waveform  agrees  well  in  magnitude  and 
fairly  well  in  phase  with  the  measured  waveform.  Correlation 
results  for  alternating  blade  root  torsion  are  presented  in 
Figures  53  and  54.  r,1he  analysis  overpredicts  alternating 
torsion  by  14  percent  in  Figure  53  and  the  predicted  waveform 
in  Figure  54  shows  a  large  3-per-rev  contribution.  This  indi¬ 
cates  a  strong  coupling  in  the  analysis  between  blade  pitching 
and  lag  motion,  probably  due  to  the  aft  shear  center  location 
(35  percent  chord)  for  this  blade.  The  Model  160  performance 
model  blade  loads  test  data  showed  this  type  of  coupling  in 
hover;  these  results  will  be  discussed  later  in  this  section. 
Correlation  with  the  steady  rotor  pitching  and  yawing  moments 
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is  shown  in  Figures  55  and  56.  A  steady  pitching  moment 
pitches  the  aircraft  nose  up  and  the  yawing  moment  rolls  the 
aircraft  in  hover.  The  results  show  the  analysis  agrees  very 
well  with  the  test  data  for  both  moments.  These  results  are 
important  sirce  they  demonstrate  the  capability  of  the  analy¬ 
sis  to  predict  the  ratio  of  rotor  hub  moment- to-b lade  bending 
moment  amplitude  in  hover.  This  ratio  is  critical  to  prop/ 
rotor  design  since  the  bending  fatigue  design  loads  are  those 
produced  by  hover  cyclic  to  produce  required  steady  rotor 
forces  and  moments  for  aircraft  control  in  hc'rer. 

Figures  57  through  66  are  the  hover  cyclic  correlation  results 
for  the  Model  160  performance  model.  As  mentioned  in  the 
introduction,  this  model  has  a  3-bladed  rotor  and  the  blades 
are  stiff  in  plane.  The  cyclic  pitch  was  phased  to  yield  pure 
longitudinal  blade  flapping  so  that  the  total  steady  hub 
moment  would  pitch  the  aircraft  with  little  aircraft  roll. 

In  terms  of  prop/rotor  nomenclature,  the  desired  effect  was 
to  produce  a  rotor  pitching  moment  with  little  or  no  yawing 
moment.  The  blade  root  flap  bending  correlation  in  Figures 
57,  58,  and  59  shows  the  analysis  agrees  well  with  the  test 
data,  both  in  alternating  value  (one-half  peak-to-peak)  and 
the  first  harmonic.  The  measured  second  harmonic  is  very  low 
compared  to  the  first  harmonic,  and  the  prediction  for  the 
second  harmonic  agrees  well  with  the  measured  data  above 
2,200  rpm.  Correlation  with  the  measured  flap  bending  wave¬ 
form  in  Figure  59  is  very  good,  both  in  amplitude  and  phase. 
The  good  phase  correlation  demonstrates  that  the  blade  dynamic 
response,  in  terms  of  phase  lag,  to  a  1-per-rev  cyclic  input 
is  correctly  accounted  for  in  the  analysis.  This  is  an  impor¬ 
tant  result  since  the  first  flapping  mode  deflection  at  the 
blade  root  is  nearly  perpendicular  to  the  blade  major  prin¬ 
cipal  axis,  which  is  inclined  nearly  40  degrees  with  the  disc 
plane  at  the  blade  root  due  to  blade  twist  and  collective 
pitch.  Thus,  in-  and  out-of-plane  motion  are  coupled  for 
blade  root  flap  bending  and  the  analysis  has  shown  the  capa¬ 
bility  to  predict  this  coupling  effect  on  flap  bending  moment. 

The  results  of  the  hover  cyclic  chord  bending  correlation  for 
the  Model  160  performance  model  are  presented  in  Figures  60, 
61,  and  62.  Predicted  alternating  chord  bending  agrees  well 
with  the  test  data.  Both  prediction  and  test  data  are  rela¬ 
tively  insensitive  to  changes  i.n  rotor  speed.  Correlation  for 
the  first  and  second  harmonics  in  Figure  61  shows  overpredic¬ 
tion  of  the  first  harmonic  and  underprediction  of  the  second 
harmonic.  Poor  correlation  of  the  first  harmonic  can  be  due 
to  the  blade  lag  root  boundary  in  the  analysis.  The  blade  was 
modeled  as  a  cantilevered  beam  such  that  the  lag  deflection 
slope  ac  the  root  boundary  was  zero.  A  detailed  analysis  of 
the  blade  root  retention,  including  bearings,  indicated  that 
significant  lag  motion  in  the  bearing  would  be  present;  and 
therefore,  the  blade  lag  root  retention  would  be  more 


accurately  defined  as  pinned  with  a  stiff  angular  spring  in¬ 
stead  of  cantilevered.  The  effect  of  changing  the  root 
boundary  condition  from  cantilevered  to  pinned  with  an  angular 
spring  would  be  small  on  first  chord  frequency  and  also  small 
on  lag  dynamic  response,  but  the  effect  on  chord  bending 
moment  radial  distribution  could  be  significant.  For  a  rota¬ 
ting  cantilevered  beam,  the  first  modal  bending  moment  distri¬ 
bution  in  a  vacuum  steadily  increases  from  zero  at  the  end  of 
the  beam  to  a  maximum  value  at  the  cantilevered  boundary.  The 
effect  of  pinning  the  beam  and  applying  an  angular  spring  is 
to  decrease  the  bending  moment  value  at  the  root  boundary  from 
that  for  the  cantilevered  condition,  even  though  tip  deflec¬ 
tions  fc1;  both  conditions  are  equal.  Therefore,  modeling  the 
blade  lag  root  retention  in  the  analysis  as  a  pin  with  angular 
spring  should  decrease  the  predicted  first-harmonic  chord 
bending f  with  little  effect  on  second-harmonic  chord  bending. 

The  waveform  correlation  for  chord  bending  in  Figure  62  shows 
that  both  prediction  and  measured  data  have  significant  higher- 
harmonic  contributions.  A  large  6-per-rev  contribution  is 
evident,  particularly  between  90  and  210  degrees  azimuth. 

Blade  lag  motion  is  coupling  with  blade  pitching  motion,  as 
shown  by  a  comparison  of  the  chord  bending  and  torsional  wave¬ 
forms  in  Figures  62  and  63,  respectively.  The  measured  data 
shows  that  chord  bending  and  torsion  are  oscillating  in  phase 
at  6  per  rev.  The  blade  first  torsional  frequency  ratio  is 
near  6  per  rev  at  this  rotor  speed  and  the  rotor  is  probably 
on  the  verge  of  stall  flutter  at  this  point.  Stall  flutter 
inception  is  characterized  by  a  large  growth  in  alternating 
blade  torsion  oscillating  at  the  torsional  natural  frequency. 
Coupling  of  the  torsional  6-per-rev  motion  with  the  lag  mode 
can  be  traced  to  the  aft  shear  center  location  for  this  blade 
(at  38  percent  chord) .  Pitching  moment  at  6  per  rev  induces 
6-per-rev  airloads  which  in  turn  force  the  chord  mode  at  6  per 
rev  due  to  induced  drag.  When  the  shear  center  is  offset  from 
the  aerodynamic  center,  the  oscillating  airloads  increase  the 
blade  angle  of  attack  at  the  same  frequency,  thus  amplifying 
6-per-rev  airloads  and  therefore  6-per-rev  chord  bending.  The 
analysis  shows  the  same  tendency  to  oscillate  at  6  per  rev, 
both  in  chord  bending  and  torsion.  In  Figure  62,  note  that 
the  phasing  of  the  predicted  6-per-rev  component  is  in  phase 
with  the  measured  data,  but  the  analysis  has  gone  into  pre¬ 
mature  stall  flutter  as  shown  by  the  large  6-per-rev  spiking 
in  the  second  quadrant.  Note  that  the  sixth  harmonic,  in  the 
predicted  chord  bending  waveform  (Figure  63)  is  in  phase  with 
the  predicted  6  per  rev  in  the  torsional  waveform.  This  re¬ 
sult  is  important  since  it  shows  the  capability  of  the  analy¬ 
sis  to  predict  pitch-lag  coupling.  A  comparison  of  predicted 
and  measured  alternating  blade  torsion  in  Figure  64  shows  good 
correlation  with  the  analysis  slightly  overpredicting .  The 
overprediction  is  due  to  premature  stall  flutter  inception  in 
the  analysis. 
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Correlation  results  for  the  steady  rotor  in-plane  moments  d 
to  hover  cyclic  for  the  Model  160  performance  model  are  she 
in  Figures  65  and  66.  Correlation  for  both  pitching  and  ya 
ing  moments  is  good,  with  the  analysis  overpredicting  the 
pitching  moment  by  10  percent  at  2,300  rpm. 
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Figure  50.  1/3-Scale  LIT  Blade,  Predicted  and  Measured 
Alternating  Chord  Bending  in  Hover  Due  to 
Cyclic  Pitch 
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Figure  51.  1/3-Scale  LIT  Blade,  Predicted  and  Measured 
First-  and  Second-Harmonic  Chord  Bending  in 
Hover  Due  to  Cyclic  Pitch 
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Figure  55.  1/3-Scale  LIT  Blade,  Predicted  and  Measured 
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Figure  57.  Model  160  Performance  Model,  Predicted  and 
Measured  Alternating  Flap  Bending  Moment  in 
Hover  for  3  Degrees  Cyclic  Pitch 
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ALTERNATING  TORSION  AT  .325R  (IN-LB) 
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Figure  64.  Model  160  Performance  Model,  Predicted  and 
Measured  Alternating  Torsion  in  Hover  for 
3  Degrees  Cyclic  Pitch 
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STEADY  ROTOR  PITCHING  MOMENT 
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Powered  Cruise 


Test  data  from  the  ] /3-scale  LIT  rotor  was  used  for  correla¬ 
tion  in  powered  cruise.  The  rotor  shaft  was  horizontal  in  a 
propeller  configuration  with  rotor  thrust  producing  the  total 
propulsive  force.  Rotor  angle  of  attack  was  varied  with  all 
other  parameters  remaining  constant.  The  results  of  the  pow¬ 
ered  cruise  correlation  are  shown  in  Figures  67  through  76. 
Alternating  blade  root  flap  bending,  chord  bending,  and  tor¬ 
sion  correlation  (Figures  67,  68,  and  69)  are  very  good,  with 
the  analysis  predicting  within  at  least  5  percent  of  the  mea¬ 
sured  bending  moment  data  and  within  15  percent  of  the  measured 
blade  root  torsion.  Correlation  with  the  first  and  second 
bending  moment  harmonics  (Figures  70  and  71)  is  also  very  good 
for  first-harmonic  flap  and  chord  bending.  The  analysis  under¬ 
predicts  second-harmonic  flap  and  chord  bending.  Interference 
due  to  wing  circulation  probably  induced  measured  higher- 
harmonic  bending  moments,  and  this  effect  was  not  included  in 
the  analysis  predictions.  Correlation  results  for  bending 
moment  and  blade  torsion  waveforms  (Figures  72,  73,  and  74) 
are  very  good.  The  analysis  predicts  the  correct  amplitude 
for  flap,  chord,  and  torsion,  but  the  phase  of  the  predicted 
first-harmonic  contribution  is  shifted  by  about  30  degrees. 

The  analysis  predicts  peak  moments  to  occur  at  90  degrees, 
whereas  test  data  peaks  at  about  120  degrees.  The  impact  of 
the  phase  shift  on  steady  hub  moments  is  large  as  shown  in  the 
hub  moment  correlation  in  Figures  75  and  76.  The  analysis 
underpredicts  steady  rotor  pitching  moment  by  a  large  margin 
but  correlates  very  well  with  rotor  yawing  moment.  A  steady 
pitching  moment  pitches  the  aircraft  nose  up  and  a  yawing 
moment  yaws  the  aircraft  in  cruise. 

The  poor  pitching  moment  correlation  can  be  traced  to  the 
source  of  its  forcing.  The  soerces  of  rotor  pitching  moment 
and  yawing  moment  are  independent  of  each  other  due  to  the 
dynamic  response  characteristics  of  the  blades.  Since  the 
blade  first-flapping  frequency  ratio  is  so  high  (w/ft  =  2.7), 
the  flapping  phase  lag  is  only  10  degrees.  Therefore  the 
blade  response  to  1-per-rev  forcing  is  nearly  in  phase  with 
the  force.  This  implies  that  rotor  inflow  induced  by  rotor 
angle  of  attack  will  produce  only  lateral  1-per-rev  flapping 
and  therefore  a  pure  rotor  yawing  moment.  Likewise,  cyclic 
phased  170  degrees  from  downwind  will  produce  only  longitu¬ 
dinal  flapping  and  therefore  a  pure  rotor  pitching  moment; 
therefore  the  source  of  the  poor  pitching  moment  correlation 
must  be  due  to  cyclic.  The  very  good  pitching  moment  corre¬ 
lation  in  hover  due  to  cyclic  for  this  rotor  demonstrated  the 
capability  of  the  analysis  to  predict  rotor  pitching  moment  due 
to  cyclic,  and  therefore  the  measured  cyclic  which  was  input  in 
the  analysis  is  suspected  to  be  incorrect.  This  is  plausible 
since  a  small  amount  of  cyclic  has  a  large  impact  on  rotor  pitch¬ 
ing  moment.  As  shown  in  Figure  55,  one  degree  of  cyclic  produces 
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150  foot-pounds  of  steady  pitching  moment  and  moment  sensi¬ 
tivity  to  cyclic  is  even  greater  in  cruise. 

Summarizing  the  powered  cruise  correlation,  the  results  are 
very  encouraging  since  the  analysis  demonstrated  a  very  good 
capability  to  predict  rotor  loads  for  a  condition  in  which  the 
flapwise  and  chordwise  modes  are  highly  coupled.  Due  to  the 
high  collective  pitch  required  for  cruise  flight  (in  this 
case,  38  degrees) ,  the  modes  become  more  highly  coupled  and 
aerodynamic  damping  plays  a  greater  roll  in  the  response  of 
both  flapping  and  lagging.  Use  of  an  uncoupled  analysis  for 
rotor  loads  predictions  would  provide  at  best  crude  approxi¬ 
mations  in  cruise,  particularly  for  chord  bending,  and  the 
cruise  correlation  results  demonstrate  the  quality  of  predic¬ 
tions  that  can  be  obtained  when  the  modes  are  fully  coupled. 
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ALTERNATING  FLAP  BENDING  AT  .22R  (IN-LB) 


TERNATING  CHORD  BENDING  AT  .22R  (IN-LB) 
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Figure  68. 


1/3-Scale  LIT  Blade,  Predicted  and  Measured 
Alternating  Chord  Bending  Due  to  Rotor  Angle 
of  Attack  in  Cruise 
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Figure  70.  1/3-Scale  LIT  Blade,  Predicted  and  Measured 
First-  and  Second-Harmonic  Flap  Bending  Due 
to  Rotor  Angle  of  Attack  in  Cruise 
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Figure  73.  1/3-Scale  LIT  Blade,  Predicted  and  Measured 
Chord  Bending  Waveform  for  5.1  Degrees  Angle 
of  Attack  in  Cruise 
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Figure  75.  1/3-Scale  LIT  Blade,  Predicted  and  Measured 
Steady  Rotor  Pitching  Moment  Due  to  Rotor 
Angle  of  Attack  in  Cruise 
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STEADY  ROTOR  YAWING  MOMENT  (FT-LB) 


Windmilling  Cruise 

Figures  77  through  83  are  the  cruise  correlation  results  for 
the  1/9-scale  conversion  model.  As  mentioned  in  the  introduc¬ 
tion,  this  model  is  a  1/9-scale,  half-span  model  of  a  Boeing- 
Vertol  stowed-rotor  aircraft  design.  The  4-bladed  rotor  is 
Froude-scaled,  as  is  the  wing,  and  the  rotor  is  unpowered. 

The  blades  are  soft  in  plane  with  a  first  chordwise  frequency 
ratio  of  0.75  u)/S2  at  cruise  rpm.  Windmilling  is  produced  by 
changing  blade  collective  pitch.  Rotor  loads  data  were  taken 
at  a  cruise  condition  which  simulated  the  end  of  transition 
flight.  On  a  full-scale  design,  the  aircraft  would  now  be 
jet-powered  with  the  rotors  windmilling  and  ready  for  stowing. 
The  model  rotor  angle  of  attack  was  varied,  with  tunnel  speed 
and  rotor  speed  being  held  constant. 

The  blade  root  flap  bending  correlation  in  Figure  77  shows 
that  the  analysis  correctly  predicts  the  slope  of  first- 
harmonic  flap  bending  with  rotor  angle  of  attack,  but  the  pre¬ 
dicted  magnitude  is  high  on  a  poini_-to-point  correlation  basis 
No  explanation  can  be  given  for  the  overprediction,  but  the 
good  correlation  of  flap  bending  moment  per  degree  of  angle  of 
attack  is  encouraging.  The  analysis  underpredicts  second- 
harmonic  flap  bending,  showing  very  little  2-per-rev  response. 
Second-harmonic  flap  bending  was  also  underpredicted  for  the 
1/3-scale  LIT  rotor  in  cruise  and,  as  for  that  rotor,  the 
large  measured  value  of  2-per-rev  flap  bending  is  suspected  to 
be  due  to  wing  interference  effects.  Correlation  with  the 
flap  bending  waveform  in  Figure  78  shows  the  analysis  overpre¬ 
dicts  in  amplitude  but  agrees  well  in  phase  with  the  measured 
first-harmonic  content.  Chord  bending  correlation  results  are 
presented  in  Figure  79.  The  analysis  correctly  predicts  the 
trend  of  first-harmonic  chord  bending  with  rotor  angle  of 
attack  tut  overpredicts  the  magnitude.  The  reason  for  this  is 
not  known.  Like  the  flap  bending  correlation,  the  analysis 
also  underpredicts  second-harmonic  chord  bending;  this  is  at¬ 
tributed  to  wing  interference  effects  on  the  measured  chord 
bending.  Correlation  with  the  measured  chord  bending  waveform 
in  Figure  80  is  good,  both  in  amplitude  and  phase. 

Correlation  with  alternating  blade  root  torsion  is  shown  in 
Figure  81.  The  analysis  overpredicts  by  a  large  margin, 
probably  due  to  an  incorrect  torsional  root  spring  descrip¬ 
tion.  Blade  pitch  on  this  rotor  was  controlled  by  a  worm  gear 
in  the  center  of  the  hub  to  which  all  4  blades  were  geared. 
This  type  of  pitch  control  was  used  instead  of  the  standard 
swashplate  configuration  to  attain  a  rapid  and  equal  rate  of 
change  of  collective  on  the  rotor.  The  torsional  root  spring 
due  to  the  gearing  could  not  be  accurately  calculated  and 
therefore  an  estimation  was  made. 
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Correlation  with  steady  rotor  in-piane  moments  due  to  rotor 
ang]e  of  attack  is  shown  in  Figures  82  and  83.  Both  pitching 
moment  and  yawing  moment  correlate  well  with  the  measured 
data. 
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FIRST  HARMONIC  FLAP  BENDING  AT  .  125R(IN-LB)  SECOND  HARMONIC  FLAP  BENDING 


FIRST  HARMONIC  CHORD  BENDING  AT  SECOND  HARMONIC  CHORD 
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Figure  79.  Model  213  I/9-Scale  Conversion  Model,  Predicted 
and  Measured  First-  and  Second-Harmonic  Chord 
Bending  at  0.125R  Due  to  Rotor  Angle  of  Attack 
in  Cruise 
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Figure  80.  Model  213  1/9-Scale  Conversion  Model,  Predicted 

and  Measured  Chord  Bending  Waveform  at  0.125R  for 
2.26  Degrees  Rotor  Angle  of  Attack  in  Cruise 


Figure  81. 


Model  213  1/9-Scale  Conversion  Model,  Predicted 
and  Measured  Alternating  Torsion  at  0.150R  Due 
to  Rotor  Angle  of  Attack  in  Cruise 
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STEADY  ROTOR  PITCHING  MOMENT  (FT-LB) 
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Figure  82. 


Model  213  1/9-Scale  Conversion  Mr-"1 
and  Measured  Steady  Rotor  Pitchi: 
Rotor  Angle  of  Attack  in  Cruise 
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Figure  83.  Model  213  1/9-Scale  Conversion  Model,  Predicted 
and  Measured  Steady  Rotor  Yawing  Moment  Due  to 
Rotor  Angle  of  Attack  in  Cruise 
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Transition 


Test  data  obtained  from  the  1/3-scale  LIT  rotor  was  used  for 
transition  correlation.  With  the  rotor  shaft  angle  held  con¬ 
stant  at  30  degrees  from  the  horizontal,  cyclic  pitch  was 
varied.  Results  of  the  blade  root  flap  bending  correlation 
are  shown  in  Figures  84  through  86.  The  prediction  correlates 
well  with  alternating  flap  bending  up  to  2  degrees  cyclic, 
pitch,  but  the  analysis  does  not  reflect  the  nonlinear  increase 
in  flap  bending  above  2  degrees  cyclic.  For  negative  cyclic 
pitch,  the  analysis  predicts  minimum  flap  bending  at  minus 
1  degree  cyclic  compared  with  a  measured  minimum  at  minus  2 
degrees.  However,  the  analysis  does  correlate  well  with  the 
sensitivity  of  flap  bending  to  negative  cyclic.  Correlation 
with  the  first  two  flap  bending  harmonics  in  Figure  85  shows 
the  same  trend  for  the  first  harmonic  as  the  alternating  value 
in  Figure  84,  ai.i  underprediction  of  second  harmonic  flap 
bending  for  positive  cyclic.  The  flap  bending  waveform  cor¬ 
relation  in  Figure  86  is  good  in  amplitude  but  the  predicted 
maximum  value  is  at  120  degrees  azimuth  compared  with  a  mea¬ 
sured  value  of  165  degrees;  the  analysis  is  underpredicting 
longitudinal  flapping  for  this  condition.  In  light  of  the 
small  phase  laq  response  (about  10  degrees)  of  the  flapping 
mode  for  this  rotor  speed,  the  measured  longitudinal  flap  bend¬ 
ing  is  abnormally  large.  The  forcing  producing  the  flap  bend¬ 
ing  moment  in  Figure  86  is  all  due  to  a  1-per-rev  variation  of 
rotor  inflow,  since  cyclic  pitch  is  zero  for  this  test  point. 
Since  aerodynamic  forcing  due  to  rotor  angle  of  attack  is  a 
maximum  on  the  advancing  blade  and  a  minimum  on  rhe  retreating 
blade  and  the  flapping  phase  lag  is  10  degrees,  the  blade 
flapping  response  and  therefore  flap  bending  moment  would  be 
expected  to  be  predominately  in  the  lateral  direction.  Down- 
wash  nonuniformity  due  to  a  strong  tip  vortex  interaction  could 
cause  a  shift  in  airloads  but,  since  uniform  downwash  was  used 
in  the  analysis  for  this  part  of  the  correlation,  any  such 
effects  would  not  be  reflected  in  the  predictions  presented 
herein. 

The  transition  chord  bending  correlation  results  are  shown  in 
Figures  87  through  89.  The  analysis  correlates  well  for  posi¬ 
tive  cyclic  pitch  but  predicts  the  same  sensitivity  of  chord 
bending  for  negative  cyclic.  The  measured  data  shows  little 
sensitivity  to  negative  cyclic  for  this  condition.  Correlation 
results  for  the  first  two  chord  bending  harmonics  in  Figure  88 
show  the  same  trend,  with  the  analysis  correlating  well  for 
positive  cyclic  pitch  values  but  overpredicting  for  negative 
cyclic  pitch.  The  predicted  chord  bending  waveform  in  Figure 
89  shows  the  same  correlation  trends  as  the  flap  bending  wave¬ 
form.  The  predicted  amplitude  correlates  well  with  the  mea¬ 
sured  data  but  the  phase  of  the  predicted  first-harmonic 
content  is  shifted  45  degrees,  with  the  maximum  value  occurring 
at  120  degrees  azimuth. 
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Correlation  with  alternating  blade  root  torsion  is  shown  in 
Figures  90  and  91.  The  analysis  underpredicts  by  50  percent 
but  correlates  well  with  the  trend  of  alternating  torsion  with 
cyclic  pitch.  Correlation  with  the  torsional  waveform  in 
Figure  91  shows  the  analysis  underpredicts  the  first-harmonic 
content,  which  dominates  the  measured  waveform.  Note  that  the 
analysis  reflects  excitation  at  the  first  torsional  frequency 
(nearly  8  per  rev) ,  as  does  the  measured  data.  The  measured 
high-harmonic  content  on  the  retreating  side  of  the  disc  indi¬ 
cates  that  the  rotor  is  on  the  verge  of  stall  flutter  incep¬ 
tion.  The  analysis  predicts  this  trend  toward  stall  flutter 
but  with  appreciably  more  damping. 

Figures  92  through  94  are  the  correlation  results  for  steady 
rotor  in-plane  moments  in  transition.  The  correlation  is  pre¬ 
sented  as  the  gradient  of  pitching  and  yawing  moment  per 
degree  of  cyclic  and  also  the  gradient  of  total  hub  moment 
with  cyclic  pitch.  The  test  data  shows  a  linear  sensitivity 
of  pitching  and  yawing  moment  with  increasing  cyclic  and  the 
analysis  correlates  very  well  with  the  gradient  of  both  in¬ 
plane  moments.  Note  that  the  predicted  and  measured  pitching 
moments  are  more  sensitive  to  cyclic  than  yawing  moment. 

Cyclic  input  was  phased  with  maximum  cyclic  angle  occurring  at 
170  degrees  azimuth  so  that  the  flapping  phase  lag  response  of 
10  degrees  would  produce  nearly  pure  longitudinal  flapping  and 
therefore  a  large  pitching  moment  contribution.  The  analysis 
correctly  predicts  the  high  sensitivity  of  pitching  moment  and 
the  low  sensitivity  of  yawing  moment  to  cyclic  pitch.  The 
analysis  correlates  well  with  the  measured  total  hub  moment 
for  positive  cyclic  values,  but  overpredicts  hub  moment  sensi¬ 
tivity  for  negative  cyclic.  This  is  to  be  expected  since  the 
flap  bending  predictions  showed  the  same  trend  and  steady  hub 
moments  are  highly  dependent  upon  blade  flapping. 
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Figure  87.  1/3-Scale  LIT  Blade,  Predicted  and  Measured 

Alternating  Chord  Bending  Due  to  Cyclic  Pitch 
in  Transition. 
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Figure  89.  1/3-Scale  LIT  Blade 
Chord  Bending  Wavef 
Pitch  in  Transition 


ALTERNATING  TORSION  AT  .22R  (IN-LB) 
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Figure  90.  1/3-Scale  LIT  Blade,  Pr« 
Alternating  Torsion  Due 
Transition 
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Figure  93.  1/3-Scale  LIT  Blade,  Predicted  and  Measured 
Rate  of  Change  of  Steady  Rotor  Yawing  Moment 
With  Cyclic  in  Transition 
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STEADY  ROTOR  IN-PLANE  HUB  MOMENT  (FT-LB) 
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HELICOPTER  LOADS  CORRELATION 


INTRODUCTION 


Three  tandem-rotor  helicopter  flights  were  chosen  for  blade 
loads  correlation:  one  flight  of  the  CH-47  with  CH-47C 
blades;  one  flight  of  the  CH-47  with  advanced-geometry  blades 
(AGB) ;  and  one  flight  of  the  Model  347  with  CH-47C  blades. 
These  flights  were  chosen  since  they  offered  a  wide  range  of 
steady  level-flight  conditions  (particularly  rotor  loading  and 
advance  ratio),  as  well  as  having  varied  rotor  configurations. 
The  CH-47C  has  a  3-bladed,  constant-chord,  constant- thickness 
rotor  of  steel  and  fiberglass  construction;  the  CH-47  AGB  con¬ 
figuration  has  a  3-bladed,  tapered-chord ,  tapered- thickness 
rotor  of  fiberglass  construction.  Test  data  from  the  Model 
347  with  4  CH-47C  blades  provided  a  3-bladed/4-bladed  rotor 
comparison.  Wind  tunnel  test  data  from  the  CH-34  isolated 
full-scale  rotor  (Reference  31)  provided  airloads  data  for 
correlation.  Table  IV  is  a  summary  of  the  helicopters  and 
flight  conditions  used  for  correlation.  All  flight  test  data 
is  for  steady,  level,  forward  flight.  Table  V  is  a  summary 
of  the  helicopter  test  data  used  for  correlation. 

CORRELATION  RESULTS 


Pitch  Link  Loads 


CH-47  and  Model  347  tandem-rotor  helicopter  test  data  were 
used  for  pitch  link  load  correlation.  The  test  data  includes 
pitch  link  loads  for  predominately  unstalled  and  predominately 
stalled  rotor  conditions,  the  deepest  stalled  condition  being 
the  aft  rotor  of  the  CH-47C,  flight  X-83,  at  133  knots.  The 
aft  rotor  for  this  condition  was  operating  at  p  =  0.306  and 
Ct/c  =  0.1)47,  which  is  the  highest  C >p/o  encountered  for  the 
total  helicopter  correlation. 

The  pitch  link  loads  measured  from  flight  test  data  are  pre¬ 
sented  as  one-half  of  the  peak-to-peak  value  for  one  rotor 
cycle.  Most  test  data  points  are  the  average  value  of  25  to 
120  rotor  cycle  values.  The  symbol  represents  the  arithmetic 
average  value  of  all  the  cycles,  while  the  vertical  line 
through  the  symbol  represents  the  highest  and  lowest  values  of 
all  the  cycles  in  that  run.  This  is  a  convenient  means  of 
plotting  flight  test  data  where  data  scatter'  is  expected  due 
to  pilot  handling  or  turbulent  conditions  but  it  also  can  be 
misleading,  especially  in  or  near  a  stalled  condition.  In 
some  cases  the  value  of  the  most  typical  cycle  (the  value  that 
occurred  most  often  during  the  run)  is  not  the  same  as  the 
arithmetic  mean  of  the  data.  A  good  example  of  this  is  the 
aft  rotor  pitch  link  load  of  CH-47C,  flight  X-83,  at  123  knots, 
shown  in  Figure  95.  The  pitch  link  load  test  data  at  this 
speed  has  a  scatter  from  2,070  to  600  pounds  with  an  arithmetic 
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mean  of  1,540  pounds.  But,  as  shown  in  Figure  96,  the  value 
of  the  cycle  that  occurs  most  often  is  between  1,800  and  1,900 
pounds,  an  increase  of  300  pounds  over  the  arithmetic  mean. 

This  indicates  that  the  aft  rotor  was  normally  in  a  stalled 
condition,  shown  by  the  large  number  of  occurrences  of  high 
pitch  link  loads,  but  due  to  pilot  techniques  or  turbulent  con¬ 
ditions  the  rotor  briefly  went  into  an  unstalled  condition. 
Therefore,  caution  must  be  taken  not  to  weigh  the  arithmetic 
mean  too  heavily  when  correlating  pitch  link  loads  in  or  near 
a  stalled  condition. 

Figures  95  through  107  are  the  aft  rotor  pitch  link  load  cor¬ 
relation  results  for  CH-47C  and  Model  347  tandem-rotor  heli¬ 
copter  test  data.  Figures  95,  97,  and  98  show  the  analysis 
correlates  well  with  the  pitch  link  load  data  with  increasing 
airspeed;  the  prediction  for  CH-47C  flight  X-83  (Figure  95)  is 
especially  good  in  predicting  stall  inception  and  the  growth 
of  loads  as  airspeed  increases.  As  mentioned  before,  the  most 
typical  value  of  pitch  link  load  is  significantly  different 
than  the  arithmetic  mean  load  for  this  flight  in  stall.  The 
most  typical  value  is  plotted  for  the  two  highest  airspeeds 
for  a  more  meaningful  correlation  with  the  analysis. 

As  shown  by  the  correlation  results  in  Figure  95,  the  analysis 
with  uniform  downwash  correlates  very  well  with  pitch  link 
loads  in  deep  stall  at  123  and  133  knots.  Uniform  downwash  was 
used  primarily  in  the  analysis  for  two  reasons.  Use  of  uni¬ 
form  downwash  instead  of  nonuniform  downwash  results  in  sub¬ 
stantial  savings  in  computer  time.  The  second  reason  for 
using  uniform  downwash  was  to  avoid  the  problems  associated 
with  a  high  Cy/o  condition  such  as  that  shown  in  Figure  95  at 
123  and  133  knots.  For  a  condition  such  as  this,  predicted 
pitch  link  loads  when  nonuniform  downwash  is  used  in  the 
analysis  have  been  found  to  be  sensitive  to  rotor  interference. 
Presumably,  blade  moment  stall,  and  therefore  pitch  link  loads, 
are  highly  dependent  upon  the  analysis  description  of  blade 
tip  vortex  proximity  to  a  passing  blade.  As  shown  in  the 
figure,  after  stall  inception  at  123  and  133  knots  where 
moment  stall  is  significant,  the  nonuniform  downwash  predic¬ 
tions  are  erratic.  Since  use  of  uniform  downwash  in  the 
analysis  correlates  well  with  measured  pitch  link  loads  in 
and  out  of  stall,  uniform  downwash  is  recommended  for  pitch 
link  load  prediction. 


Figures  99,  100,  and  101  are  the  aft  rotor  pitch  link  load 
waveform  correlation  results  for  CH-47C  flight  X-83.  The 
three  highest  airspeeds  in  the  flight  are  shown  and  the  test 
data  at  each  airspeed  was  chosen  from  a  cycle  near  the  most 
typical  alternating  value  for  that  airspeed.  At  111  knots  the 
analysis  shows  definite  stall  indications,  as  does  the  test 
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data,  and  at  123  and  133  knots  the  analysis  predicts  the  large 
spike  on  the  retreating  side  of  the  disc  plane.  The  analysis 
not  only  correlates  well  with  the  amplitude  of  the  spike  but 
also  predicts  within  30  degree  ■>  the  phasing  of  the  spike  in 
the  rotor  azimuth. 

For  the  aft  rotor  pitch  link  load  correlation  for  the  CH-47 
with  advanced-geometry  blades.  Figure  97  shows  the  analysis 
correlates  well  with  increasing  airspeed.  As  mentioned  before, 
the  aft  rotor  for  this  flight  attained  a  maximum  C-j/a  of 
0.C890,  a  predominately  unstalled  condition.  The  analysis 
correctly  predicts  this  predominately  unstalled  trend  as  shown 
by  the  smooth  increase  in  predicted  load  with  airspeed. 

Figures  102,  103,  and  104  show  the  waveform  correlation  results 
at  116,  125,  and  131  knots.  The  test  data  is  predominately 
1-per-rev,  indicating  that  the  rotor  is  operating  below  stall. 
The  analysis  correctly  predicts  the  amplitude  of  the  1-per-rev 
load  and  the  absence  of  a  spike  in  the  prediction  shows  that 
the  analysis  correctly  predicts  that  the  rotor  is  operating 
out  of  stall.  The  higher  harmonics  in  the  test  data  are  due 
to  rotor  interference  and  cannot  be  reflected  in  the  analysis 
results  when  uniform  downwash  is  used.  These  correlation  re¬ 
sults  are  important  since  they  demonstrate  the  capability  of 
the  analysis  to  predict  aft  rotor  pitch  link  loads  for  blades 
that  are  tapered  in  chord  and  thickness. 

Figure  98  shows  the  aft  rotor  pitch  link  load  correlation  re¬ 
sults  with  airspeed  for  the  Model  347  with  CH-47C  blades.  The 
analysis  correlates  very  well  with  the  test  data,  with  both 
test  data  and  analysis  showing  the  rotor  is  operating  predom¬ 
inately  unstalled.  The  waveform  correlation  in  Figures  105, 
106,  and  107  show  that  both  test  data  and  prediction  are 
primarily  1-per-rev  and  correlation  is  good  in  amplitude  and 
phase.  The  test  data  also  contains  2-per-rev  and  higner  har¬ 
monics  due  to  rotor  interference.  Use  of  uniform  downwash 
neglects  higher  harmonic  forcing  due  to  interference  and  the 
predicted  waveform  does  not  show  any  significant  higher  har¬ 
monics.  But  since  the  first  harmonic  dominates  the  measured 
pitch  link  load,  this  correlates  well  with  the  test  data. 

Figures  108  through  119  are  the  forward  rotor  pitch  link  load 
correlation  results.  For  the  CH-47C  forward  rotor  in  Figure 
108,  the  prediction  increases  sharply  at  123  knots.  The  analy¬ 
sis  is  showing  premature  stall  inception,  as  demonstrated  in 
the  pitch  link  load  waveforms  in  Figures  109,  110,  and  111. 

At  111  and  123  knots,  both  test  data  and  prediction  are  pri¬ 
marily  1-per-rev  with  the  analysis  showing  a  slight  stall 
indication  in  the  fourth  quadrant  at  123  knots.  At  133  knots 
the  analysis  shows  definite  stall  characteristics  as  demon¬ 
strated  by  the  large  spikes  in  the  fourth  quadrant.  The  test 
data  at  133  knots  is  showing  a  slight  indication  of  stall  in 
the  fourth  quadrant  but  is  still  primarily  1-per-rev.  The 
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rotor  is  operating  at  a  Gp/a  of  0.0937  ,  which  is  on  the  verge 
of  stall  inception,  and  the  analysis  has  predicted  stall  to 
occur  prematurely  by  about  10  knots.  Figures  112  through  119 
are  the  forward  rotor  pitch  link  load  correlation  results  for 
the  CH-47  AGB  and  the  Model  347  with  CH-47C  blades.  The  pre¬ 
dictions  and  test  data  for  both  aircraft  flights  show  the 
forward  rotor  is  unstalled  and  the  waveforms  show  that  the 
pitch  link  load  is  primarily  1-per-rev.  The  analysis  overpre¬ 
dicts  the  growth  of  pitch  link  load  with  airspeed  for  the 
CH-47  AGB  (Figure  112)  and  the  analysis  correlates  well  with 
the  growth  of  pitch  link  load  for  the  Model  347  (Figure  116) . 

As  in  the  case  for  the  CH-47C  forward  rotor,  the  pitch  link 
load  waveforms  for  both  of  these  aircraft  show  higher  harmonics 
due  to  interference.  But  note  that  the  first  harmonic  domi¬ 
nates  the  amplitude  of  the  measured  waveforms.  The  capability 
of  the  analysis  to  predict  pitch  link  loads  before  stall 
(primarily  1-per-rev  loads)  and  in  stall  (impulsive  loading  in 
the  fourth  quadrant)  demonstrates  that  the  use  of  uniform 
downwanh  in  the  analysis,  thereby  neglecting  rotor  inter¬ 
ference,  is  justified. 
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Figure  98.  Model  347  With  CH-47C  Blades /  Aft  Rotor 
Predicted  and  Measured  Alternating  Pitch 
LinK  Loads  With  Airspeed 
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Figure  100.  CH-47C  Aft  Rotor,  Predicted 
Pitch  Link  Waveform  at  123 
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ROTOR  AZIMUTH  (DEG) 
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GW  =  39,000  LB 


Figure  109.  CH-47C  Forward  Rotor,  Predicted  and  Measured  Pitch 
Link  Load  Waveform  at  111  Knots 


AIRCRAFT  B361 


Figure  110.  CH-47C  Forward  Rotor,  Predicted  and  Measured  Pitch 
Link  Load  Waveform  at  123  Knots 
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Figure  112.  CH-47C  With  Advanced-Geometry  Blades,  Forward 
Rotor  Predicted  and  Measured  Alternating  Pitch 
Link  Loads  With  Airspeed 


Figure  113.  CH-47C  With  Advanced-Geometry  Blades,  Forward  Rotor 

Predicted  and  Measured  Pitch  Link  Waveform  at  117  Knots 
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Figure  116.  Model  347  With  CH-47C  Blades,  Forward  Rotor 
Predicted  and  Measured  Alternating  Pitch 
Link  Loads  With  Airspeed 
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Flap  Bending 


Figures  120  through  131  are  the  correlation  results  for  aft 
rotor  alternating  flap  bending.  The  test  data  presented  in 
these  figures  is  for  the  CH-47  and  Model  347  (both  tandem- 
rotor  configuration) .  The  alternating  bending  moments  shown 
are  one-half  of  the  peak-to-peak  value  for  one  rotor  cycle. 

The  symbol  represents  the  arithmetic  average  value  of  all  the 
cycles  in  the  run  while  the  vertical  line  through  the  symbol 
represents  the  highest  and  lowest  value  of  all  the  cycles  for 
that  run. 

The  analysis  with  uniform  downwash  shows  a  consistent  trend  to 
increase  mid-span  and  outboard  alternating  flap  bending  with 
airspeed.  This  predicted  trend  is  consistent  with  the  test 
data,  but  the  analysis  overpredicts  with  increasing  airspeed. 
Model  347  flight  313  shows  the  best  outboard  bending  moment 
correlation,  while  inboard  bending  moment  correlation  for  both 
the  Model  347  and  the  CH-47  tends  to  be  underpredicted.  The 
primary  reason  for  overprediction  outboard  of  mid-span  and 
underprediction  inboard  is  the  use  of  uniform  downwash  in  the 
analysis.  When  uniform  downwash  is  used,  downwash  strength  on 
the  retreating  side  of  the  rotor  azimuth  is  overpredicted, 
thereby  creating  a  smaller  blade  section  angle  of  attack  and 
reduced  lift  on  the  retreating  side  of  the  rotor  disc.  This 
results  in  overprediction  of  1-per-rev  forcing  which  in  turn 
results  in  overprediction  of  1-,  2-,  and  3-per-rev  flap  bend¬ 
ing  moments.  Figure  120  shows  why  1-,  2-,  and  3-per-rev  flap 
bending  are  overpredicted  when  uniform  downwash  is  used  in  the 
analysis.  This  figure  shows  the  predicted  flap  mode  shapes  in 
a  vacuum  for  the  CH-47C  rotor  blade.  Note  that  the  first  flex¬ 
ible  mode  with  a  frequency  ratio  of  2.57  at  230  rpm  is  roughly 
symmetrical  about  its  peak  bending  moment  at  mid-span.  Any 
1-per-rev  forcing  will  result  in  1-,  2-,  and  3-per-rev  flap 
bending  moments  with  a  radial  distribution  similar  to  the  first 
flexible  mode.  Note  that  all  the  flap  bending  predictions 
with  uniform  downwash  have  this  same  radial  distribution,  re¬ 
flecting  dominant  excitement  of  the  first  flexible  mode. 

Use  of  nonuniform  downwash  in  the  analyses  more  closely  pre¬ 
dicts  blade  lift  radial  distribution  around  the  azimuth  and 
also  introduces  blade  tip  vortices  lying  in  or  near  the  disc 
plane.  For  the  forward  rotor  of  a  tandem-rotor  configuration, 
blade  vortices  are  predominantly  due  to  that  of  a  preceding 
blade  of  the  same  rotor  (called  blade-to-blade  interference)  , 
while  for  the  aft  rotor  of  a  tandem-rotor  helicopter  blade 
vortex  interference  is  predominantly  due  to  tip  vortices  from 
the  forward  rotor  trailing  downstream  (called  rotor-on-rotor 
interference) .  Vortex  strikes  or  near-strikes  tend  to  impart 
impulsive  forcing  to  the  blade  due  to  the  high  velocities 
associated  with  tip  vortices,  therefore  exciting  the  higher 
flap  bending  modes.  These  modes  (primarily  second  flexible 
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through  fourth  flexible  modes)  are  lightly  damped  and  therefore 
contribute  significantly  to  the  flap  bending  moment  when  vortex 
interference  is  involved.  Note  from  Figure  120  that  the  sec¬ 
ond,  third,  and  fourth  flexible  modes  for  the  CH-47C,  4.75  o/O, 
7.90  uj/fi ,  and  12.19  ui/ft,  respectively,  all  peak  at  10  to  13 
percent  radius  and  80  to  85  percent  radius.  This  indicates 
that  rotor  blade  flap  bending  moments,  when  operating  with 
significant  vortex  disturbances,  would  have  a  large  outboard 
and  inboard  peak  mainly  due  to  higher  harmonics  (4  to  12  per 
rev) .  Also  note  that  the  third  flexible  mode  suppresses  bend¬ 
ing  moments  near  mid-span  due  to  the  first  flexible  mode,  thus 
indicating  that  use  of  nonuniform  downwash  in  the  analysis 
would  tend  to  suppress  first-,  second-,  and  third-harmonic 
bending  at  mid-span.  The  analysis  was  run  with  nonuniform 
downwash  for  the  highest  speed  condition  for  each  helicopter 
to  show  the  impact  of  higher  harmonics  due  to  vortex  inter¬ 
ference  on  predicted  blade  flap  bending.  The  results,  shown 
in  Figures  122,  123,  127,  and  131,  reflect  the  outboard  and 
inboard  peak  and  correlate  better  with  the  flap  bending  test 
data  than  the  predictions  with  uniform  downwash.  A  sharper 
inboard  peak  and  outboard  peak  would  be  predicted  with  non- 
uniform  downwash  if  the  analysis  had  12-harmonic  capability 
instead  of  10  harmonics.  Note  that  the  fourth  flexible  mode 
in  Figure  120  has  a  frequency  ratio  of  12.19  per  rev;  and 
therefore,  12-per-rev  bending  moment  contributions  due  to  this 
mode  at  10  and  80  percent  span  will  not  be  reflected  in  the 
analysis . 

Figures  132  through  142  are  the  correlation  results  for  the 
forward  rotor  alternating  flap  bending  on  the  tandem-rotor 
helicopters.  The  uniform  downwash  predictions  for  these  data 
follow  the  same  trend  as  the  uniform  downwash  predictions  for 
the  aft  rotor.  Mid-span  and  outboard  flap  bending  increases 
with  increasing  airspeed  and  overpredicts.  As  with  the  aft 
rotor  predictions,  use  of  uniform  downwash  overpredicts  1-per- 
rev  forcing  and  neglects  vortex  interference.  The  only  dif¬ 
ference  is  that  forward  rotor  interference  is  mainly  blade-to- 
blade  while  aft  rotor  interference  is  predominantly  rotor-on- 
rotor.  Figures  134,  138,  and  142  show  nonuniform  downwash 
predictions  for  the  highest  speed  case  for  each  helicopter 
rotor.  The  predictions  with  nonuniform  downwash  correlate 
better  than  those  with  uniform  downwash,  reflecting  the  out¬ 
board  ana  inboard  peak  due  to  higher  harmonics  induced  by 
interference . 
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Figure  121.  CH-47C  Aft  Rotor,  Predicted  and  Measured 

Alternating  Flap  Bending  Moment  at  V  =  111  Knots 
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ALTERNATING  FLAP  BENDING  MOMENT  X  10  (IN-LB) 
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Figure  122.  CH-47C  Aft  Rotor,  Predicted  and  Measured 

Alternating  Flap  Bending  Moment  at  V  =  123  Knots 
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Figure  124.  CH-47C  With  Advanced-Geometry  Blades,  Aft 

Rotor  Predicted  and  Measured  Alternating  Flap 
Bending  Moment  at  100  Knots 
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ALTERNATING  FLAP  BENDING  MOMENT  xlO  3 (IN-LB) 


Figure  125.  CH-47C  With  Advanced-Geometry  Blades,  Aft 

Rotor  Predicted  and  Measured  Alternating  Flap 
Bending  Moment  at  117  Knots 
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Figure  130.  Model  347  With  CH-47C  Blades,  Aft  Rotor  Pre¬ 
dicted  and  Measured  Alternating  Flap  Bending 
Moment  at  V  =  140  Knots 
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Model  347  With  CH-47C  Blades,  Aft  Rotor  Pre¬ 
dicted  and  Measured  Alternating  Flap  Bending 
Moment  at  V  =  152  Knots 


Figure  132.  CH-47C  Forward  Rotor,  Predicted  and  Measured 

Alternating  Flap  Bending  Moment  at  V  =  111  Knots 
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136.  CH-47C  With  Advanced-Geometry  Blades#  Forward 
Rotor  Predicted  and  Measured  Alternating  Flap 
Bending  Moment  at  117  Knots 
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Chord  Bending 


Figures  143  through  154  are  the  correlation  results  for  heli¬ 
copter  aft  rotor  chord  bending.  CH-47  and  Model  347  flight 
test  data  was  used  for  correlation.  Uniform  downwash  was  used 
primarily  in  the  analysis  for  chord  bending  predictions  to  re¬ 
duce  computer  time.  However,  one  nonuniform  downwash  predic¬ 
tion  was  made  for  each  flight  to  show  the  impact  of  downwash 
representation  on  chord  bending. 

The  aft  rotor  correlation  results  show  that  the  analysis  con¬ 
sistently  underpredicts  when  uniform  downwash  is  used.  This 
is  to  be  expected  due  to  the  lag  damper  representation  in  the 
analysis.  The  lag  damper  in  the  analysis  is  an  ideal  torsional 
viscous  damper  and  therefore  the  energy  absorbed  by  the  damper 
is  dependent  upon  the  lag  velocity.  The  lag  damper  on  the 
CH-47  and  Model  347  is  a  torsional  viscous  damper  with  a  con¬ 
stant  coulomb  (preload)  damper  superimposed.  As  shown  in 
Figure  155,  the  CH-47C  lag  damper  has  a  large  preload  which  is 
independent  of  lag  velocity  except  for  valve  leakage.  There¬ 
fore,  to  equate  the  lag  damper  value  in  the  analysis  to  the 
actual  lag  damper,  a  1-per-rev  lag  angle  is  assumed  and  the 
lag  damper  value  in  the  analysis  is  determined  by  equating 
energy.  The  energy  absorbed  by  the  lag  damper  in  the  analysis 
is  forced  to  be  equal  to  that  absorbed  by  the  actual  lag  damper 
by  the  following  equation: 


C  =  +  C 

^ANALYSIS  "“^MAX 

where  Cr  ,  „  „  =  viscous  torsional  damping  rate  in  the 

^ANALYSIS  analysis 

C  =  actual  viscous  torsional  damping  rate 

Mp  =  preload  moment  in  damper 

oj  =  lag  natural  frequency 

^MAX  =  l“Per_rev  amplitude  of  lagging  motion 


This  method  of  representing  a  lag  damper  with  a  preload  would 
yield  a  fair  approximation  of  the  amplitude  of  the  lag  bend¬ 
ing  moment  if  the  lag  motion  was  nearly  all  1  per  rev.  But 
even  then  i_he  phasing  would  be  incorrect  since  the  force  ex¬ 
erted  by  the  actual  damper  preload  (acting  like  a  spring)  is 
90  degrees  out  of  phase  with  the  force  exerted  by  the  viscous 
torsional  damper.  In  addition,  the  use  of  uniform  downwash  in 
the  analysis  partially  neglects  1-,  2-,  and  3-per-rev  airload¬ 
ing  induced  by  blade  tip  vortex  proximity  to  a  passing  blade, 
thus  in  part  neglecting  1-,  2-,  and  3-per-rev  lagging  due  to 
induced  drag.  Therefore  use  of  nonuniform  downwash  in  the 
analysis  to  account  for  airloading  due  to  vortex  proximity 
should  increase  1-,  2-,  and  3-per-rev  lag  motion  and  yield 
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higher  chord  bending'  motion,  even  though  the  rate  at  which  the 
damper  is  doing  work  in  the  analysis  is  based  on  1-per-rev  lag 
motion . 

As  shown  in  Figures  144,  145,  and  149,  use  of  nonuniform  down- 
wash  in  the  analysis  substantially  increases  predicted  chord 
bending  and  the  results  correlate  well  with  bottom-of-scatter 
data.  It  should  be  noted  that  the  effect  of  rotor  interfer¬ 
ence  on  chord  bending  moment  is  limited  to  excitation  of  the 
first  flexible  mode,  primarily  due  to  the  high  frequency  ratio 
of  the  second  flexible  mode  as  shown  in  Figure  156.  Comparing 
this  figure  to  the  flap  bending  moment  mode  shapes  in  Figure 
120,  chord  bending  for  the  CH-47C  has  only  one  flexible  mode 
with  a  frequency  between  1  and  12  per  rev,  while  flap  bending 
has  four  flexible  modes  up  to  12  per  rev.  Flight  test  experi¬ 
ence  has  shown  that  the  frequency  spectrum  of  rotor  interfer¬ 
ence,  as  reflected  in  bending  moments,  is  between  1  and  12  per 
rev.  Thus  the  second  flexible  mode  is  lightly  forced  and 
high-harmonic  (6-  through  12-per-rev)  chord  bending  is  low  in 
amplitude.  The  aft  rotor  chord  bending  correlation  for  the 
Model  347  in  Figures  150  through  153  shows  the  analysis  under¬ 
predicts  by  a  large  margin  even  with  nonuniform  downwash  in 
the  analysis. 

Comparing  Figures  145  and  152,  the  test  data  for  the  Model  347 
is  1.5  times  as  large  as  the  CH-47C  at  50  percent  span  at 
nearly  the  same  airspeed.  The  high  chord  bending  moments  on 
the  Model  347  are  due  to  a  resonant  condition  between  4-per- 
rev  in-plane  forcing  and  the  third  and  fourth  flexible  modes 
of  the  coupled  rotor/drive  system  which  are  nearly  4  per  rev. 
This  condition  is  discussed  extensively  in  Reference  32.  As 
shown  in  Figure  154,  alternating  chord  bending  is  dominated 
by  the  fourth-harmonic  contribution.  Since  the  analysis  does 
not  have  the  capability  to  input  hub  motion  and  simulate  the 
rotor/drive  system  coupling,  the  predicted  4-per-rev  chord 
bending  is  much  lower  than  measured  test  data. 

Figures  157  through  168  are  the  correlation  results  for  heli¬ 
copter  forward  rotor  chord  bending.  These  figures  show  pre¬ 
diction  trends  similar  to  that  for  the  aft  rotor'.  The  analysis 
underpredicts  when  uniform  downwash  is  used  due  to  the  lack  of 
higher-harmonic  lagging  induced  by  blade-to-blade  interference. 
When  nonuniform  downwash  is  used  in  the  analysis  the  chord 
bending  prediction  increases  substantially  and  correlates 
better  with  the  test  data  (Figures  159  and  163) .  The  higher- 
harmonic  lag  motions  are  now  included  and  the  lag  damper  is 
absorbing  more  energy  per  rotor  cycle,  thus  increasing  chord- 
wise  bending  moments.  Note  the  small  increase  in  predicted 
chord  bending  in  Figure  163  when  nonuniform  downwash  is  used; 
this  indicates  that  predicted  interference  effects  are  small 
and  the  predicted  lag  motion  is  primarily  1  per  rev.  The  Model 
347  forward  rotor  correlation  results  in  Figures  164  through 
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167  show  the  analysis  to  be  underpredicting  by  a  large  margin 
with  uniform  and  nonuniform  downwash.  As  shown  in  Figure  168, 
the  test  data  has  a  large  4-per-rev  contribution  due  to  the 
rotor/drive  system  coupling,  as  discussed  earlier,  and  the 
analysis  does  not  have  the  capability  to  simulate  drive  system 
coupling . 

In  summary,  the  chord  bending  predictions,  especially  for  the 
aft  rotor  with  nonuniform  downwash,  correlate  well  with  the 
test  data  when  it  is  considered  that  the  lag  damper  represen¬ 
tation  in  the  analysis  does  not  have  the  capability  to  simulate 
damper  preload.  It  is  extremely  difficult  to  represent  the  lag 
damper  effect  on  chord  bending  by  equating  energy  as  discussed 
earlier,  especially  when  the  lag  motion  contains  higher  har¬ 
monics.  Inputting  a  lag  damper  rate  into  the  analysis  for  a 
preloaded-type  damper  is  at  best  an  approximation.  A  para¬ 
metric  study  of  the  effect  of  lag  damper  value  on  predicted 
chord  bending  would  be  helpful,  but  the  study  would  be  compli¬ 
cated  by  the  dependence  of  lagwise  motion  (amplitude,  harmonic 
contribution,  and  phasing)  on  specific  aircraft  and  flight  con¬ 
ditions.  Interference,  which  contributes  to  higher-harmonic 
lagwise  motion,  is  highly  dependent  upon  tandem-rotor  separa¬ 
tion,  C-p/a ,  advance  ratio,  and  tip  path  plane  tilt.  The  most 
promising  solution  for  chord  bending  prediction  is  to  treat 
the  lag  damper  in  the  analysis  as  a  torsional  viscous  damper 
and  a  constant  coulomb  damper  in  the  form  of  two  table  lookup 
routines.  One  table  would  account  for  the  preload  (constant 
coulomb  damper)  shown  in  Figure  155  and  the  other  table  would 
account  for  the  viscous  damper.  With  this  method,  the  force 
exerted  by  the  damper  when  the  blade  is  in  a  specific  azimuth 
position  will  be  entirely  dependent  upon  aerodynamic  forcing 
and  dynamic  response  at  that  instant  instead  of  being  depen¬ 
dent  upon  an  average  lag  damper  value  for  the  entire  rotor 
cycle . 
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Figure  146 


CH-  47C  With  Advanced-Geometry  Blades,  Aft 
Rotor  Predicted  and  Measured  Alternating 
Chord  Bending  Moment  at  100  Knots 
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AIRCRAFT  347-B164 
FLIGHT  313-23 
GW  42,079  LB 
CG  6  IN.  AFT 
Hd  1,607  FT 
Nr  222  RPM 
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AIRCRAFT  347-B164 
FLIGHT  313-26 
GW  =  42,000  LB. 

CG  =  6  IN.  AFT 
H-  =  1,600  FT. 

Nr  =  220  RPM 
V  =  152  KTS . 


0  180  360 


ROTOR  AZIMUTH  (DEG) 


Figure  154.  Model  347  With  CH-47C  Blades,  Aft  Rotor 

Predicted  and  Measured  Chord  Bending  Wave¬ 
form  at  152  Knots 
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ALTERNATING  CHORD  BENDING  MOMENT  X  10  (IN-LB) 


AIRCRAFT  B36 1 
FLIGHT  83-72 
GW  38,915  LBS 
CG  4  IN.  AFT 
Hd  7200  FT 
Nr  230  RPM 


Figure  157.  CH-47C  Forward  Rotor,  Predicted  and  Measured 
Alternating  Chord  Bending  Moment  at  V  =  111 
Knots 
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ALTERNATING  CHORD  BENDING  MOMENT  X  10  '  (IN-LB) 
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AIRCRAFT  B361 
FLIGHT  272-5 
GW  =  46,000  LB 
CG  =  4  IN.  AFT 
Hd  =  6,000  FT 
Nr  =  245  RPM 
V  =  100 


AIRCRAFT  B36 1 
FLIGHT  272  -  7 
GW  =  46000  LB. 
CG  =  4  IN.  AFT 
H,  =  6000  FT. 
Np  =  245  RPM 
V  =  117 


r/R 

Figure  161.  CH-47C  With  Advanced-Geometry  Blades,  Forward 
Rotor  Predicted  and  Measured  Alternating  Chord 
Bending  Moment  at  117  Knots 
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ALTERNATING  CHORD  BENDING  MOMENT  xlO_J  (IN-LB) 


AIRCRAFT  B361 
FLIGHT  272-9 
GW  =  46000  LB. 
CG  =  4  IN.  AFT 
H  ,  =  6000  FT. 

=  245  RPM 
V  —  126 
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Figure  162. 


CH-47C  With  Advanced-Geomi 
Rotor  Predicted  and  Measu: 
Bendina  Moment  at  126  Kno 


AIRCRAFT  B361 
FLIGHT  272  -  10 
GW  =  46000  LB. 
CG  =  4  IN.  AFT 
H,  =  6000  FT. 

=  245  RPM 
V  =  132 
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AIRCRAFT  347-B164 
FLIGHT  313-19 
GW  42,274  LB 
CG  6  IN.  AFT 
Hd  1,687  FT 
Nr  222  RPM 
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Figure  164 


Model  347  With  CH-47C  Blades,  Forward  Rotor 
Predicted  and  Measured  Alternating  Chord 
Bending  Moment  at  V  =  87  Knots 
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AIRCRAFT  347-B164 
FLIGHT  313-23 
GW  42,079  LE 
CG  6  IN.  AFT 
Hd  1,607  FT 
Nr  222  RPM 
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Figure  165 


Model  347  With  CH-47C  Blades,  Forward  Rotor 
Predicted  and  Measured  Alternating  Chord 
Bending  Moment  at  V  =  122  Knots 


ALTERNATING  CHORD  BENDING  MOMENT  X  10  (IN-LB) 


AIRCRAFT  347-B164 
FLIGHT  313-25 
GW  42,008  LB 
CG  6  IN.  AFT 
Hd  1,410  FT 
Nr  219  RPM 


TEST 


THEORY 


Figure  166. 


Model  347  With  CH-47C  Blades,  Forward  Rotor 
Predicted  and  Measured  Alternating  Chord 
Bending  Moment  at  V  =  140  Knots 


ALTERNATING  CHORD  BENDING  MOMENT  X  10  (IN-LB) 


AIRCRAFT  347-B164 
FLIGHT  313-26 
GW  41,744  LB 
CG  6  IN.  AFT 
Hd  1,657  FT 
Nr  218  RPM 


THEORY 


Figure  167. 


Model  347  With  CH-47C  Blades,  Forward  Rotor 
Predicted  and  Measured  Alternating  Chord 
Bending  Moment  at  V  =  152  Knots 
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AIRCRAFT  347-B164 
FLIGHT  313-26 
GW  =  42,000  LB. 

CG  =  6  IN.  AFT 
H.  =  1,600  FT. 

Nr  =  220  RPM 
V  =  152  KTS . 
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Figure  168.  Model  347  With  CH-47C  Blades,  Forward  Rotor 
Predicted  and  Measured  Chord  Bending  Wave¬ 
form  at  152  Knots 
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Airloads 


Airloads  correlation  was  performed  with  measured  differential 
chordwise  pressure  distributions  obtained  from  a  test  of  a 
CH-34  full-scale  rotor  conducted  in  the  NASA-Ames  40-foot  by 
80-foot  wind  tunnel.  Blade  physical  and  geometric  properties 
and  the  airloads  test  data  were  obtained  from  Reference  31. 

The  predicted  and  measured  airloads  are  presented  as  lift-per- 
unit-span  time  histories  for  one  rotor  cycle.  Correlation 
with  uniform  downwash  is  made  at  5  radial  stations  spaced  from 
mid-span  to  the  blade  tip.  Test  data  for  three  airspeeds  (110, 
150,  and  175  knots)  at  a  constant  shaft  tilt  of  5  degrees  for¬ 
ward  was  used  for  correlation  to  show  the  effect  of  airspeed. 

In  addition,  test  data  at  a  110-knot  condition  with  10  degrees 
forward  shaft  tilt  was  used  to  show  the  effect  of  shaft  tilt. 

Results  of  the  correlation,  presented  in  Figures  169  through 
188,  show  the  analysis  correlates  with  the  steady  and  peak-to- 
peak  measured  airloads.  The  correlation  is  particularly  good 
at  0.55R  (Figures  169,  174,  179,  and  184)  and  0.75R  (Figures 
170,  175,  180,  and  185) .  At  85  percent  and  95  percent  span  at 
175  knots  (Figures  186  and  187) ,  the  analysis  underpredicts 
the  peak-to-peak  airloads.  Underprediction  is  probably  due  to 
a  strong  blade-to-blade  interference  effect  which  cannot  be 
simulated  when  uniform  downwash  is  used  in  the  analysis.  Note 
that  the  analysis  overpredicts  airloads  by  a  large  margin  at 
99  percent  span  for  all  airspeeds  (Figures  173,  178,  183,  and 
188)  .  The  effect  of  the  tip  vortex  on  airloads  is  very  strong 
over  the  outboard  5  percent  of  the  blade  and  washout  of  the 
blade  tip  cannot  be  simulated  with  uniform  downwash.  Figures 
169  through  173  can  be  compared  with  Figures  174  through  178 
for  correlation  with  forward  shaft  tilt  at  110  knots.  The 
analysis  correlates  well  at  0.55R  and  0.75R  but  overpredicts 
at  0.85R  and  0.95R,  as  well  as  at  the  blade  tip.  Overpredic¬ 
tion  at  0.85R  and  0.95R  is  surprising  since  increasing  forward 
shaft  tilt  lessens  blade-to-blade  interference,  therefore 
suggesting  that  uniform  downwash  predictions  should  be  better 
at  higher  shaft  tilts.  Just  the  opposite  is  shown  to  be  true 
at  0.95R  (Figures  172  and  177);  correlation  is  better  at  the 
lower  shaft  tilt. 

A  comparison  of  the  predicted  and  measured  azimuthal  variation 
of  airloads  shows  that  both  the  analysis  and  measured  data 
show  a  large  decrease  in  lift  on  the  advancing  blade.  How¬ 
ever,  the  analysis  shows  a  phase  shift  of  about  30  degrees — 
the  predicted  minimum  lift  occurring  30  degrees  before  the 
measured  minimum  at  120  degrees.  The  uniform  downwash  repre¬ 
sentation  in  the  analysis  is  not  sufficient  to  describe  the 
downwash  nonuniformity  in  the  region  of  rapidly  changing  vortex 
strength.  In  this  region  of  the  azimuth,  the  measured  airloads 
also  show  evidence  of  tip  vortex  interference  from  the  pre¬ 
ceding  blade.  Note  in  Figures  171  and  172  the  high  harmonic 
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airloading  at  about  90  degrees  azimuth  for  5  degrees  shaft  tilt 
and  110  knots.  As  mentioned  before,  vortex  interference  is 
expected  to  decrease  as  forward  shaft  tilt  is  increased,  and 
this  trend  is  evident  when  Figures  171  and  172  are  compared 
with  Figures  176  and  177.  In  the  last  two  figures,  the  forward 
shaft  tilt  has  been  increased  to  9  degrees,  and  the  high- 
harmonic  airloading  at  90  degrees  azimuth  is  no  longer  present. 

In  summary,  the  good  peak-to-peak  airloads  correlation  indi¬ 
cates  that  the  analysis  with  uniform  downwash  is  adequate  for 
predicting  peak-to-peak  bending  moments  on  a  single-rotor 
helicopter.  However,  the  higher  harmonic  content  of  bending 
moments  will  not  be  reflected  in  the  peak-to-peak  value  if 
blade-to-blade  interference  effects  are  significant  (near  zero 
or  aft  shaft  tilt) .  It  also  should  be  remembered  that  rotor- 
on-rotor  interference  effects  are  significant  even  at  large 
forward  shaft  tilts  for  tandem-rotor  configurations  and  there¬ 
fore,  use  of  uniform  downwash  is  not  recommended  for  bending 
moment  predictions  for  tandem-rotor  configurations.  Use  of 
nonuniform  downwash  in  the  analysis  should  improve  airloads 
and  bending  moment  predictions  for  both  single-rotor  and 
tandem-rotcr  helicopters. 
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Figure  170.  CH-34  Predicted  and  Measured  Lift/Unit  Span 

at  0.75R  at  110  Knots 
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o  60  120  180  240  300 

ROTOR  AZIMUTH  (DEG) 

Figure  174.  CH-34  Predicted  and  Measured  Lift/Unit  Span 

at  0.55R  at  110  Knots 
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Figure  176.  CH-34  Predicted  and  Measured  Lift/Unit  Span 

at  0.85R  at  110  Knots 


179.  CH-34  Predicted  and  Measured  Lift/Unit  Span 

at  0.55R  at  150  Knots 
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Figure  181.  CH-34  Predicted  and  Measured  Lift/Unit  Span 

at  0.85R  at  150  Knots 


175  K 


e  187 


Figure  188.  CH-34  Predicted  and  Measured  Lift/Unit 

at  0.99R  at  175  Knots 


Part  V.  Impact  and  Potential 
of  New  Methods 


H.  R.  Alexander 
R.  B.  Taylor 
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IMPACT  OF  NEW  METHODOLOGY  ON  PREDICTION  CAPABILITY 


AEROFLASTIC  STABILITY 


The  stability  methodology  featured  in  Part  I  differs  from  that 
of  earlier  programs  in  several  important  aspects.  These  fall 
into  two  categories: 

(a)  The  new  capability  to  analyze  a  deflected  hingeless  rotor 
blade  fully  coupled  in  the  flap,  lag,  and  pitch  degrees 
of  freedom. 

(b)  An  upgraded  aerodynamic  representation  which  will  treat 
nonaxial  flow  through  the  rotor  without  the  usual  small 
angle . 

The  effectiveness  of  this  new  capability  has  been  demonstrated 
in  correlation  with  a  number  of  difficult  cases  which  had  not 
previously  been  correlated  satisfactorily.  Success  in  these 
areas  lends  credibility  to  design  predictions  which  will  be 
made  using  the  C-39  program  and  to  the  parametric  trend 
studies  already  accomplished  and  presented  in  Part  II  of  this 
volume . 

There  is  now  good  reason  to  believe  that  nonclassical  as  well 
as  classical  instability  mechanisms  can  be  analyzed  in  the 
design  phase  of  aircraft  featuring  large  rotors.  In  partic¬ 
ular  the  stability  of  wing-rotor  combinations  may  be  addressed 
in  all  phases  of  flight  from  hover  through  transition  to  high¬ 
speed  cruise. 

ROTOR  LOADS 


Results  of  the  rotor  loads  correlation  with  the  fully  coupled 
loads  analysis  clearly  showed  an  improvement  in  prediction 
capability  over  loads  analyses  used  in  the  past.  Correlation 
over  a  large  range  of  flight  conditions  with  a  number  of  rotor 
blades  whose  geometric  and  dynamic  characteristics  differ 
widely  demonstrated  that  coupling  of  the  pitch-flap-chord 
modes  and  the  use  of  unsteady  aerodynamics,  in  addition  to 
nonuniform  downwash  representation,  are  the  primary  sources 
of  improved  prediction  capability.  Specific  results  and  im¬ 
plications  of  the  correlation  as  applied  to  prop/rotor  and 
helicopter  loads  predictions  are  discussed  below. 

In  the  area  of  prop/rotor  loads,  the  coupled  analysis  showed 
the  capability  to  predict  blade  bending  moments  with  a  much 
greater  accuracy  than  previously  experienced  with  uncoupled 
analyses.  A  summary  of  the  prop/rotor  bending  moment  correla¬ 
tion  results  is  shown  in  Table  VI.  The  results  are  presented 
as  a  ratio  of  predicted-to-measured  bending  moment  sensitivity 
to  cyclic  pitch  rotor  angle  of  attack  or  rotor  rpm,  whichever 
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is  applicable  to  the  particular  model  and  flight  condition. 
First-harmonic  bending  moments  are  used  for  the  comparison 
since  prop/rotor  bending  moments  are  predominantly  1  per  rev. 
The  analysis  predicts  flap  bending  within  +5  percent  and  -11 
percent  and  chord  bending  within  +23  percent  and  -6  percent. 

The  quality  of  this  correlation  is  important  to  prop/rotor 
design  since  it  implies  that  the  coupled  analysis  can  be  used 
to  accurately  predict  fatigue  design  bending  loads.  Past  pre¬ 
dictions  of  prop/rotor  design  loads  have  been  a  mixture  of 
uncoupled  analysis  results  coupled  with  test  data  and  a  large 
amount  of  engineering  judgment.  Use  of  the  coupled  analysis 
for  prop/rotor  fatigue  design  loads  should  considerably  reduce 
the  amount  of  engineering  judgment  required. 

In  the  area  of  helicopter  loads  correlation,  the  analysis 
showed  the  importance  of  nonuniform  downwash  and  unsteady 
aerodynamics  on  loads  predictions.  Use  of  nonuniform  downwash 
increased  the  higher-harmonic  content  of  blade  flap  bending  due 
to  rotor  interference  and  improved  the  flap  bending  moment 
radial  distribution  predictions.  Use  of  nonuniform  downwash 
also  improved  chord  bending  correlation,  but  the  consistent 
underprediction  implies  poor  representation  of  the  lag  damper 
in  the  analysis.  The  chord  bending  correlation  results  show 
that  lag  damper  preload  is  not  properly  accounted  for  by  using 
an  analytical  torsional  viscous  damper  in  the  analysis.  Addi¬ 
tional  work  should  be  done  with  the  analysis  to  include  lag 
damper  preload  effects. 

The  highlight  of  the  helicopter  loads  correlation  was  the 
significant  improvement  of  pitch  link  load  predictions  for 
predominantly  unstalled  and  predominantly  stalled  conditions. 
Use  of  unsteady  aerodynamics  in  the  analysis  to  include 
hysteresis  effects  on  blade  section  lift  and  pitching  moment 
enabled  the  analysis  to  accurately  predict  pitch  link  load 
amplitude  at  the  onset  of  stall  and  the  growth  of  pitch  link 
loads  as  stall  increased.  Correlation  results  when  nonuniform 
downwash  was  used  in  the  analysis  were  erratic  for  predomin¬ 
antly  stalled  conditions,  indicating  strong  interference 
effects.  Since  pitch  link  load  predictions  were  good  both  in 
and  out  of  stall  when  uniform  downwash  was  used  in  the  analy¬ 
sis,  use  of  uniform  downwash  is  recommended  for  pitch  link 
load  predictions.  The  quality  of  the  pitch  link  load  correla¬ 
tion  indicates  that  the  analysis  can  be  used  to  determine 
fatigue  design  pitch  link  loads  for  stalled  and  unstalled 
conditions  for  helicopter  blade  design. 
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TABLE  VI 

SUMMARY  OF  PROP/ROTOR  BENDING  MOMENT  CORRELATION  RESULTS 
RATIO  OF  PREDICTED  TO  MEASURED  BENDING  MOMENT  SENSITIVITY  TO  A  FOR 
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POTENTIAL  EFFECTS  OF  NEW  CAPABILITY  ON 
DESIGN  PRACTICE  AND  CRITERIA 


AEROELASTIC  STABILITY 

The  impact  of  the  new  methodology  on  design  practices  may  be 
expected  to  arise  from  increasing  confidence  in  our  capability 
to  predict  and  prevent  instabilities  that  had  previously  been 
beyond  the  range  of  routine  analysis.  Thus  rotor-airframe 
designs  of  advanced  configurations  may  now  be  investigated 
with  confidence.  This  may  be  expected  to  lead  to  a  gradual 
relaxation  of  conservative  conventions  currently  applied  in 
blade  design  since  novel  ideas  may  be  investigated  paramet¬ 
rically  to  select  those  which  are  sufficiently  promising  to 
initiate  hardware  programs  or  test  investigations.  Currently, 
certain  designs  which  may  appear  attractive  peripherally  will 
be  aborted  at  the  analysis  stage  without  the  expense  of  exten¬ 
sive  test  programs. 

In  pa  -icular,  the  new  analysis  provides  a  tool  for  investi¬ 
gating  and  understanding  pitch-lag-flap  flutter  which  is  cur¬ 
rently  a  matter  of  concern  in  hingeless  rotor  designs.  It  may 
be  expected  that,  with  additional  understanding,  design  cri¬ 
teria  and  procedures  will  be  developed  which  will  reduce  this 
type  of  mechanism  to  the  same  status  as  other  forms  of  aero- 
elastic  behavior. 

ROTOR  LOADS 


The  results  of  the  loads  correlation  indicated  two  areas  of 
significant  improvement  of  rotor  loads  predictions  for  which 
the  analysis  can  be  used  to  develop  design  criteria.  In  the 
area  of  prop/rotor  design,  correlation  of  alternating  flap 
and  chord  bending  moments  and  steady  rotor  moments  in  hover 
due  to  cyclic  pitch  showed  significant  improvement  over  pre¬ 
dictions  made  in  the  past  with  uncoupled  analyses.  Since 
hover  cyclic  produces  fatigue  design  blade  loads  for  prop/ 
rotor  design,  the  analysis  can  be  used  as  a  tool  to  develop 
an  improved  prop/rotor  design  criteria. 

In  the  area  of  helicopter  loads,  the  good  correlation  of  pitch 
link  loads  indicates  that  the  analysis  can  be  used  to  develop 
a  sounder  pitch  link  design  criterion  for  both  unstalled  and 
stalled  conditions.  The  analysis  can  be  used  in  conjunction 
with,  or  possibly  replace,  pitch  link  load  flight  test  data 
for  pitch  link  design.  As  the  amount  of  correlation  of  flight 
test  data  with  the  analysis  increases,  confidence  in  the  pitch 
link  load  prediction  capability  of  the  analysis  should  also 
increase,  and  the  analysis  can  then  be  used  to  replace  the 
extrapolation  of  existing  flight  test  data  for  pitch  link  load 
design  of  heavier  and  faster  helicopter  designs. 
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